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COMPILER AND OBJECT TIME SYSTEM 



PART 7 
CHAPTER 1 
INTRODUCTION TO FORTRAN 

FORTRAN (FORmula TRANslation) is a problem-oriented language designed to permit 
scientists and engineers to express mathematical operations in a form with which 
they are familiar. The term FORTRAN is used interchangeably to designate both 
the FORTRAN "language and the FORTRAN Compiler or translator. A program written in 
the FORTRAN language is called a FORTRAN source program. 

A FORTRAN source program is composed of commands describing the functions per- 
formed and computational elements expressed in a notation similar to that of 
mathematics. The source program is compiled by the FORTRAN Compiler into code 
which is then automatically assembled into a binary object module. The Linker 
program links the resultant relocatable module with any other required module or 
modules to create an executable program image called a load module. Execution is 
started by using the appropriate commands as described in Chapter 7-9 of this 
manual . 

There is no one-to-one correspondence between a FORTRAN statement and a machine- 
language instruction. Some FORTRAN statements generate many machine instructions 
while others yield none. 

1.1 DOS/BATCH FORTRAN 

DOS/BATCH FORTRAN conforms to the specifications for American National Standard 
FORTRAN. It is also highly compatible with IBM 1130 FORTRAN. The following 
additions to American National Standard FORTRAN are included in DOS/BATCH FORTRAN: 

a. An array subscript may consist of any valid arithmetic expression. If 
the result of the expression is not an integer, it is truncated and 
converted to integer format before being used as a subscript. 

b. Hexadecimal, octal and Radix-SjZi' constants are allowed within DATA 
statements. 

c. Hollerith constants may "be delimited by single quote (apostrophe) 
characters. 

d. General mixed-mode expressions are allowed for all data types, 
including complex. 

e. . Assigned GOTO statements may contain optional label lists.' 
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f. The following statements have been added to handle random access I/O; 



DEFINE FILE 

FIND 

READ (a'b) 

WRITE (a'b) 



g. The IMPLICIT statement is provided to allow redefinition of the rules 
for data type classification of variables by default. (In tlie absence 
of any data type specification/ a variable name beginning with one of 
the letters I through N indicates an integer variable; A through H and 
through Z indicate real variables.) 

h. The data types 

LOGICAL* 1 
INTEGER* 2 
REAL* 4 
REAL* 8 

are provided. 

j . The I/O conditions END=n and ERR=n can be specified on any READ or WRITE 
statement to transfer control upon detection of an end-of-file or error 
condition. 

k. Default values for field width and number of characters to the right of 
the decimal point are provided on numeric and logical format conversions. 

1. Variable format expressions are provided. 

m. A simplified type of free form formatted input (termed "short field 
termination") is provided. 



1.2 SOFTWARE AND HARDWARE ENVIRONMENTS 

The FORTRAN Compiler produces code suitable for assembly by MACRO. The user can 
compile and assemble his FORTRAN program in one step (which* is normally done) , or he 
can compile the program into assembly language for later assembly (indicated by 
means of a switch option to the compiler). FORTRAN programs can call assembly 
language subroutines. The FORTRAN Library routines can be called by either FORTRAN 
or assembly language code. 

DOS /BATCH FORTRAN programs can be compiled and run on any hardware configuration 
that supports the DOS /BATCH operating system, and that has a minimum of 16K of 
memory. DOS/BATCH FORTRAN supports all standard hardware options supported by the 
operating system. 
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1.3 HOW TO USE THIS PART 

This Part is primarily a reference document describing the DOS /BATCH FORTRAN 
language. It is not intended to teach the reader how to program in FORTRAN. The 
bulk of this Part describes details of the language statements and descriptions 
of the file handling, library, internal structures, and Object Time System. Later 
chapters provide concise summaries of material in tabular form. 

Although designed as a reference tool, this Part contains a number of programming 
examples to show the use of statements in a context. Internal details supplied 
will aid the more advanced FORTRAN programmer. 

1.4 FORTRAN STATEMENT STRUCTURE 

A statement is the basic functional unit of the FORTRAN language. An executable 
statement performs a calculation, does I/O, or directs control of the program; a 
nonexecutable statement provides the Compiler with information regarding variable 
structure, array allocation, storage requirements, etc. 

All categories of FORTRAN statements are coded in a standard format for ease of 
program preparation. Figure 7-1 shows a FORTRAN coding form (there are many types 
available) with four miscellaneous lines written on it. FORTRAN coding forms lend 
themselves to paper tape input as well as to punched card input preparation. The 
form is divided into four major sections: 

a. Statement number field - columns 1-5. 

b. Line continuation field - column 6. 

c. Statement field - column 7-72. 

d. Sequence field (ignored by the FORTRAN Compiler) - columns 73-80. 

,1.4.1 Formatting a FORTRAN Line 

Punched card input for the DOS/BATCH FORTRAN Compiler is formatted as shown on 
the coding form in Figure 7-1. This is the standard format for FORTRAN punched 
cards on most computer systems. 
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o 
fa 

c 

•H 
TJ 
O 
U 



When programs are prepared on-line with the Editor, statements can be typed in the 
same manner as described for cards or they can use the TAB character as follows: 



A TAB causes the internal character count to be set to 6 or 7 
depending on whether or not the character following the TAB is numeric. 
If the character is not numeric, it is assumed to be the start of a 
statement proper. If the character following the tab is numeric, the 
line is assumed to be a continuation line. 

A 1- to 5-digit statement label may precede the TAB character. 



In order to obtain the formatting shown on the coding form in Figure 7-1, the 
user would type the following to EDIT; 



123-»jX=12.+3. 


or 


-+|a=b+c 


or 


HlAl'ABCy 


or 


H COMMENT LINE 


or 



123AAAx=12.+3 

AAAAAAA+B+C 

AAAAAlAl'ABC ' (continuation line) 

CAAAAACOMMENT LINE 



where : 



H 



indicates typing of the TAB character (CTRL/TAB) 
indicates typing of the SPACE bar. 



As specified in Section 7-2.1, the TAB character is not included in the legal FORTRAN 
character set; therefore, the TAB character can be used only at the beginning of a 
line to indicate spacing to the continuation or statement field. FORTRAN ignores 
spaces within each field (except within Hollerith constants) ; spaces may therefore 
be inserted at the user's convenience to increase program legibility. 

1.4.2 Statement Number Field 

A statement number consists of one to five digits in columns 1 to 5. Leading zeros 
and all blanks in this field are ignored. Any statement to which reference is made 
by another statement must have a statement number. A label consisting only of one 
or more zero characters is not allowed. 

Statement numbers must be unique. If two statements within a program unit have the 
same statement number, the error is noted at assembly time. There is no limit to 
the value of the statement number. 
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1.4.3 Continuation Lines 

Any FORTRAN statement may be continued onto more than one physical source line. A 
nonzero or nonblank character is placed in column 6 of all continuation lines (or 
cards). When using an editor, the TAB character followed by a nonzero numeric 
character may be used to indicate a continuation line. 

DOS/BATCH FORTRAN normally allows up to five continuation lines; however, this 
number may be changed to any number in the range to 99 through the use of the 
/CO switch (see Appendix J) . 

1.4.4 Comment Lines 

A comment may be inserted into a FORTRAN program by making the first character of 
the line a C. If the C character is detected in column 1 of a card or as the first 
character of any input record, all information on that line is disregarded. A 
comment line must not precede a continuation line. Continuation lines cannot be 
used to continue comments. All comment lines begin with the letter C. 

1.5 PROGRAM UNIT STRUCTURE 

# 

A FORTRAN program unit is a FORTRAN main program or a subprogram. Except for BLOCK 
DATA subprograms, which are a special class (described in Section 7-6.4.8), program 
units are constructed of FORTRAN statements in the following order: 

a. Subprogram statement, if the program unit is a subprogram. 

b. IMPLICIT statement, if used. 

c. Specification statements (other than IMPLICIT). Those specification 
statements that determine the storage required by a variable or 
array must precede statements that initialize the values of the 
variable or array. DATA statements must follow all other specification 
statements. 

d. Arithmetic statement function definitions, if used. 

e. Executable program statements. Each program unit (except BLOCK DATA 
subprograms) must have at least one executable statement. 

f. END statement. 
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PART 7 

CHAPTER 2 

FORTRAN STATEMENT COMPONENTS 



2.1 FORTRAN CHARACTER SET 



The character set from which statements can be constructed is as follows. 



the letters A through Z 
the digits fi through 9 

blank 
= equal sign 
+ plus sign 

minus sign 
* asterisk 
$ dollar sign 



/ slash 

( left parenthesis 
) right parenthesis 
, comma 

decimal point 



The TAB character can be used only at the beginning of a line to indicate spacing 
to the continuation or statement field. Other printing characters (except +) 
may appear within a FORTRAN program only as part of a comment line or a Hollerith 
constant (text string) . 

While the PDP-11 Monitor provides device independence for input to the FORTRAN 
Compiler, it is mandatory that FORTRAN source programs be in ASCII code. Programs 
intended for compilation by the DOS /BATCH FORTRAN Compiler and prepared on cards, 
paper tape, or magtape from other systems must conform to this requirement. 

2.2 FORTRAN ' CONSTANTS 



A constant is a self-defining value that does not change from one execution of 
the program to the next. In general, constants can be numeric, logical, or 
character string entities. 

Types of constants allowed in DOS/BATCH FORTRAN include: 

integer \ 

real I 

double precision \ s 

octal 1 

hexadecimal / 
complex 



numeric constants 
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logical - logical constants 

Hollerith - character string constants 

Storage formats for each type of constant are described in Chapter 7-12. 

2.2.1 Integer Constants 1 

An integer constant is a string of numeric characters, from one to five digits in 
length; it is interpreted as a decimal integer. Integers are 16-bit signed quantities. 
A negative integer is preceded by a minus sign; a positive integer is optionally pre- 
ceded by a plus sign. An integer constant must lie within the range -32768 to 
+32767 (which is -2 15 to.2 15 -1). 

If the result of any operation exceeds 16 bits, the Object Time System displays an 
error diagnostic. 

Examples of integers acceptable to DOS/BATCH FORTRAN follow. 

3564 

+14 

-26357 

Commas and decimal points are not allowed in integer constants. 

2.2.2 Real Constants 

A real (or floating-point) constant is a string of numeric characters including a 
decimal point. A real constant may contain any number of digits, but only the 
leftmost eight digits are significant. A negative value is preceded by a minus sign; 
for positive values, the plus sign is optional. 

A real constant may alternately be expressed in scientific notation with a decimal 
exponent, represented by the letter E, and a signed integer constant. The decimal 
point with the real number is optional if there is an exponent. The exponent field 
cannot be blank, but may be Eero. For positive exponents, the plus sign is 
optional . 

A real constant has 24 bits of precision, or somewhat more than seven decimal digits. 
The magnitude must lie within the range tf. 14*ljZf 38 to 1.7*ljzf 38 . Real constants occupy 
two words of PDP-11 storage (see Chapter 7-12.) 



1 See Appendix J for the effect of the /ON switch on the storage allocation of integer 
values . 
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Examples of real constants follow. 



1J . 

.579 






* 0.0 






-ljtf.685 






5.0E3 


(i.e. , 


5000,) 


4.2E+3 


(i.e. , 


4200.) 


2.J0E-3 


(i.e. , 


.002) 


5E0 


(i.e. , 


5.0) 



2.2.3 Double-Precision Constants 

A double-precision constant is a string, of numeric characters including a decimal 
point and the letter D followed by an integer exponent. The decimal point is 
optional. A double-precision constant may contain any number of digits, but only 
the leftmost 16 digits are significant. A double-precision constant is always 
expressed in scientific notation with a decimal exponent (represented by the letter 
D and a signed integer constant) . The exponent field cannot be blank, but may be 
zero. The plus sign preceding a positive constant and/or exponent may be omitted. 

A double-precision constant has 56 bits of precision, or about 16 decimal digits. 

— 38 
The magnitude of a double-precision constant must lie between .14*10 and 

38 
1.7*ljd • Double-precision constants occupy four words of PDP-11 storage. 

Examples of double-precision constants are shown below. 

24.6713264837465D0 

2.5D2 (i.e., 25,0. to 16 digits precision) 

3.6D-2 (i.e., 036 to 16 digits precision) 

5.6D0 (i.e., 5.6 to 16 digits precision) 

2.2.4 Octal Constants 

The general form for an octal constant is a string of one to six octal digits (0 to 
7 inclusive) preceded by a quotation mark (") . For example: 

"120 

"0 

"177777 

Octal constants of this form may appear anywhere an integer constant is valid. An 
octal constant may also be a string of one to six octal digits preceded by the letter 
0. Octal constants of this form are valid only in DATA, STOP and PAUSE statements. 
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Examples are shown below. 



0120 

00 

0177777 



The maximum value of an octal constant is 177777. 



2.2.5 Hexadecimal Constants 

A hexadecimal constant is a string of one to four hexadecimal digits (0 to 9 r and 
A to F, inclusive) preceded by the letter Z. The hexadecimal digits follow. 

Hexadecimal Decimal 

Digit Equivalent 



1 1 

2 2 

3 3 

4 4 

5 5 

6 6 

7 7 

8 8 

9 9 
A lp 
B 11 
C 12 
D 13 
E 14 
F 15 

A hexadecimal constant is valid only in a DATA statement. The maximum value of a 
hexadecimal constant is FFFF. Some examples: 

ZA34 

ZC1 

ZFFFF 



2.2.6 Complex Constants 

A complex constant is written as an ordered pair of real constants separated by a 

comma and enclosed in parentheses. The first constant represents the real part of 

the complex number; the second represents the imaginary part. Each may be signed. 

The parentheses are part of the constant and must appear. A complex constant is 

represented internally by four consecutive words of PDP-11 storage. Examples 

follow. 

(.7jef712,-K 70712) 
(8.763E3,2.297) 
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2.2.7 Logical Constants 1 

The two logical constants .TRUE, and .FALSE, are represented internally by the 
integer values -1 and 0, respectively. The delimiting periods are part of the 
constant and always appear. Any nonzero value is considered true when tested by a 
logical IF statement; -1 is equivalent to the octal number 177777. These values 
may be entered via DATA statements as .TRUE, and .FALSE. (See Section 7-5.7). (See 
Section 7-7.1/5 concerning input and output of logical values.) 

Both arithmetic and logical operators may be used with logical constants, which are 
treated as one-word integer values. 

2.2.8 Hollerith Constants 

A Hollerith (or literal) constant is a string of up to 255 characters, which may be 
represented in one of two ways. 

a. An integer constant indicating the number of characters, 
followed by an H, followed by a character string: 

5HW0RDS 
4H1#23 

b. A character string enclosed in single quote characters: 

'WORDS 1 
, 1#23' 

The single .quote (apostrophe) character may be included in the character 
string if it is immediately preceded by another single quote character. 
For example: 

, DQN ,, T' 

Hollerith constants may be entered in DATA or input statements as a string of one 
or two ASCII characters per integer variable, one to four characters per real 
variable, and one to eight characters per complex or double-precision variable. 
Hollerith constants are stored in memory as byte strings, filling up to word 
boundaries. If a Hollerith constant has an odd number of characters, a blank is 
appended to the end of the constant. 



l See Appendix J for the effect of the /ON switch on the storage allocation of 
logical values. 
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A Hollerith constant is treated as an integer quantity when it appears in an 
arithmetic context. For example: 

I = 'ABC +34 
is equivalent to 

I = 'AB' +34 

The third and any subsequent characters are ignored. The integer value of 'AB' 
can be determined by examining the storage format, which is shown here. 



15 







ASCII value 
of B 



ASCII value 
of A 



2.2.9 Radix-Sjzf Constants 

Radix-5J3 is a special character data representation in which up to three characters 
from a limited set of characters may be packed into a single PDP-11 word. Radix-5^ 
constants are allowed only in DATA statements. They are represented by a leading 
R followed by up to six characters and may be used to initialize only variables 
of type REAL. 



A Radix-50 constant terminates with the first non-Radix-50 character or the sixth 
Radix-5^ character, whichever comes first. Examples are shown below. 



R123456 
RJftStftftfA 
Rl.$ 



If i, j, and k are three characters from the Radix-5jtf set, they are combined into a 
single word value using the formula: 



(<i * 50e + j) * 508 + k) 



The Radix-50 characters and their code values follow. 



Character 

Space 
A - Z 



Radix- 5 $ Value 
(Octal) 

1-32 
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Radix-50 Value 
Character (Octal) 

$ 33 

34 ' 

(not used) 35 

0-9 36-47 



Note that leading and embedded blanks are significant if they are within the 
6-character limit, since blank is a valid Radix-5jtf character. Since the Radix-50 
code value blank is zero, trailing blanks may be omitted without affecting the 
Radix- 50 value. 

2.3 FORTRAN VARIABLES 

A simple FORTRAN variable is a named quantity whose value may change during execu- 
tion of a program, a variable name is a string of one to six alphanumeric characters, 
the first of which must be alphabetic. Variable names longer than six characters 
are rejected by the FORTRAN Compiler, and cause an error diagnostic to be printed. 

Legal Names Illegal Names 

ALPHA 2A (begins with a number) 

MAX MAXIMUM (more than 6 characters) 

A34 



A variable has a principal attribute type, just as constants are of one type or 
another. The type of variable indicates the type of value it can be assigned, its 
precision, and its storage requirements. 

Type INTEGER, REAL, LOGICAL, DOUBLE PRECISION or COMPLEX is assigned by means of an 
explicit type declaration statement (see Section 7-5.2) or by an IMPLICIT statement 

(see Section 7-5.1). INTEGER or REAL may also be assigned by virtue of the initial 

letter of the variable name. 

In the absence of any type declarations or IMPLICIT statements, a variable name that 
begins with I, J, K, L, M, or N is assumed to represent an integer variable. All 
other initial letters indicate a real variable. 

Integer' Variables Real Variables 

MAX REAL 

IDAT DATA1 

NJOB X 
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2.3.1 Integer Variables 

Integer variables undergo arithmetic calculations with automatic truncation of any 
fractional part. For example, if the current value of K is 5 and the current value 
of J is 9, J/K yields 1 as a result. Integer values can be converted to real values 
by the FLOAT function or by an arithmetic assignment statement. Integer arithmetic 
operations do check for overflow and print an error message if the result of an 
integer operation overflows 16 bits. 

2.3.2 LOGICAL* 1 Variables 

Data items in L0GICAL*1 format are eight bits long (1 byte) , and are always accessed 
by means of a variable name defined as BYTE or LOGICAL* 1 type. The range of numbers 
from +127 to -128 can be represented. Any arithmetic operation is performed internally 
by taking two 8-bit operands, extending the sign to a full word, and performing the 
desired operation. No diagnostics are issued for byte overflow; the result is trun- 
cated to an 8-bit quantity. 

Logical and masking operations are performed with one byte of data at a time. 
L0GICAL*1 (or BYTE) array elements are stored in adjacent bytes. 

2.3.3 Array Variables 

The extent of a variable refers to the number of values that are represented by a 
single variable name. A scalar variable represents a single quantity; e.g., BETA 
represents a single numeric value. 

An array variable represents a one, two, or three dimensional array of quantities. 
A given variable name cannot be used as both a simple variable and an array name 
in the same program unit. An array element is denoted by the array name followed 
by a subscript list enclosed in parentheses. The subscript list can be any one, two, 
or three arithmetic expression (s) separated by commas. 

An array can be defined with a DIMENSION or COMMON statement or as part of a type 
declaration statement. Array dimensioning is discussed in Section 7-5.3. 

Any legal arithmetic expression can be used as a subscript of any array element 
once the array has been defined. Noninteger subscripts are converted to integers 
by truncation of the fractional part before use. 

The extent of any array variable is determined by the number of dimensions it is 
assigned. For example, the following are legal references to array elements. 
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one dimensional array (Al) 

two dimensional array REL(24,3) 

three dimensional array TR4 (X*2,5,2*X) 

Arrays are stored in contiguous storage locations that are addressed in ascending 
order with the first subscript varying most rapidly. For instance, the 2-dimensional 
array B(J,K) is stored in the following order: 

B(1,1),B(2,1),...,B(J,1) 
B(l,2) ,B(2,2),...,B(J,2)' 



B(1,K), ... ,B(J,K) 

See Section 7-5.3.1 for a detailed description of array storage. 

References to an array variable must contain the number of subscripts with which 
the array was defined, the exception being those statements that can refer to an 
entire array given only the array name. 

2.4 FORTRAN EXPRESSIONS 

FORTRAN expressions can be arithmetic, relational, or logical. Relational or logical 
expressions have logical values that can be treated as integers. All expressions 
are combinations of elements and operators. The different types of expressions are 
discussed below. 

2.4.1 Arithmetic Expressions 

An arithmetic expression is a combination of arithmetic operators and elements. The 
arithmetic operators follow in descending order of priority. 

Operator Meaning 

** Exponentiation (A** 2) 

+,- unary operators (algebraic signs: +3,-1) 

*,/ multiplication (A*B) , division (1/2) 

+,- addition (A+B) , subtraction (D-4) 
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An arithmetic expression cannot contain adjacent arithmetic operators. An element 
preceded by a unary operator must therefore be separated from a preceding binary 
operator by parentheses, as in A + (-B) . 

In cases where two operations within an expression have the same priority (*,/ and 
+,-) the operations are performed from left to right. For example: 

-A**2+B*C-1 

is performed as 

(-(A**2)) + (B*C) -1 

Additional computations (such as sine, cosine, or square root) can be specified 

via library function references. A function reference acts as a basic element in 

an expression since all functions return a single value. (For example, SQRT(4.0) 
returns the value 2.0 in an expression.) 

An arithmetic expression may consist of a single element: 

B4 

22.7 

NEX(3) 

or a combination of arithmetic operators and basic elements. 

-X+3 
A/2+DAT 
TAN (PI*M) 
B**(-2) 

2.4.2 Use of Parentheses 

Any arithmetic expression can be enclosed in parentheses and considered as a basic 
element (it is evaluated before the remainder of the expression). Parentheses can 
be used to change the order (and the result) of evaluation. For example: 

4+3*2-6/2 = 4+ (3*2) -(6/2) = 4+6-3 = 7 
(4+3)*2-6/2 = (7*2) -(6/2) = 14-3 =11 
(4+3*2-6) /2 « (4+ (3*2) -6) /2 * (4+6-6) /2 - 4/2 - 2 
((4+3)*2-6)/2 - ((7*2)-6)/2 = (14-6) /2 =8/2=4 
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Within parentheses, the order of operator evaluation is as follows: 

a. function evaluation 

b. exponentiation 

c. unary minus 

d. multiplication, division 

e. addition, subtraction 

Parentheses cannot be used to imply multiplication. That is: 

A(2) 
is the second element of the array A, whereas 

A*2 . 

indicates twice the value of the variable A. 

2.4.3 Mixed Mode Arithmetic Expressions 

With few exceptions, any type of constant, variable or function can be combined 
with any other type in an arithmetic expression. The type of the resultant 
expression when any two types are combined is found in Table 7-1. 

Logical, octal and Hollerith constants are treated as integer constants when com- 
bined with other elements in arithmetic expressions. The result of such a combination 
is then changed to the appropriate resultant mode. 

Arithmetic operations can be performed between complex and other data types. Real, 
integer, and double-precision values (represented by V) are converted, if necessary, 
to the complex value (V,0.0) prior to the arithmetic operations; V is a real value 
in this representation. 
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Table 7-1 
Mixed Mode Arithmetic Results 







Byte 
Logical*l 


Integer 


Logical 


Real 


Double 
Precision 


Complex 


















Byte 
Logical*l 

Integer 
Logical 

Real 

Double 
Precision 

Complex 




Byte 

Integer 
Logical 

Real 

Double 
Pre- 
cision 

Complex 


Integer 

Integer 
Integer 

Real 

Double 
Pre- 
cision 

Complex 


Logical 

Integer 
Logical 

Real 

Double 
Pre- 
cision 

Complex 


Real 

Real 
Real 

Real 

Double 
Pre- 
cision 

Complex 


Double 
Precision 

Double 
Precision 

Double 
Precision 
Double 
Precision 

Double 
Pre- 
cision 

Complex 


Complex 

Complex 
Complex 

Complex 

Complex 
Complex 


The general rule for mixed modes is that the result has the type of the most sig- 
nificant operand. 


The following special rules apply for determining the type resulting from expressions 
of the form A**B: 

1. If B is integer, the expression is of the same type as A; 

2. If A and B are both Real, the expression is real; 

3. If A or B is Double-precision, the expression is Double-precision; 

4. No other possibilities for A and B are legal. 


When a real number is converted to Double-precision format, repeating decimals are 
not carried through the additional precision. 
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The only illegal mixed mode operations are 



B**X 
I**R 
I**D 
C**R 
C**D 
X**B 
X**G 



where 

B is a Byte value. 

I is an Integer value. 

R is a Real value. 

C is a Complex value. 

D is a Double-precision value. 

X is a value of any data type. 

2.4.4 Logical Expressions 

A logical expression combines constants, logical variables, logical function 
references and arithmetic expressions using relational or logical operators. 
Logical masks can be represented by using octal constants. The logical constants 
.TRUE, and .FALSE, can be used in any logical expression and have the values -1 
and 0, respectively. The result of a logical expression is the logical value true or 
false and uses one word of storage space. (Any nonzero integer value is considered 
true when tested by a logical IF statement. ) 

A relational expression consists of two arithmetic expressions separated by a 
relational operator. The arithmetic expressions can be of any complexity and any 
combination of integer, real, or double-precision types. The relational operators 
follow (where A and B are arithmetic expressions) . 

Relational 
Operator Example Meaning 

.GT. A.GT.B A is greater than B 

.GE. A.GE.B A is greater than or equal to B 

.LT. A.LT.B A is less than B 

.LE. A.LE.B A is less than or equal to B 

.EQ. A.EQ.B A is equal to B 

.NE. A.NE.B A is not equal to B 
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Relational expressions can be written comparing values of different precision 
(having different data types) . All lower-precision values are converted to the 
highest precision format. The following should be noted: 

a. When combined with, or related to, a complex value, a real, integer, 

or double-precision value (represented by V) is converted to the complex 
value (V,0.0) . V is a real value in this representation. 

b. When a real number is converted to double-precision format, repeating 
decimals are not carried through the additional precision. This causes 
a real 0.333..., for instance, to be less than a double-precision 
0.333... . 

Only the relational operators .EQ. and ,NE. can be used with complex expressions 
(complex quantities are equal if their corresponding parts are equal) . 

The value of a relational expression is always either true or false (-1 or 0,*no 
other value) . The following are relational expressions. 

A.GT.B 

A* * 2+C . LT . NET *B+4 

AB/C.LE.AB/D 

Parentheses can be used in relational expressions. For example: 

(A+B)*C.EQ. (D+l)**2 
(A(1)+1).GT.(B(1)) 

A logical expression may also contain logical operators which operate on relational 
expressions. The logical operators are as follows. 



Logical 
Operator 



Example 



Meaning 



.NOT. 



.AND. 



.OR. 



.NOT. A 



A.AND.B 



A.OR.B 



Has the value .TRUE, only if A is 
.FALSE.; has the value .FALSE, if A 
is .TRUE. 

Has the value .TRUE, only if A and B 
are both .TRUE.; has the value .FALSE, 
if either A or B is .FALSE. 

Has the value .TRUE, if either A or B 
or both are .TRUE.; is .FALSE, only 
if both A and B are .FALSE, (inclusive 
OR) . 
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A logical expression may consist of basic logical elements or a combination of 
logical elements. For example, the following are all logical expressions (relational 
expressions are a subtype of logical expressions) . 

.TRUE. 

X.GE.3.14159 
TV (1) .AND. INDEX 
B(M) .OR. (K.EQ. LIMIT) 

No two relational or logical operators may appear in sequence, except in the case 
where .NOT. appears as the second of two logical operators. Any logical expression 
can be preceded by the unary operator .NOT. as follows. 

.NOT.T 

.NOT. (X+7.GT.Y+Z) 

B(M) .AND.. NOT. TV. OR. R which is the same as (B(M) .AND. (.NOT. TV) ) .OR.R 

However, where .NOT. appears next to an arithmetic operator, the arithmetic expres- 
sion must appear in parentheses. For example? 

.NOT. (+1.0) .LE.DA(2) 

as opposed to the following form, which is incorrect: 

.N0T.+l.jd.LE.DA(2) 

Logical and relational operators have priorities in the following order when the 
precedence is not established by parentheses. 

.GT. .GE. .LT. .LE. .EQ. .NE. .NOT. .AND. .OR. 

Relational operators are evaluated in an expression from left to right subject 
to the priority rules above. In a logical expression, all arithmetic evaluation 
is performed prior to any logical evaluation. Thus, the logical expression 

.N0T.Z**2+Y*MA .GT. K-2 .OR. PA .AND. X .EQ. Y 

is interpreted as * 

(.N0T.((CZ**2) + (Y*MA)) .GT.(K-2))) .OR. (PA.AND. (X.EQ.Y) ) 
I l ~ ' . 
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Parentheses may be used within logical expressions to make an expression more 
readable or to override established priorities. For example, the following two 
expressions are evaluated differently. 

A. AND. (B.OR.C) 
(A.AND.B) .OR.C 

2.5 OPERATOR SUMMARY 

Operators in each type are shown in order of descending priority. 



T yp e 



Operator 



Operates Upon 



arithmetic 



** 

+ ,- 
*,/ 

+ ,- 



exponentiation 
unary plus, minus 
multiplication, division 
addition, subtraction 



numeric constants 
variables and 
expressions 



relational 



logical 



• GT. 
.GE. 
.LT. 
-LE. 
.EQ. 
.NE. 

.NOT. 



greater than 

greater than or equal to 

less than 

less than or equal to 

equal to 

not equal to 

•NOT. A is true if 
and only if A is 
false; it is false 
if A is true. 



logical variables, 
logical constants 
and arithmetic 
expressions (all 
relational operators 
have equal priority) 

logical variables 
and logical constants 



.AND. 



.OR. 



A.AND.B is true if and 
only if A and B are both 
true; it is false if either 
A or B is false. 

A.OR.B is true if 
either A or B or both 
are true; it is false 
if both A and B are 
false. 
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PART 7 

CHAPTER 3 

ASSIGNMENT STATEMENTS 

There are three types of assignment statements. 

1. Arithmetic assignment statements 

2. Logical assignment statements, and 

3. ASSIGN statements. 

A value can be assigned to a variable at any point in a source program. During 
program execution, the most recent assignment determines the value of that variable 
in subsequent statements . 

3.1 ARITHMETIC ASSIGNMENT STATEMENT 

The arithmetic assignment statement provides the means by which the results of 
computations can be stored. 

The format of the arithmetic assignment statement is 

variable = expression 

where variable is any legal simple or subscripted variable name, expression is any 
legal arithmetic expression, and = is the replacement operator. The equal sign does 
not indicate equality as in an algebraic statement, but rather a replacement. The 
value of the expression is assigned to the variable name, regardless of any previous 
value of the variable. For example, 

I ■ 1+1 

takes the current value if I; adds 1 to it, and uses the new value as the current 
value of I for any subsequent operations. 

The expression to the right of the equal sign is evaluated and converted when 
necessary to conform to the type of the variable to the left of the equal sign. 
That is, if a real expression is assigned to an integer variable, the value of the 
expression is converted (truncated) to an integer before assignment. For example: 

AND = Y*(X**2+Z) 
I = I*N 
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X(J) = (J)-B(J) 
J - A**2+C 



The expression to be assigned must be capable of yielding a value that conforms to 
the type attribute of the variable that is to receive that value. Type conversions 
are performed in accordance with the rules stated for mixed mode expressions in 
Section 7-2.4.3 and Tables 7-1 and 7-2. 



Table 7-2 
Conversions Rules for Assignment Statements 
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3.2 LOGICAL ASSIGNMENT STATEMENT 

A logical assignment statement is of the form 

variable = expression 

where variable is a predefined logical variable name or a logical array element 
name, expression is a logical expression, and = is the replacement operator. As 
in the arithmetic assignment statement, the equal sign indicates that the value of 
the expression (true or false) is given to the variable specified. For example: 

P ■ .TRUE. 

VAL = X.GT.4.0R.B.EQ.X 

A - R.LE.X**2 

In the case where mixed mode arithmetic expressions are being logically compared, 
Table 7-3 shows the mode in which the comparison is effected. For example, 

A - x**2 .GE. I 

causes X**2 to be evaluated, the value of the integer variable I converted to real 
format, and the comparison made. 

3.3 ASSIGN STATEMENT 

The ASSIGN statement is of the form 

ASSIGN n TO ivar 

where n is a statement number and ivar is an integer variable. The ASSIGN statement 
permits a subsequently-executed assigned GOTO statement (see Section 4.1.3) to 
transfer control to the statement identified by the variable ivar. The ASSIGN 
statement and the associated GOTO statement (s) must appear in the same FORTRAN 
program unit. The statement number assigned must be that of an executable state- 
ment (assignment, I/O, or control statement) . 

An integer variable used in an ASSIGN statement must be redefined in an arithmetic 
assignment statement before it can be used in any context other than the GOTO 
statement (and after which it can no longer be used in the GOTO statement) . For 
example, 

ASSIGN 10 TO COUNT 
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associates the variable name COUNT with statement number 10, after which the state- 
ment 

C0UNT=C0UNT+1 

is considered invalid by the FORTRAN system. The above statement can be performed, 
however, if preceded by the statement. 

COUNT=2jtf 

which assigns the value of 20 to the variable COUNT. But the variable COUNT can no 
longer be used to refer to statement number 10. Incorrect usage of a variable 
defined in an ASSIGN statement causes a warning diagnostic to be printed by the 
Compiler. 
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PART 7 

CHAPTER 4 

CONTROL STATEMENTS 

4.1 GOTO STATEMENTS 

GOTO statements allow program execution to jump from one section of code to another, 
either directly or based on the value of a variable. 

There are three types of GOTO statements. 

a. Unconditional GOTO statements 

b. Computed GOTO statements 

c. Assigned GOTO statements 

4.1.1 Unconditional GOTO Statement 

The unconditional GOTO statement transfers control directly to the specified 
statement number. The form of this statement is 

GOTO n 

where n is. the statement number of an executable statement. For example: 

GOTO 10)3 

When the above statement is encountered, program execution continues at the state- 
ment identified by statement number 100. 

An unconditional GOTO statement may appear anywhere in the executable part of a 
source program, except as the terminal statement of a DO Loop. The GOTO statement 
and the statement to which it transfers control must both appear in the same FORTRAN 
program unit. 

4.1.2 Computed GOTO Statement 

The computed GOTO statement allows the user to transfer control to one of several 
statement numbers, depending upon the value of a variable. The form of the state- 
ment is as follows: 
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GOTO (n ,n , n , . . . , n k ) , ivar 

where n. are statement numbers and ivar is an integer variable whose value is 

between +1 and k. The comma between the statement number list and ivar is 

optional. Control is automatically transferred to the statement whose number is 

n. in the list, 
xvar 

GOTO (75,150,100,85,275)L 

If L is equal to 1 when the statement is executed, control transfers to statement 
number 75 in the program; if L equals 5, control transfers to statement number 
275. If L has a value less than 1 or greater than 5 at the time this statement 
is executed, control passes to the next statement in the program. There is no 
restriction on other uses of the integer variable L within the program. 

A computed GOTO statement may appear anywhere in the executable part of a source 
program except as the terminal statement of a DO loop. The statement numbers in 
parentheses must identify executable statements that exist in the same program 
unit as the GOTO statement. 

4.1.3 Assigned GOTO Statement 

The assigned GOTO statement allows the user to transfer control based on an integer 
variable name. 

Assigned GOTO statements may take one of two forms. 



GOTO K 



GOTO K, (n ]L ,n 2 ,n 3 ,...,1^) 



where K is an integer variable name and n. are statement numbers. Control is 
transferred to the statement whose number is currently associated with the variable 
K via an ASSIGN statement. 

As described in Section 7-3.3, an ASS IGN^ statement defines an integer variable as 
a statement number. Thus, when the statement 

ASSIGN itf TO LOOP 
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is executed, the programmer can transfer control to statement 10 by inserting the 
statement: 

GOTO LOOP 

or, alternatively: , 

GOTO LOOP ,( lit, 20, 100) 

which transfers control to whichever statement number is currently associated with 
LOOP. A statement number must have been assigned to LOOP prior to execution of the 
GOTO statement; if the optional statement list is included, the number assigned to 
LOOP must be one of the listed numbers. Otherwise, control passes to the next 
sequential statement. 

The GOTO statement and the executable statement to which it transfers control must 
both appear in the same FORTRAN program unit. 

4.2 IF STATEMENTS 

IF statements allow program execution to jump from one section of code to another 
or conditionally execute a statement based on the value of an expression. There 
are two types of IF statements. 

a. Arithmetic IF statements 

b. Logical IF statements 

Care should be taken in mixed mode comparisons or comparisons in which one or more 
values have undergone a type conversion. Mixed mode values to be compared are set 
to the more significant of the two types being compared. When this is done, repeating 
binary fractions do not repeat to the full precision of the new mode. For example, 
this may cause a real 0.1 (1.0E-1) to be unequal to a double-precision 0.1 (1.0D-1) . 

4.2.1 Arithmetic IF Statement 
An arithmetic IF statement is of the form 
IF (expression) n ,n ,n 
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where expression is an arithmetic expression in parentheses and n ,n ,n represent 
the statement numbers to which control is transferred if the expression is less than, 
equal to, or greater than zero, respectively. All three statement numbers must be 
present (although two may refer to the same statement) . For example: 

ALPHA = -7 



IF (ALPHA) lj2f, 2j6, lj6 

Control transfers to statement number "Lp, since the expression is less than zero. 
It would also transfer to statement 1)6 if the expression were greater than zero. 
Control would transfer to statement number 20 only if the expression were exactly 
equal to zero. 

Expressions having complex or logical results cannot be used in arithmetic IF 
statements. 

4.2.2 Logical IF Statement 

A logical IF statement is of the form 

IF (expression) statement 

where expression is a logical expression in parentheses and statement is an executable 
FORTRAN statement other than another IF statement or a DO statement. If the 
expression has the value .TRUE, the contained statement is executed; otherwise, the 
next sequential statement is executed. 

Examples of logical IF statements are shown below. 

IF (T.OR.S) X=Y+1 

IF (Z.GT.X(K)) CALL SWITCH (S,Y) 

IF (M.EQ. INDEX) GOTO 15 

4.3 DO STATEMENT 

The DO statement simplifies the coding of repetitive procedures. The DO statement 
is of the form 



DO n i = m ,m ,m 
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where n is a statement number, i is a unsubscripted integer variable, and m m m 

1 '2 3 
are positive integer constants or unsubscripted integer variables. If m is not 

specified, it is understood to be 1. A zero or negative m value is not permitted. 

Statements following the DO statement, up to and including statement n, are 

executed repeatedly for values of i starting with m, and incremented by m until 

1 3 

i is greater than in . 
2 

The statements that are executed as a result of a DO statement are called the 

range of the loop. The variable i is called the inde x variable. The terms m ,m ,m 

12 3 

are the initial , terminal , and increment parameters of the index, respectively. The 

initial parameter value must be less than or equal to the terminal parameter value. 

When the DO statement is executed, its range is first executed with the index equal 
to the initial parameter ' s value . After each execution of the range, the increment 
value is added to the value of the index and the result compared with the terminal 
parameter value. If the index value is not greater than the terminal parameter 
value, and integer' overflow has not occurred, the range is executed again with the 
new index value. The range of a DO loop is always executed at least once. 

For example : 

DO 20 I = 5,100,2 (final iteration has I - 99) 

DO 100 m = 10,100,2 (final iteration has M = 10) 

After the last execution of the range, control passes to the statement immediately 
following the loop. This is tne normal exit from a DO loop. Exit can also be 
accomplished by the execution of a control statement within the range. (Control 
cannot be transferred into the range of a DO loop from outside that loop except 
in the case of an extended range as described below. ) 

The values of the terminal and increment parameters, and of the index variable, of 
the DO loop cannot be altered within the range of the DO statement. The index vari- 
able is, however, available for use as an ordinary variable within the range of a 
DO loop, with a GOTO or IF, the index variable retains its current value and is 
available for use as a variable. The value of an index variable is not defined on 
normal edit from the DO loop. 

The terminal statement of a DO loop cannot be a GOTO, DO, RETURN, or arithmetic 
IF statement. A logical IF statement is allowed as the last statement of the 
range, provided it does not contain any of the statements mentioned above. For 
example : 
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DO 5 K = 1,4 

5 IF (X(K) .GT.Y(L)) Y(K)=X(K) 

6 ... • 

In this case, the loop ends when control passes to statement number 6. Statement 5 
is executed four time whether or not the Y(K)=X(K) statement is executed. Note that 
if statement 5 were 

5 IF (X~(K) .GT.Y(L)) GOTO 10 

improper execution may result. 

The range of a DO statement may include other DO statements, provided the range of 
each DO statement is contained (nested) entirely within the range of the next - 
outermost DO statement. That is, the range of any two DO statements must intersect 
completely or not at all. More than one DO loop within a nest of DO loop can end on 
the same terminal statement. 

The range of a DO loop need not be merely a section of straight line code following 
the DO statement. A control statement that causes execution of instructions else- 
where in the program is permissible as long as control returns to the range of the 
originating loop. The statements executed between the pair of control statements 
(leaving and returning to the DO loop) are called the extended range of the loop. 

DO loops can be nested up to 10 levels deep. In calculating this depth, an implied 
DO in an I/O statement (see Section 7-9.1.3) counts as one level, whose range is 
a single statement; and n implied DO loops within one I/O statement count as n 
levels whose ranges are all within the single statement. For example, 

50 WRITE (5,250) (BUF,F (I) ,1=1, 5) , (DIFF (I) ,1=1,3) 

contains two complete implied loops resulting in one level of nesting at a time. If 
the above statement were contained in an unnested DO loop, the maximum nesting level 
generated by line 50 would be two. 

NOTE 

A common typing mistake is to replace a comma in 
a DO statement with a period. For example: 

DO 5 I = 1.5 

This statement is interpreted by the FORTRAN Compiler 
as. an arithmetic assignment statement: 

D05I =1.5 
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4.4 CONTINUE STATEMENT 

The CONTINUE statement is of the following form. 

CONTINUE 

No processing is done when this statement is executed. It is used primarily as a 
target line number for control transfers, particularly as the terminal statement of 
a DO loop. For example: 

DO 7 K = START, END 



IF (X(K)) 22,13,7 



7 CONTINUE 

A positive value of X(K) begins another execution of the range or causes a normal 
exit from the loop. 

In this example, the CONTINUE statement provides a target address for the IF state- 
ment and ends the range of the DO loop. 

The CONTINUE statement can be inserted anywhere in the executable portion of a 
program with no effect. 

4.5 PAUSE STATEMENT 

The PAUSE statement takes one of the following forms. 

PAUSE 



PAUSE n 



where n is an octal constant from one to six digits long (up to 177777) optionally 
preceded by the letter 0. The PAUSE statement temporarily interrupts program 
execution. 



7-33 



4.6 STOP STATEMENT 



The STOP statement is of the form 



STOP 



STOP n 

where n is an octal constant from one to six digits long (optionally preceded by 
the letter 0) . The STOP statement prints the constant (zeros, if no constant is 
specified) , terminates execution of the current program, and returns control to the 
Monitor. 

4.7 END STATEMENT 

The END statement is of the following form. 

END 

This statement must be the last statement in every FORTRAN program unit. Program 
compilation is terminated when the END statement is encountered in a main program 
unit. 

An END statement forces a call to the EXIT subroutine in a main program. A 
function or subroutine subprogram forces a return to the calling program if a 
RETURN statement is not present. 

If a compiler input file is exhausted before an END statement is encountered, a 
diagnostic is printed and an END statement forced. 
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PART 7 

CHAPTER 5 

SPECIFICATION STATEMENTS 

Specification statements are nonexecutable statements generally placed at the 
beginning of a program to inform the FORTRAN Compiler about storage allocation, 
variable types, subprograms, or data. 

Specification statements must be placed in the program prior to the first usage 
of any variable they describe. Some specification statements determine storage 
allocation and must appear prior to any executable statements. Those statements 
that must precede all executable statements are called declaratives and include: 

IMPLICIT 

EXTERNAL 

type declaration 

COMMON 

EQUIVALENCE 

DIMENSION 

DATA 

The IMPLICIT declaration, if present, should occur first as its action may affect 
any other declaratives. DATA statements, if used must appear last, since other 
declaratives may have an effect on the data elements they define. 

5.1 IMPLICIT STATEMENT 

The IMPLICIT statement is of the form 

IMPLICIT type a l ,a 2 ,(b 1 -b 2 ) / a 3 ,(b 3 -b 4 ),... 

where type indicates the variable type being defined (see Section 7-5.2 for a 
complete list) . The elements a . represent single alphabetic characters separated 
by commas. The elements (b -b ) represent a range of characters denoted by the 
first and last characters of the range and separated by a minus sign. A range of 
characters is always enclosed in parentheses. For example: 

IMPLICIT INTEGER (A-D) , (X-Z) 

This statement instructs the FORTRAN Compiler that any variable name in the program 
beginning with the letters A, B, C, D, X, Y, or Z is an integer variable. 
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The effect of the IMPLICIT statement is to extend the basic rule that all names 
beginning with I, J, K, L, M or N indicate integer variables; all others are real 
variables. Any variable name that does not appear explicitly in a type declaration 
statement and whose first character is one of those listed in an IMPLICIT statement 
is classified according to the type appearing in that IMPLICIT statement. 

The initial state of the FORTRAN Compiler is as though the statements 

IMPLICIT REAL(A-H) , (0-Z) 
IMPLICIT INTEGER (I-N) 

appeared at the beginning of each program. This state is in effect unless a type 
declaration or IMPLICIT statement changes the interpretation of any variable name. 
Definition of a variable name in a type declaration statement overrides any type 
specification determined by an IMPLICIT statement. 

The IMPLICIT statement, if present in a subprogram, must precede all other state- 
ments (except the SUBROUTINE and FUNCTION statement) . The IMPLICIT statement 
becomes effective when it is encountered in the program. Previous specifications 
of variables having implicitly defined values will cause a diagnostic. For 
example , 

COMMON Yl 

IMPLICIT INTEGER (X-Z) 



results in a W (warning) diagnostic. 

Most of the FORTRAN Library function names (FLOAT, IFIX, etc.) follow the default 
type rules. That is, FLOAT produces a real result, IFIX produces an integer 
result, and so on. The use of an IMPLICIT statement may change the default type, 
causing undesirable results. For example, in the program, 

IMPLICIT INTEGER (A-Z) 

REAL X 

1=1 

X=FLOAT(I) 

END 

the function FLOAT is treated as though it produced an integer result leading to 
incorrect values being generated. To have the FLOAT function return a real value 
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in this context, an explicit type declaration should be included in the program 
before or after the IMPLICIT statement, as follows. 



REAL FLOAT 



Those library functions whose names do not follow the default type rules (those 
that return complex or double -precis ion values) are recognized by the compiler 
under any conditions. The data type of these names is not affected by an IMPLICIT 
statement but may be changed by an explicit type declaration. For example: 

INTEGER CABS 

This statement could cause incorrect operations at execution time if the intrinsic 
function CABS is called. 

5.2 TYPE DECLARATION STATEMENT 

A type declaration statement provides the compiler with a description of the 
storage format of the variables defined, and is of the form 

type var lf var 2 ,... 

where type indicates the name of the variable type being defined. Acceptable 
types : 

INTEGER (1- or 2-word integer numbers) 

INTEGER* 2 (same as INTEGER) 

REAL (2-word real numbers) 

REAL*4 (same as REAL) 

DOUBLE PRECISION (4-word real numbers) 

REAL* 8 (same as DOUBLE PRECISION) 

COMPLEX (4-word values, each having a 2-word real 

part and a 2-word imaginary part) 

LOGICAL (true = -1, false = 0) 

BYTE (8-bit format) 

LOGICAL* 1 (same as BYTE) 

Allocation of two bytes for 1-word integers is effected only by the /ON switch. 
INTEGER* 2 causes four bytes per element when the /ON switch is not set (see 
Section 7-12.1) . 
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Allocation of two bytes for LOGICAL variables is effected only by the /ON switch. 
LOGICAL allocates four bytes per element when the /ON switch is not set (see 
Section 7-lj0f. 2) . 

The terms var ,var , ... indicate the variable narae(s) that are to be allowed 
storage space in the data type format specified. 

Type declarations must precede all executable statements. For example: 

COMPLEX RET,X,TST 

BYTE IAMB , LEN , R ( 20) , S ( 20) 

In the first statement, three variable names are specified as complex variables. 
In the second statement, two variable names are specified as 8-bit byte variables, 
and the arrays R and S are dimensioned as having byte-sized elements each. 

Type declarations can be used to define dimension specifications for arrays. 
However, an array defined in a type declaration cannot be redimensioned by any 
other specification statement, and vice versa. Adjustable arrays in subprograms 
can also be defined in type declarations. 

If a variable name is defined in a type declaration, the type so specified over- 
rides any type specification determined by the initial letters of the variable. 
The priority for defining a variable type is, therefore, as follows: 

a. type declarations. 

b. IMPLICIT variable type. 

c. real or integer, depending on the initial letter (where I through N 
indicates an integer) . 

5.3 DIMENSION STATEMENT 

The DIMENSION statement is of the form 



DIMENSION S ,S ,S ,... 
1 ^ 3 



where S is an array specification of the form 



var (n ) 
var(n lf n 2 ) 



or 



var(n ,n 2 ,n 3 ) 
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• 



where var is the array variable name, and n ,n , and n are unsigned integer 
constants (or variables, optionally, in the case of adjustable array dimensions 
in subprograms) indicating the maximum size of the array in each of the three 
possible dimensions. 

An array can have one, two, or three dimensions. For example: 

DIMENSION A (20), B(10,10), C(5,5,5) 

Where more than one array is defined in a single DIMENSION statement, the array 
specifications are separated by commas. 

DIMENSION statements must precede all executable statements. 

The appearance of the DIMENSION statement causes the FORTRAN Compiler to allocate 
space for the appropriate number of variables. Array definition can also be 
performed by a COMMON statement or any type declaration statement. For example: 

•COMMON X(10,4) ,Y,Z 
INTEGER A (7, 32) 

Once defined, an array cannot be redimensioned by any other specification statement. 

Unless an array is specified in a type declaration statement, the type of the array 
is determined by the first letter of the array variable name. 

An error message is always generated at run-time if an array subscript is less than 
one. Optionally, using the /CK switch option (see Appendix J) , the system will 
check to see that array subscripts do not exceed the maximum size specified. 

A subprogram can use adjustable arrays, the size of which may depend upon the 
parameters of the subprogram call (see Sections 7-6.4.1 and 7-^6.4.2). 

5.3.1 Array Storage 

Arrays of one, two, or three dimensions are each stored as a linear sequence of values 
in memory. Arrays are always stored in me: 
element of a 1-dimensional array is A(n) . 



in memory. Arrays are always stored in memory column by column. Therefore, the n th 



In order to find the linear element of a 2-dimensional array A (I, J) where a given 
element is A(i,j), the formula 
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n * i + I*(j-1) 

is evaluated. A(i,j) is stored at the n element. 

In the 3-dimensional array (A(I,J,K) where a given element is A(i,j,k), an indi- 
vidual element is stored at the n element, where n is 

n - i + I*(j-1) + I*J (k-1) . 

It is sometimes desirable to equivalence a 2- or 3-dimensional array to a 1-dimensional 
array. The above formula is useful for determining the algorithms for manipulating 
the arrays after such a change. However, for any given array, elements of that 
array must be referenced using the number of subscripts specified in the defining 
statement (although in some contexts an entire array can be referenced by indicating 
the array name only) . See Figure 7-2. 

5.4 COMMON STATEMENT 

The COMMON statement causes the specified variables to be stored in an area of memory 
available to other programs and subprograms. Common block areas can be of any size 
and any number of them can exist within the space available on the system. One 
unnamed common block (known as blank common) is allowed. The variables are placed 
in the individual common blocks in the order in which they appear in the block 
description. 

The COMMON statement takes the form 

COMMON /blockl/a,b,c/block2/d,e,f 

where blockn is a 1- to 6-character alphanumeric name (the first character of 
which must be alphabetic), delimited by slashes, and considered to be the block 
name. The sequences a,b,c and d,e,f represent the names of the variables assigned 
to the common block (s) . The form 

COMMON a,b,c,. . . ,n 

is used to indicate the blank common block. 

One or more common blocks can be defined in the same COMMON statement. 

For example: 

COMMON /A/X,Y,Z/B/L,M,N 
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This statement defines two common blocks, each containing three variables. 

Where blank common is the first common block being described in the COMMON statement, 
no block name or specification is necessary. Subsequent block names in the same 
COMMON statement are delimited by slashes. However, when elements of blank common 
are indicated other than in the first position in a COMMON statement, two con- 
secutive slashes (any separating spaces are ignored) are used. For example: 

COMMON /R/X,Y//B,C,D 
COMMON B,C,D/R/X,Y 

The two COMMON statements above are identical in effect. 

All COMMON statements within a program unit must precede all executable statements. 

A common variable cannot appear in more than one common block in a single program 
unit. However, since labeled common block names are used only by the compiler and 
are not present in the functional program at run time, such block names can be used 
as variable or subroutine names within the program. 

Block entries are linked sequentially throughout the program, beginning with the 
first COMMON statement. For example, the statements 

COMMON/A/ALPHA/B/BET , GAM/C/ET 
COMMON/B/DEL , EP , ZE/C/THE , IO/D/KA , KA 

have the same effect as the statement 

COMMON/A/ALPHA/B/BET , GAM , DEL , EP , ZE/C/ET , THE , IO/D/KA , LA 

Storage allocation for blocks of the same name begins at the same location for all 
program units linked together in a single load module. For example, if a program 
contains the statement 

COMMON A,B/R/X,Y/Z 

as its first COMMON statement, and a subprogram contains the statement 

COMMON /R/U,V,W//D,E 
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as its first COMMON statement, the quantities represented by X and U, by Y and V, 
and by Z and W (provided they are of corresponding word lengths) are stored at the 
same locations in common block R; a similar correspondence holds for the quantities 
A and D, and B and E, in the blank common block. Responsibility for the ,values of 
the variables used and the correspondence of word lengths of the variables rests 
with the programmer. 

The size of a common block is the sum of the storage required for all elements 
introduced through COMMON and EQUIVALENCE statements. Where the size of the 
various common blocks referenced by several program units in a single load module 
differs, the length of each common block is resolved at linkage time to be the 
length of the longest block of the given name declared by any program unit. 



The variable names that follow the block name indicate scalar or array variables 
assigned to that block. Where an unsubscripted array name appears in a COMMON 
statement, the array must have been previously defined in a DIMENSION or type 
declaration statement. Where a subscripted array name appears in a COMMON state- 
ment and the array has not been defined, the array is defined as being the indi- 
cated size. Where the array has been dimensioned previously, an error diagnostic 
is printed. For example: 

DIMENSION T (10) 
COMMON ABS( 5, 10,5), T 

The COMMON statement defines the array ABS and allocates storage space for that 
array in blank common. All elements of the array T are also stored in blank 
common. (Arrays in common storage are stored according to the rules specified 
in Section 7-5.3.1.) 

When both BYTE (LOGICAL*l) and other types of variables are included in the same 
COMMON block, care must be taken to assure that non-BYTE elements that follow 
BYTE elements are allocated on a word boundary. Non-BYTE elements may only follow 
an even number of BYTE elements. (Failure to assure this may result in a "word 
reference to odd address" hardware error when the program is executed.) 

5.5 EQUIVALENCE STATEMENT 

The EQUIVALENCE statement is of the form 

EQUIVALENCE (var ,var ,...), (var. ,var .,...),.. . 
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where var indicates a variable name used in the current program unit. Each of the 
variables within a set of parentheses is assigned to the same storage location. 
For example, 

EQUIVALENCE (LAN, MAR) 

specifies that the variables LAN and MAR are stored in the same location, and 
have the same value. Where variables being equivalenced are not of the same 
type, the equivalencing operation assumes that each variable begins in the same 
word (see the description of storage formats in Chapter 12) . 

EQUIVALENCE statements must precede all executable statements. 

The equivalence relationship is transitive; that is, the two statements below have 
the same effect. 

EQUIVALENCE (A,B) , (B,C) 
EQUIVALENCE (A,B,C) 

The EQUIVALENCE statement is used to conserve core storage during execution or to 
obtain a portion of a numeric value in a different format. Rather than reserving 
space for a variable that is only used for a short time in the program, a later vari- 
able can be equivalenced to use this space when the earlier value is no longer of 
interest. A statement such as the following 

EQUIVALENCE (TRI,ISO) 

allows ISO to represent the storage location formerly occupied by the real variable 

TRI. 

To make full use of the space required by a real variable, that real variable can 
be equivalenced to two integer variables (assuming use of the /ON switch 1 which 
causes 1-word integers to be used throughout the program) . Similar techniques 
can be used to equivalence two real, or two to four integer, variables to complex 
or double-precision variables. For example: 

DIMENSION N(2) 
EQUIVALENCE (TRI,N) 



See Appendix J. 
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These two statements equivalence the first word of TRI to N(l) and the second word 
of TRI to N(2), assuming use of the /ON switch. 

5.5.1 Equivalencing Array Variables 

The subscripts of any array variables appearing in an EQUIVALENCE statement must be 
integer constants. For example: 

EQUIVALENCE (X,A(5) ), (BET (2,2) ,Y (2,4,1) ,NET) 

A previously-dimensioned array name used in an EQUIVALENCE statement without sub- 
scripts implies reference to the first element of the array. For example, 

DIMENSION Y(20) 
EQUIVALENCE (Y,X) 

has the same effect as 

DIMENSION Y(20) 
EQUIVALENCE (Y(1),X) . 

Where a subscript is present in an EQUIVALENCE statement, the array element speci- 
fication must contain an integer constant for each of the array dimensions as 
originally defined. No default value is supplied for missing dimension elements, 
and a diagnostic is generated. For example, 

DIMENSION REAL (2, 10, 10) 
EQUIVALENCE ( REAL (1,1), TRAN ) 

is illegal and would result in an error message. (See also the description of array 
storage in Section 7-5.3.1). 

5.5.2 EQUIVALENCE and COMMON Interaction 

The same variable name may appear in both COMMON and EQUIVALENCE statements. How- 
ever, two quantities in common cannot be made equivalent to each other. Quantities 
placed in a common block by an EQUIVALENCE statement may cause the end of the 
common block to be extended. For example, 

COMMON/R/X,Y,Z 
DIMENSION A (4) 
EQUIVALENCE (A,Y) 
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causes the common block R to extend from X to A (4) as follows. 
x 

Y-A(l) 
Z=A(2) 

A(3) 

A(4) 

The end of the common block has been extended by two elements. No EQUIVALENCE 
statements are allowed to extend the beginning of a common block, however. For 
example, the sequence 

COMMON/R/X,Y,Z 
DIMENSION A (4) 
EQUIVALENCE (X,A(4)) 

is not permitted, since it would require A(l) , A (2) and A (3) to extend beyond the 
starting location of common block R. 

An EQUIVALENCE statement is not permitted to violate any previous EQUIVALENCE 
storage assignments or extend any array beyond its maximum predefined dimensions. 
For example: 

DIMENSION A(10),B(10) 

EQUIVALENCE (A(10) ,B(1) ) , (A(5) ,B(3) ) 

EQUIVALENCE (A(5,l$) ,B (5) ) 

Both of the preceding EQUIVALENCE statements are unacceptable. 

5.5.3 EQUIVALENCE and BYTE Arrays 

In using EQUIVALENCE with items in a BYTE (LOGICAL* 1) array, only those items 
aligned on a word boundary can be specified in an EQUIVALENCE statement. Individual 
BYTE variaoles are always aligned on word boundaries except within COMMON blocks 
(see Section 7-5.4). In a BYTE array, only the elements having odd subscripts are 
aligned on word boundaries. For example: 



BYTE A,B(9) 
REAL X 



correct : EQUIVALENCE ( A , X , B ( 3 ) ) 

incorrect: EQUIVALENCE (A,X ,B (4) ) 
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5.6 EXTERNAL STATEMENT 

The EXTERNAL statement is of the form 

EXTERNAL a,b,C,.. . 

where a,b,c,... represent the names of FUNCTION or SUBROUTINE subprograms to be used 
as parameters to a subprogram call. 

FUNCTION and SUBROUTINE subprogram names can be used as arguments of other sub- 
programs. When used in this manner, subprogram names are distinguished from 
ordinary variables by their appearance in an EXTERNAL statement. For example: 

EXTERNAL SIN, COS 
CALL TRIGF (SIN, 1,5, ANSWER) 
. CALL TRIGF (COS, 18. 7, ANSWER) 



END 



SUBROUTINE TRIGF (FUNC,ARG, ANSWER) 



ANSWER - FUNC(ARG) 

RETURN 

END 

The EXTERNAL statement must appear before any executable statements in any program 
or subprogram. This statement is used in programs or subprograms using the name 
of a function or subroutine as a parameter in a subprogram call. This statement 
does not take a line number and should not contain any name previously defined. 

5.7 DATA STATEMENT 

The DATA statement is of the form 

DATA var/value/, var /value/, .. . 
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where var is a list of variable names, and value is an ordered list of numeric or 
string values that are to be assigned to the respective elements of the variable 
list. 

The DATA statement is used to supply initial or constant values for variables. 
DATA statements must precede all executable statements and must follow all other 
declaratives. Elements of the variable and value lists are separated by commas. 
More than one set of variable and value lists can be contained in a single DATA 
statement. For example: 

DATA A,B,C/3.4,5.9,7./ 

DATA I,J/7,3/X,Y,Z/43.5,2.01,8.4/ 

Dummy arguments must not appear in the variable list. Whenever several variables 
are to be assigned the same value, the item in the value list can be preceded by 
an integer constant indicating the number of repetitions of that variable. The 
integer and data values are separated by an asterisk. For example: 

DATA A,B,C/3*3.17/ 

The three real variables are each assigned the value 3.17. 

Radix- 50 constants (See Section 7-2.2.9) may only be used to initialize variables 
of type REAL. 

Hollerith (or character) constants may be used to initialize variables of any type. 
Each Hollerith constant may be used to initialize exactly one variable. If the number 
of characters in a Hollerith constant is fewer than needed to initialize the entire 
variable element, the compiler will append additional blank characters to the right 
of the constant to completely fill the variable. If the number of characters in 
the constant is more than needed, the leftmost characters will be used and the 
remainder ignored. 

The specified values are compiled into the object program and are assumed by the 
variables when program execution begins. Variables in a labeled common block can 
only be defined in a BLOCK DATA subprogram (see Section 7-6.4.8). Variables in 
the blank common block cannot be initialized. 
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Variables may also be subscripted array elements or unsubscripted array names. 
When the unsubscripted name of a predefined array is given in the list, a data 
value must be specified for every element of the array. Data elements are stored 
in the array in the order used for data transmission and storage of arrays; that is, 
in order of increasing subscripts with the first subscript varying most rapidly 
(column by column storage, see Section 7-5.3.1). Allocation of data to an array 
stops when 

a. the data (value) list is exhausted; or 

b. data items have been allocated to the entire array, in which case 
additional data items are allocated to any additional items in the 
variable list. 

For example, 

DATA XfYjZ/'A', 'BCDE 1 , 'FGHIJKL'/ 



produces the following in memory. 

X+l 



blank 


A 



X+3 



X+2 



blank 


blank 




Y+l 


Y 




C 


B 


Y+3 


Y+2 




E 


D 


Z+l 


Z 




G 


F 


I 


Z+3 


Z+2. 




I 


H 


1 



The data items following each list of variables must have a one-to-one corres- 
pondence with the variables of the list, and must agree in type, since each item 
of the data specifies the value given to the corresponding variable. 
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Data items assigned may be numeric, Hollerith, octal, hexadecimal, Radix-50, or 
logical constants. For example, 

DATA ALPHA , BETA , IVAL/5 . , 16 . E- 2 , . TRUE . / 

specifies the value 5.0 for ALPHA, the value 0.16 for BETA, and the value -1 for 
IVAL. Any item of data may be preceded by an integer constant followed by an 
asterisk. This notation indicates that the item is to be repeated. Example: 

DATA A(l) ,A(2) ,A(3)/3*0./ 
specifies the value 0.0 for array elements A(l) to A(3) . For example, 

DIMENSION A(2, 2), B(3) 

DATA A,B/2*1 .0,3*2.0,3.0,4.0/ 

will initialize 

A(l,l) and A(2,l) to 1.0 
A(l,2), A(2,2) and B(l) to 2.0 
B(2) to 3.0, and B(3) to 4.0 

As another example, 

DIMENSION R(4) ,1(3) ,L(5) 

DATA L,I,R/4*0177777,2* , AB', 2*4,4*1.0/ 

results in 

L(l) to L(4) = 177777 (octal) 
L(5),I(1) = 'AB' (Hollerith) 
1(2), 1(3) = 4 
R(l) to R(4)= 1.0 
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PART 7 

CHAPTER 6 

FORTRAN FUNCTIONS AND SUBROUTINES 

6.1 PROCEDURES 

FORTRAN contains provisions for the creation and use of coded procedures. Pro- 
cedures consist of FORTRAN (or assembly language) routines that perform operations 
for, and independently of , the main FORTRAN program. Procedures perform operations 
that may be required more than once in the course of a program and that might 
otherwise require repetition of code within a program. Their use, then generally 
saves programming effort and reduces storage used/ 

There are two main types of procedures: functions and subroutines. Functions 
include the following. 

1. library functions (provided by the FORTRAN OTS) 

2. statement functions (one-line, user-defined function definitions) 

3. external functions (user-defined function subprograms) 

Subroutines include the following. 

1. system subroutines (provided by the FORTRAN OTS) 

2. subroutines (user-defined subroutine subprograms) 

Library functions and system subroutines are provided by the FORTRAN Object Time 
System and need only be referred to or called by the user program. (See Section 
7-6.2 and Chapter 7-14 for a description of the available routines.) External 
functions and subroutines can be compiled independently of the main program; 
they are linked to the main program by the Linker Program (LINK) . 

The term subprogram refers to both external functions and subroutines. Sub- 
programs can communicate with the main program and among themselves by means of 
parameters specified in the subprogram call or by means of Variables stored in 
common blocks. Calls to subprograms cannot be recursive; that is, a call to any sub- 
program cannot result in a subsequent call to that subprogram, directly or indirectly. 
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6.2 FORTRAN LIBRARY FUNCTIONS 

The FORTRAN Library functions are described in Table 7-3. In order to use a library 
function in any routine, it is only necessary to use the symbolic name of the 
function, together with one or more data references (arguments) upon which the 
function is to act, in a FORTRAN statement. The value obtained from the execution 
of the function is assigned to the function's name, which can then be used as a 
variable in the execution of the FORTRAN statement. For example, 

R = ABS(X-l) 

causes the absolute value of X-l to be calculated and assigned to the variable R. 

The data type of each library function is predefined as described in Table 7-3. 
Arguments passed to these functions may consist of subscripted or simple variable 
names, constants, arithmetic expressions or other intrinsic functions. Arguments 
to these functions must correspond to the type indicated in Table 7-3. If the 
argument types do not match, no conversion is performed, nor is any diagnostic 
given. See also Section 7-5.1 for the effect of the IMPLICIT statement on function 
values. 

Details on the algorithms used in the individual library functions are contained in 
Chapter 7-13. 

6.3 ARITHMETIC STATEMENT FUNCTIONS 

An arithmetic statement function is used to define one-line functions for the user. 
These functions can be referenced only by the program unit in which they are 
defined. 

To define an arithmetic statement function, the following form is used: 

name (arg , . . . ) = expression 

where name is the function name, which must be a legal FORTRAN variable name. 
Function type is determined either implicitly, by the initial letter of the function 
name, or explicitly, through a data type specification. The term (arg , ... ) 
represents the list of dummy variables used in the function definition; at least 
one argument must be present. Expression is an arithmetic expression that defines 
the computation to be performed by the function. 
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Table 7-3 



FORTRAN Library Arithmetic Functions 









Argument 


Result 


Form 


Definition 




Type 


Type 


ABS(X) 


Real absolute value 




Rea,l 


Real 


IABS(I) 


Integer absolute value 




Integer 


Integer 


DABS (X) 


Double precision absolute value 




Double 


Double 


CABS (Z) 


Complex to Real, absolute value 










where Z=(X,Y) CABS (Z)= ^2^2^/2 




Complex 


Real 


FLOAT (I) 


Integer to Real conversion 




Integer 


Real 


IFIX(X) 


Real to Integer conversion 
IFIX(X) is equivalent to INT(X) 




Real 


Integer 


SNGL(X) 


Double to Real conversion 




Double 


Real 


DBLE(X) 


Real to Double conversion 




Real 


Double 


REAL(Z) 


Complex to Real conversion, 
obtain real part 




Complex 


Real 


AIMAG(Z) 


Complex to Real conversion, 
obtain imaginary part 




Complex 


Real 


CMPLX(X,Y) 


Real to Complex conversion 
CMPLX(X,Y)=X+Y*i 




Real 


Complex 




Truncation functions return the largest 








integer £ | argument | , carrying the 


sign 








of the argument. 








AINT (X) 


Real to Real truncation 




Real 


Real 


INT(X) 


Real to Integer truncation 




Real 


Integer 


IDINT (X) 


Double to Integer truncation 




Double 


Integer 




Remainder functions divide the first 








argument by the second and return 1 


the 








remainder from that division 








AMOD(X,Y) 


Real remainder 




Real 


Real 


MOD (I, J) 


Integer remainder 




Integer 


Integer 


DMOD(X,Y) 


Double precision remainder 




Double 


Double 




Maximum value functions return the 


largest 








value in a list of at least two arguments. 






•AMAXJZf(I,J, ...) 


Real maximum from Integer list 




Integer 


Real 


AMAX1(X,Y,...) 


Real maximum from Real list 




Real 


Real 


MAX#(I,J,...) 


Integer maximum from Integer list 




Integer 


Integer . 


MAX1(X,Y,...) 


Integer maximum from Real List 




Real 


Integer 


DMAXKX^Y,...), 


Double maximum from Double list 




Double 


Double 
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Table 7-3 (Cont.) 



FORTRAN Library Arithmetic Functions 



Form 



AMINjef(I,J,...) 

AMINl(X,Y,...) 

MIN$(I,J,...) 

MIN1(X,Y,...) 

DMIN1(X,Y,...) 



SIGN(X,Y) 

ISIGN(I,J) 

DSIGN(X,Y) 



Definition 



Argument 
Type 



Result 
Type 



Minimum value functions return the 
smallest value in a list of at least 
two arguments. 

Real minimum from Integer list 
Real minimum from Real list 
Integer minimum from Integer list 
Integer minimum from Real List 
Double minimum from Double list 



Integer 


Real 


Real 


Real 


Integer 


Integer 


Real 


Integer 


Double 


Double 



The transfer of sign Functions return 
the value of the first argument to which 
the sign of the second has been attached. 

Real transfer of sign 

Integer transfer of sign 

Double precision transfer of sign 



Real 


Real 


Integer 


Integer 


Double 


Double 



DIM(X,Y) 
IDIM(I,J) 



Positive difference functions return 
the first argument minus the lesser 
of the two arguments (never less than 
zero) . 

Real positive difference 
Integer positive difference 



Real 
Integer 



Real 
Integer 



EXP(X) 
DEXP (X) 
CEXP (Z) 



Exponential functions return the value 
of e raised to the power of the argument . 



Real 

Double 

Complex 



Real 

Double 

Complex 



ALOG(X) 
ALQGIJZJ (X) 
DLOG(X) 
DLOG10(X) 
CLOG 



SQRT (X) 

DSQRT(X) 

CSQRT(Z) 



SIN(X) 

DSIN(X) 

CSIN(Z) 



Returns log (x) 

Returns log. (X) 

Returns log (X) 

Returns log-_(X) 

Returns log of complex argument 



Real 

Real 

Double 

Double 

Complex 



Real 

Real 

Double 

Double 

Complex 



Square root of Real argument Real Real 

Square root of Double precision argument Double Double 
Square root of Complex argument Complex Complex 



Real sine 

Double precision sine 

Complex sine 



Real 

Double 

Complex 



Real 

Double 

Complex 
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Table 7-3 (Cont.) 
FORTRAN Library Arithmetic Functions 







Argument 


Result 


Form 


Definition 


Type 


Type 


COS (X) 


Real cosine 


Real 


Real 


DCOS (X) 


Double precision cosine 


Double 


Double 


CCOS (Z) 


Complex cosine 


Complex 


Complex 


TANH(X) 


Hyperbolic tangent 


Real 


Real 


ATAN (X) 


Real arctangent 


Real 


Real 


DATAN(X) 


Double precision arctangent 


Double 


Double 


ATAN2(X,Y) 


Real arctangent of (X/Y) 


Real 


Real 


DATAN2(X,Y) 


Double precision arctangent of (X/Y) 


Double 


Double 


CONJG(Z) 


Complex conjugate, if Z=X+Y*i 
CONJG(Z)=X-Y*i , 


Complex 


Complex 


RAN(I,J) 


Returns a random number of uniform 
distribution over the range j6 to 1. 
I and J must be integer variables and 
should be set initially to j6. Resetting 
I and J to regenerates the random 
number sequence. Alternate starting 
values for I and J will generate 
different random number sequences. 


Integer 


Real 



Arguments of the sine and cosine functions represent angles expressed in radians. 
Arithmetic statement function definitions must occur after* all declaratives (speci- 
fication statements) and before any executable statements. 

An arithmetic statement function definition must be contained in a single statement. 
The expression that defines the function may include dummy arguments, other variables, 
array elements, non-Hollerith constants,' and references to intrinsic functions, 
external functions, and previously-defined statement functions. For example, 

ACOSH(X,A) = (EXP(X/A))/2.0 

is an acceptable arithmetic statement function definition. X and A are dummy argu- 
ments. Since they serve only to indicate the number, type, and order of actual 
arguments, they may also be used as variables els'ewhere in the program unit. 
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Array references cannot be used as dummy arguments in a statement function definition. 
For example, 

DIMENSION X(10) 

F(X) = Y+X(D+X(2)+X(3) 

would generate a function reference to the array X instead of an array reference. 
No diagnostic is issued. However, in the following definition, 

G(Y) = F(Y)*X(Y+2) 

X is legal as an array name since it is not a dummy argument. 

6.4 EXTERNAL SUBPROGRAMS 

6.4.1 Subprogram Arguments 

Variable names and identifiers used in any given program unit are completely inde- 
pendent of variable names in other program units. Relationship between data 
references in different program units (and, thereby, communication of data between 
those program units) is established by means of arguments. The variable names and 
identifiers that appear in a function reference or a subroutine CALL statement are 
called actual arguments; those that appear in a FUNCTION or SUBROUTINE statement 
that defines an external subprogram are called dummy arguments. The rules con- 
cerning dummy arguments are identical for both function and subroutine sub- 
programs . 

Subprograms are written much like a FORTRAN main program. However, the first 
statement of a subprogram must be either a FUNCTION statement with a list of 
arguments or a SUBROUTINE statement with or without an argument list. The argu- 
ment list contains the dummy arguments of the subprogram. When the subprogram is 
called, the dummy arguments are replaced by the actual arguments contained in the 
calling statement. 

A subroutine may return zero, one, or more values to the calling program by means 
of its arguments. A function returns a single value, which is assigned to the 
function's name. The values of a function's arguments may, however, be altered 
within the function to effectively return additional values to the calling program 
unit. Actual arguments that are passed to a subprogram may be any of the following. 

arithmetic or logical expressions 
alphanumeric strings 
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array names or elements 
subprogram names 

(Expressions may contain constants, variables, array elements, and library function 
references, alone or in conjunction with operators.) 

Actual arguments must agree in number, order, and type .with the dummy arguments of 
the subprogram to which they are passed. No type conversion is performed where a 
difference in data type occurs, nor is any diagnostic given. 

Dummy arguments that represent array names must be defined within the subprogram, 
either i n a DIMENSION statement or in one of the type specification statements 
that can contain dimensioning information. Array dimensions within a subprogram, 
although they need not be exactly the same as those of the corresponding array (s) in 
the calling program unit, must not exceed those dimensions. Dummy arguments may 
be used as dimension specifications within a subprogram; the actual dimensions may 
then be passed to the subprogram by the calling program unit (see Section 7-6.4.2). 

Dummy arguments must not appear in EQUIVALENCE, DATA, or COMMON statements within a 
subprogram. 

A function or subroutine subprogram may be used as. an actual argument of a sub- 
program, provided that name appears in an EXTERNAL statement in the calling program 
unit. The EXTERNAL statement identifies such a name as being an external subprogram 
name, rather than a variable name. 

For example, the subprograms TRIG and TAN: 

SUBROUTINE TRIG (ARG,FUNC, RES) 

RES = FUNC(ARG) 

RETURN 

END 

FUNCTION TAN(ARG) 

TAN = SIN (ARG) /COS (ARG) 

RETURN 

END 

If the calling program unit contains an EXTERNAL statement naming the library 
functions SIN and COS, and the user-written function TAN, three different calls 
to the subroutine TRIG can be made to produce three different results. For 
example: 
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EXTERNAL SIN, COS, TAN 



CALL TRIG (ANGLE, SIN, SINE) 



CALL TRIG (ANGLE, COS, COSINE) 



CALL TRIG (ANGLE, TAN, TANGNT) 

Because each call transmits a different function name to the subroutine TRIG, the 
first call returns the sine of the first argument, the second call returns the 
cosine, and the third, the tangent. 

6.4.2 Adjustable Dimensions 

In order to obtain maximum flexibility from a subprogram, the user may wish to vary 
the dimensions of one or more arrays with each subprogram call. An array within 
a subprogram may have adjustable dimensions, provided it refers to an array defined 
in the calling program unit and passed as an argument. 

In this case, both the array name and its dimension values must be expressed as 
dummy arguments in the FUNCTION or SUBROUTINE statement that begins the subprogram. 
The array name (which may be of any data type) and the dimension variables (which 
mus€ be simple integer variables) are used to define the array within the sub- 
program as follows. 

SUBROUTINE TRAN (A, I , J,RET,DET) 
DIMENSION A-(I, J) 



RETURN 
END 

The dimension arguments passed in this way must not exceed the actual dimensions 
of the corresponding array, as defined in the calling program unit, although they 
may be smaller. 

The arguments representing the array dimensions must be defined prior to the subpro- 
gram call; the dimension values will not change "during the execution of the 
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subprogram even though the dummy arguments may be redefined or made undefined within 
the subprogram. The dimension specifications for an adjustable array must be supplied 
as arguments within the subprogram call; they cannot reside in common. 

6.4.3 FUNCTION Subprograms 

A FUNCTION subprogram is a computational procedure that returns a single value through 
its name; a function can also return values through its arguments. A FUNCTION subpro- 
gram is invoked by using the function name with associated arguments i n an 
arithmetic expression. A FUNCTION subprogram begins with a FUNCTION statement and 
ends with an END statement. Control is returned to the calling statement by means 
of one or more RETURN statements. (The END statement acts like an implied RETURN 
statement.) The mode of the function (the type of value returned) is determined 
by the name of the function or a type designation inserted into the initial FUNCTION 
statement . 
The FUNCTION statement format: 

FUNCTION name ( arg. ,. ..arg ) 
In 



or 



type FUNCTION name (arg , ...,arg ) 

This statement declares the program unit that follows to be a FUNCTION subprogram. 
The function name may be any legal FORTRAN variable name that is not used as a 
dummy argument. This name must not appear in any nonexecutable statement in the 
calling program or function subprogram, except as a scalar variable in a type declara- 
tion. At some point in the subprogram the function name must appear as a scalar 
variable to which a value is given by an assignment statement. The value assigned 
to the function name is the value returned by the function reference in the calling 
program. 

The argument list in the FUNCTION statement must consist of at least one dummy 
argument. Dummy arguments follow the rules established for subprogram parameters 
as described in Section 7-6.4.1 and 7-6.4.2. 

The data type of the value returned by the function is normally determined to be 
either INTEGER or REAL according to the default type rule applied to the function 
name. However, a different type can be defined for the returned value by placing 
a type designation term at the beginning of the FUNCTION statement or placing a 
type declaration statement within the subprogram. 
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A FUNCTION statement must not have a statement number. The only FORTRAN statements 
that must not appear in a FUNCTION subprogram are SUBROUTINE, BLOCK DATA and 
another FUNCTION statement. 



An example of a FUNCTION subprogram is shown in Figure 7-3. 



0001 



0002 
0003 



0004 



0001 



0002 

0003 
0004 
0005 
0006 
0007 
0008 



FRAGMENTARY FORTRAN PROGRAM 
FOR ADJUSTABLE ARRAYS 

REAL X(10,10), Y(25,25) 



XT=TRACE(X,10) 
YT=TRACE(Y,25) 



END 



ROUTINES CALLED: 
TRACE 

BLOCK LENGTH 

MAIN 1500 (005670) * 

**COMPILER CORE** 

PHASE USED FREE 
DECLARATIVES 00366 17629 
EXECUTABLES 00446 17549 
ASSEMBLY 00881 19977 



10 



BLOCK 
TRACE 



REAL FUNCTION TRACE ( A, N) 

COMPUTE THE TRACE OF N * N MATRIX A 

DIMENSION A(N,N) 

SUM=0.0 

DO 10 1=1 ,N 

SUM=SUM + A (I, I) 

TRACE=SUM 

RETURN 

END 

IjENGTH 
82 (000244)* 



**COMPILER CORE** 

PHASE USED FREE 
DECLARATIVES 00366 17629 
EXECUTABLES 00446 17549 
ASSEMBLY 00929 19929 



Figure 7-3 
Sample FUNCTION Subprogram 
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6.4.4 SUBROUTINE Subprograms 



A SUBROUTINE subprogram is a computational procedure that returns zero, one, or more 
values to the calling program. Values and data types are associated with the sub- 
routine arguments, but not with the subroutine name. A SUBROUTINE subprogram 
begins with a SUBROUTINE statement and ends with an END statement. Control is 
returned to the calling statement by means of one or more RETURN statements. 



The SUBROUTINE statement is of the following form. 



SUBROUTINE name 



or 



SUBROUTINE name(arg, , . . . ,arg ) 
1 n 

This statement declares the program that follows to be a SUBROUTINE subprogram. The 
subroutine name may be any legal FORTRAN variable name, but must not be used as a 
variable name within the calling program or subprogram itself, or as the name of a 
labeled COMMON block used in the same load module. 

The argument list in the SUBROUTINE statement consists of any number of dummy 
arguments (or none) . Dummy arguments follow the rules established for subprogram 
parameters as described in Section 7-6.4.1 and 7-6.4.2. Dummy arguments to a 
SUBROUTINE subprogram must not appear in an EQUIVALENCE or COMMON statement within 
the subprogram. 

A SUBROUTINE subprogram can redefine any of the actual arguments or use any of the 
dummy arguments specified to return results to the calling program. For example: 

SUBROUTINE COMPUT (A,B,ANS) 

ANS=A+B 

A=A+1.0 

RETURN 

END 

In this simple subroutine, the calling statement provides numeric values for A and 
B and a variable name in which the value of ANS is returned. If the calling state- 
ment is 

CALL COMPUT (AI, 4.0, RET) 

where the value of AI is 7.0, the value returned to the calling program in RET is 
11.0, and AI is 8.0. 
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A SUBROUTINE subprogram cannot contain a FUNCTION, BLOCK DATA or another SUBROUTINE 
statement. 

6.4.5 CALL Statement 

A subroutine is invoked by means of a CALL statement having one of the following 
forms : 

CALL name 



CALL name(arg„ , . . .arg ) 
1 n 

where name is the name of the subroutine referenced.. The argument list of the 
subroutine call, where included, should agree in number, order, and type with the 
argdments in the subroutine definition. The actual arguments specified in the 
subroutine call replace the dummy arguments during execution of the subroutine. 

The arguments of a CALL statement may be expressions, array names or elements, 
scalar variables, alphanumeric strings, or subprogram names. The subroutine name 
cannot be referenced as a basic element in an expression in the calling program 
or in the subroutine itself. 

A call to a subroutine may contain fewer, or more, actual arguments that the sub- 
routine definition, as long as 

a. an unsupplied call argument is not a dummy array in the subroutine, and 

b. some indication of the number of real arguments (such as an argument 
that states the number of other arguments supplied on a given subroutine 
call) is given to the subroutine so that no attempt is made to reference 
unsupplied arguments. 

6.4.6 RETURN Statement 

The RETURN statement is of the form 

RETURN 

and causes an exit from the subprogram to the calling program. It is generally the 
last executable statement in the subprogram. Any number of RETURN statements can 
appear in either a FUNCTION or SUBROUTINE program. 
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The RETURN statement must not occur in a main program. 

En END statement must occur following each subprogram. If a subprogram contains an 
END statement and no RETURN statement, executing the END statement causes a return to 
the calling program. 

6.4.7 BLOCK DATA Subprograms 

The BLOCK DATA subprogram is used to enter initial values for variables assigned to 
labeled COMMON blocks. The BLOCK DATA subprogram is of the following form. 

BLOCK DATA 

(Specification Statements) 

END 

The first statement must be a BLOCK DATA statement, which does not allow either a 
line number or any arguments. No executable statements may appear in a BLOCK DATA 
subprogram. The entire BLOCK DATA subprogram can contain only EQUIVALENCE, DATA, 
DIMENSION, COMMON, IMPLICIT and type declaration statements. A complete set of 
specifications is given for an entire common block by a single BLOCK DATA sub- 
program, although a single BLOCK DATA subprogram can initialize any number of named 
common blocks. For example: 

BLOCK DATA 

COMMON/R/S , Y/C/Z , W , Y 

DIMENSION Y(3) ' 

COMPLEX Z 

DATA Y/1E-1,2*3E2/W,Z/11.877E0, (-1.41421,1.41421)/ 

END 

6.5 NULL ARGUMENTS IN CALLS TO SUBROUTINES OR FUNCTIONS 

An argument in a call to a subroutine or function may be left unspecified by 
writing a null, or blank, expression as the actual argument. For example, 

CALL S(A,,B) 

is a call with three actual arguments of which the second argument is null. 

The use of null arguments is only meaningful when calling routines that are 
specifically designed to accept optional arguments. Such routines can only be 
written in assembly language, not in the PDP-11 FORTRAN language. 
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In general, the use of a null argument in calling a routine that is not designed to 
detect and handle it will result in a fatal error at execution time. 

The implementation of null arguments is discussed in Chapter 7-15. 
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PART 7 

CHAPTER 7 

FORMAT STATEMENTS 

Data transmission statements govern the transfer of data between internal storage 
and either peripheral devices or other locations in the computer memory. These 
statements are divided into four categories. 

a. Data description statements: Format (described in this Chapter) and 

DEFINE FILE (Section 7-8.4.1) 

b. Input/Output statements READ, WRITE and PRINT (Chapter 7-8) 

c. Device control statements FIND, BACKSPACE, REWIND and END FILE 
for magnetic tape and Section 7-8.4.2, 7-8.6.1, 7-8.6.2, 
disk usage: and 7-8.6.3) 

d. Statements specifying data ENCODE and DECODE (Section 7-8.7) 
transmission between 

internal storage areas: 

Chapter 7-7 describes the construction of FORMAT statements . Chapter 7-8 describes 
the remaining data transmission statements . 

7.1 FORMAT STATEMENTS 

The FORMAT statement describes the format in which one or more records are trans- 
mitted. FORMAT statements are nonexecutable statements used in conjunction with 
input/output statements or with ENCODE and DECODE statements. 

FORMAT statements have the following form. 

FORMAT(S ,S.,...S /S' , S' ,....S /...) 
x z n i z n 

A statement number must be specified. (FORMAT statements are referenced by the data 
transmission statements with which they operate.) Each S is a data field speci- 
fication; the slash (/) character indicates the start of a new record. The closing 
parenthesis at the end of the FORMAT statement also terminates a record. 

Each record description may consist of one or more field specifications, a field 
being a consecutive series of characters within the record. Field specifications 
are separated by commas within the FORMAT statement (consecutive commas are not 
allowed) . 
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FORMAT statements must be placed in the executable part of the source program. Unless 
the FORMAT statement contains only alphanumeric data for direct transmission, it is 
used in conjunction with the list of a data transmission statement. The FORMAT 
statement indicates the form in which the variables specified in the data trans- 
mission statement are to be read or written. 

During data transmission, the object program scans the indicated FORMAT statement. 
If one or more specifications for a numeric field, for example, are present, and 
the data transmission statement contains numeric data items to be transmitted, those 
items are transmitted according to the format specification. Execution of the data 
transmission statement is terminated when all items in the data list have been 
handled. The FORMAT statement can contain specifications for more items than are 
contained in a given data transmission statement, in which case excess format 
specifications are ignored. Conversely, the FORMAT statement can contain speci- 
fications for fewer items than are contained in the data transmission list, in 
which case the set of format specifications beginning with the last left parenthesis 
or slash is repeated as often as necessary to supply the needed number of data 
formats. 

A field specification within a FORMAT statement is of the form 

Cw 



or 



Cw.d 



where C is a format code, w indicates the field width and d (where present) indi- 
cates the number of characters within the field that occur following the decimal 
point. The following types of format conversion codes may appear in a FORMAT 
statement. 

a. numeric: Iw, Ow, Fw.d, Ew.d, Dw.d, Gw.d 

b. numeric with scale factor: nPFw.d, nPEw.d, nPDw.d, nPGw.d 

c. logiccil: Lw 

d. alphanumeric and editing Aw, nH, '...*, nX, Tw 
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There are several ways to specify the numeric values for repeat count, (see Section 
7-7.19) field width and number of characters following the decimal point in formats. 
They are used in the following order of precedence. 

1. A default value is assumed. The default repeat count is 1. The default 
values for width and number of decimals are discussed in Section 7-7.14. 

2. If an explicit value is specified, it is used in place of the default. 
Explicit values may be specified in the format either as a literal integer 
or as a variable format expression (Section 7-7.17). 

3. In the case of formatted input, if the field thus specified contains a 
comma or extends past the end of the input record, then the field is 
said to be short field terminated and the field width (but not the num- 
ber of decimal characters) is reduced accordingly. See Section 7-7.15. 



7.2 I FORMAT CONVERSIONS 



A specification of the form 



Iw 



is used to transmit decimal integer values. The corresponding element in the data 
transmission list must be an integer or byte Variable. 

On output, a field w characters long is reserved for the value to be output. The 
value is transmitted into this field right- justified. Where a minus sign is 
expected in the output, one position must be included in w for that sign. Since 
a plus sign is suppressed, no allowance need be made for it. Where the field 
specified for the output value is too small, the entire output field is filled 
with asterisks. For example: 

Format Internal Value External Representation 

13 
14 
15 
12 
13 
17 

On input, a field of w positions is read as it appears. The external data must be 
a decimal integer and must not contain a decimal point or exponent field. Blanks 
are interpreted as zeroes. Leading blanks are ignored. An all blank field is 
interpreted as a zero on input. An integer overflow occurs if the value being 
read is beyond the range -32768 to 32767. For example: 



284 




284 


-284 




-284 


174 




AA174 


3244 




** 


-473 




*** 


29.443 


not 


permitted, error 
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Format 

14 
13 
19 
19 



External Value 

2788 
-26 

AAAAAA312 
312AAAAAA 



Internal Representation 

2788 

-26 

312 

not permitted, error; 

system attempts to input 

312000000 



7.3 FORMAT CONVERSIONS 



A specification of the form 



Ow 



is used to transmit octal integer values. The corresponding element in the data 
transmission list must be an integer or byte variable. 

On output a field w characters long is reserved for the value. The value is 
transmitted into this field right-justified. Where the field specified is too 
small, the entire output field is filled with asterisks. For example: 



Format 


Internal (Decimal) 


06 


32767 


06 


-32768 


02 


14261 


04 


27 


05 


13.52 



External Representation 

A77777 
100000 
** 

AA33 
not permitted, error 



On input, a field of w positions is read as it appears. The external values must 
not contain a sign, a decimal point, an exponent field, or the digits 8 or 9. 
Leading blanks are ignored, embedded and trailing blanks are interpreted as zeros; 
an all blank input field is treated as a value of 0. The value is read into the 
first word of the data element specified in the input statement. For example: 



Format 



External Value 



05 


32767 


04 


16234 


06 


13AAAA 


03 


97 



Internal Octal Representations 

32767 
1623 
130000 
not permitted, error 



7.4 F FORMAT CONVERSIONS 

A specification of the form 
Fw.d 
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is used to transmit real or complex values. The corresponding element in the data 
transmission list must be a real or complex variable. 

On output a field w characters long is reserved for the value , which is transmitted 
into this field right- justified. The total field length reserved (w) must be 
large enough to include a decimal point, at least one digit to the left of the 
decimal point, a sign (only minus signs are printed on output), and d decimal digits 
following the decimal point; that is, w>d+3. Where the field length specified for 
the output value is too small, the entire output field is filled with asterisks. 

For F conversion output, the data value after scaling must be within the range 

-38 38 

0.14 x 10 to 1.7 x 10 . For example: 



Format 

F8.5 

F9.3 

F2.3 

F10.4 

F5.2 

F5.2 



Internal Value 

2.3547188 

8789.7361 

51.44 

-23.24352 

325.013 

-.2 



External Representation 

A2. 35472 

A8789.736 

** 

AA-23.2435 
***** 

-0.20 



On input, a field of w positions is read as it appears. The external data must 
be a decimal number with or without a decimal point or exponent field. Blanks 
are interpreted as zeros. Leading blanks are ignored. An all blank field is 
interpreted as zero. Where no decimal point is encountered within the first w 
positions, those digits are accepted as the input value and a decimal point is 
inserted d positions from the right. A decimal point found in the input value over- 
rides the format specification. For example: 



Format 

F8.5 
F8.5 
F8.5 
F5.2 



External Value 

123456789 
12345.67 
24.77E+2 
1234567.88 



Internal Representation 

123.45678 
12345.67 
2477.0 
123 . 45 



7.5 E FORMAT CONVERSIONS 



A specification of the form 



Ew.d 



is used to transmit real or complex values, with a decimal point and E exponent 
field. The corresponding element in the data transmission list must be a real or 
complex variable. 
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On output a field w characters long is reserved for the value, which is transmitted 
into this field right- justified. The total field length reserved (w) must be large 
enough to include a decimal point, at least one digit to the left of the decimal 
point, a sign (only minus signs are printed on output), d decimal digits following 
the decimal point, and four positions for the exponent representation; that is to 
say, w>d+7 in most cases. Where the field length specified for the output value 

is too small, the entire output field is filled with asterisks. For E conversion 

38 
output, the data value after scaling must be within the range 0.14 x 10 to 

38 
1.7x10 . For example : 

Format Internal Value External Representation 

E9.2 475867 ,.22 2 A0.48EA06 

E12.3 475867 ..222 AAA0.476EA06 

E5.3 56.. 12 ***** 



On input, a field of w positions is read as it appears. The external data must be 
a decimal number with or without decimal point or exponent field. Blanks are 
interpreted as zeros. Leading blanks are ignored. An all blank field is inter- 
preted as zero. Where no decimal point is encountered within the first w positions, 
those digits are accepted as the input value and a decimal point is inserted d 
positions from the right? a zero exponent is assumed. A decimal point found in 
the input value overrides the format specification. For example: 

Format External Value Internal Representation 

734.432E8 734.432E8 

AA1022.43E-6 1022.43E-6 

52.3759663 52.3759663 

A1971.3213E12 1971.3213 

123456789 1234567. 8E0 
210.5271D+10 210.5271E10 

7.6 D FORMAT CONVERSIONS 

A specification of the form 

Dw.d 

is used to transmit real values, with a decimal point and D exponent field. The 
corresponding element in the data transmission must be a real or double-precision 
variable . 



E9.3 


E12.4 


E15.3 


E10.4 


E8.1 


E12.5 
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On output, the effect of this field specification is identical to that of an equiva- 
lent Ew.d specification, except that a D exponent indicator is used in place of the 
E indicator. For example: 



Format 

D.14.3 
D23.12 
D9.6 



Internal Value 

0.363 
5413.87625793 
1.3 



External Representation 

AAAAA0.363DA00 

AAAAA0.541387625793DA04 

********* 



On input, the effect of a Dw.d specification is identical to an equivalent Ew.d 
specification, except that the internal representation appears as in the following 
examples. 



Format 

D10.2 
D10.2 
D15.3 



External Value 

12345AAAAA 

AA123.45AAA 

367.4981763D+04 



Internal Representation 

12345000. 0D0 
123.45D0 

3.674981763D+06 



7.7 G FORMAT CONVERSIONS 



A specification of the form 



Gw.d 



is used to transmit real or complex data. 

Oh output a field w characters long is reserved for the value, which is trans- 
mitted into this field right- justified. The form of the output conversion depends 
upon the individual variable in the data transmission list. The output conversion 
format is a function of the magnitude of the data being converted, as described 
in Table 7-4. 
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Table 7-4 



Floating-Point Magnitudes and Resulting G Format Conversions 



Data Magnitude 


Effective Conversion 


m < 0.1 


Ew.d 


0.1 <_ m < 1.0 


F(w-4).d, 4X 


1.0 <_ m < 10.0 


F (w-4) . (d-1) , 4X 


10 d " 2 < m < 10 d_1 


F(w-4).l, 4X 


lO*" 1 < m < 10 d 


F(w-4) .0, 4X 


m > 10 d 


Ew.d 



(The 4X specification, which is inserted (in effect) , by the Gw.d specification, 
indicates that four blanks are to follow the numeric field representation. . The 
X specification is described in Section 7-7.16). 

The total field length reserved (w) includes a decimal point, at least one digit 
to the left of the decimal point, a sign (only minus signs are printed on output) , 
d decimal digits following the decimal point, and four positions for the exponent 
representation; that is to say, w>d+7 in most cases. Where the field length 
specified for the output value is too small, the entire output field is filled 
with asterisks. For example: 



Format 

G13.6 
G13.6 
G13.6 
G13.6 
G13.6 
G13.6 
G13.6 
G13.6 
G13.6 



Internal Value 

0.01234567 

0.12345678 

1.23456789 

12.34567890 

123.45678900 

1234.56789000 

12345.67890000 

123456.78900000 

1234567.89000000 



External Representation 

A0.123457E-01 
A0.123457AAAA 
Ml. 23457 AAAA 
AA12.3457AAAA 
AA123.457AAAA 
AA1234.57AAAA 
AA12345.7AAAA 
AA123457.AAAA 
A0.123457EA07 



For comparison, consider the following example of the same values output under the 
control of an equivalent Fw.d specification. 
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Format 



Internal Value 



External Representation 



F13.6 
F13.6 
F13.6 
F13.6 
F13.6 
F13.6 
F13.6 
F13.6 
F13.6 



0.01234567 

0.12345678 

1.23456789 

12.34567890 

123.45678900 

1234.56789000 

12345.67890000 

123456.78900000 

1234567.1 



AAAAA0. 012346 
AAAAA0. 123457 
AAAAA1. 234568 
AAAA12. 345679 
AAA123. 456789 
AA1234. 567890 
A12345. 678900 
123456.789000 
************* 



On input, a field of w positions is read as it appears. The external data must be 
a decimal number with or without a decimal point or exponent field. The rules for 
input conversion are the same as those for F format input conversion. 

7.8 COMPLEX I/O 

Complex data values are transmitted as two independent real qualities. The following 
specification consists of one repeated real specification or two successive real 
specifications of the form. 



Two Real Specifications 



Repeated Real Specifications 



Fw.d,Fw.d 
Ew.d,Ew.d 
Gw.d,Gw.d 



2Fw.d 
2Ew.d 
2Gw.d 



Nothing is output between the two parts of a complex number unless specified by 
the format. 

Output Conversion: 



Format 



Internal Value 



External Representation 



2F8.5 
E9.2,E5.3 



2.3547188,3.456732 
47587.222,56.123 



A2.35472A3. 45673 
A0.48EA06***** 



Input Conversion: 



Format 

F8.5,F8.5 
E9.1,E9.3 



External Value 

1234567812345.67 
734.432E8123456789 



Internal Representation 

123.45678,12345.67 "* 
734 . 432E8 , 1234567 . 8E-1 
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7.9 SCALE FACTORS 



Scale factors of the form 



nP 



can be specified for D, E, F, arid G conversions. The number n is a signed or 
unsigned integer specifying the scale factor. 

For E and D type output conversions, the basic real constant part of the output 
quantity is multiplied by 10 (n) and the exponent is reduced by n. For G type 
output conversion, the effect of the scale factor is suspended unless the magnitude 
of the value to be converted is outside the F conversion range (see Table 7-4) . 

For the F, E, G and D type input conversions where no exponent is found in the 
external value, and F type output conversions, the scale factor causes the external 
value to equal the internal value times 10 (n) . 

For F, E, G and D type input conversions, the scale factor has no effect if there 
is an exponent in the external value. 

Where no scale factor is specified, a scale factor of zero is assumed. Once a 
scale factor has been specified, it is carried onto all subsequent D, E, F, and G 
type conversions within the same format specification unless another scale factor is 
encountered. A zero scale factor can only be reinstated by an explicit 0P speci- 
fication. 

Output Conversions: 



Format 

3PE15.3 

1PE15.3 

1PE15.2 

-1PE15.2 



External Value 

-273.612 
-273.612 
-273.612 
-273.612 



External Representation 

AM-273.612EA00 
AAAAA-2.736EA02 
AAAAAA-2.74EA02 
AAAAAA-0.03EA04 



Input Conversions; 



Format 

3PE10.5 

3PE10.5 

-3PE10.5 



External Value 

AAA37.614 

AA37.614E2 

AAAAA37.614 



Internal Representation 

.037614 
3761.4 
37614.0 
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7.10 L FORMAT CONVERSIONS 
A specification of the form 

Lw 

is used to transmit the value of a logical variable. 

On output, a field w characters long is reserved for the value, which is trans- 
mitted into this field right-justified. The letter T or F is transmitted following 
w-1 blanks. For example: 

Format Internal Value External Representation 

L5 - .TRUE. AAAAT 

LI .FALSE. F 

On input, the first nonblank character in the data field must be a T or an F; the 
value of the logical variable is stored as true or false, respectively. If the 
data field is blank or empty, a value of false is stored. Any other value in the 
data field causes an error diagnostic. 

7.11 A FORMAT CONVERSIONS 

Alphanumeric data can be transmitted in a manner similar to numeric data by use of 
a specification of the following form. 

Aw 

Conversions of this type are generally used where the alphanumeric data is generated 
or operated upon by the user program. 

The field specified has a length of w characters. The alphanumeric characters are 
transmitted as the value of a variable in the data transmission list. The variable 
may be of any type. The value of w is limited to the maximum number of characters 
that can be stored in the space allotted for a single variable of the type speci- 
fied. Table 7-5 indicates the number of alphanumeric characters that can be stored 
in a single variable of a given type. 
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Table 7-5 



Alphanumeric Data Storage 





Number of Characters 


Variable Type 


per Variable 


Integer 


2 


Byte 


1 


Logical 


2 


Real 


4 


Doubl e-pr ec i s ion 


8 


Complex 


8 



For example, the sequence 

READ(2,5) V 
5 FORMAT (A4) . 

causes four characters to be read and stored in memory as the value of the real 
variable V. 

If the value of w exceeds the amount specified in Table 7-5, the leftmost characters 
are lost on input; on output, w characters appear right- justified in the external 
output field. 

If the value of w is less than the number of characters that can be stored in the 
space allotted to a variable, on input the characters are left- justified and 
blank-filled on the right of each list item. On output the leftmost w characters 
in the variable are transmitted to the output field. 

7.12 ALPHANUMERIC DATA WITHIN FORMAT SPECIFICATIONS 

Alphanumeric data can be transmitted directly into or from a format specification 
by means of Hollerith conversion (H) or of single quote delimiting characters. 
Direct transmission of alphanumeric data is performed where such data is not 
operated upon by the program; i.e., headings and explanations to be printed along 
with output data. 
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The H conversion is indicated by a specification of the form 
nH 

where n indicates the number of characters to be transmitted. On output, the follow- 
ing n characters are sent to the output device. On input, n characters are read 
and inserted into the Hollerith conversion field of the format specification. For 
example: 

C OUTPUT HEADING TO LINE PRINTER: 

WRITE (5,10) N 
10 FORMAT (15H OUTPUT HEADING, 13) t 

C INPUT HEADING FROM KEYBOARD FOR LATER USE; UP TO 30 CHARS: 

READ (6,20) 
20 FORMAT (30H ) 



In line 10, the leading space in the line to be output is the carriage control 
character (see Section 7-7.14). The 14-character heading is output, followed by 
the integer value of N printed in a field of three characters. The comma following 
the Hollerith field is optional. 

In line 20, the Hollerith field indicated is blank; however, any alphanumeric 
characters within this field would be replaced by those characters read from the 
keyboard. If fewer than n characters are read, they overlay the first portion of 
the Hollerith field , with the remainder of the field being blank filled. 

An alternative form is to enclose Hollerith data in single quotes (apostrophes) . 
The result on both input and output is exactly as described for H conversion. An 
apostrophe within such a Hollerith string must be represented by two successive 
single quotes to distinguish it from the terminating single quote of the string. 
For example: 

FORMAT ('DON' ' 'T TYPE ON KEYBOARD UNTIL OUTPUT IS FINISHED') 
When output is to the keyboard or line printer, the line 

DON'T TYPE ON KEYBOARD UNTIL OUTPUT IS FINISHED 

is printed. 

Note that variable format expressions (Section 7-7.17) cannot be used with H 
conversion. 
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7.13 Q FORMAT SPECIFICATION 

The Q format specification is used to read the actual number of characters input 
on the current READ operation. It must correspond to an input list item of 
integer type into which the character count will be placed. For example, 

READ(4,10)X,Y,Z,K 
10 FORMAT (3F13. 5, Q) 

reads input values into the real variables X, Y, and Z, and places the number of 
characters read in the integer variable K. 

Q format is ignored on output. 

7.14 DEFAULT FIELD SPECIFICATIONS 

The w.d may be omitted from the numeric and logical conversions, in which case 
the following default values are used. 

Format Default 



E 15.7 

F 15.7 

G 15.7 

t> 25.16 

I 7 

7 

L 2 

A 8 



For example, the format 

FORMAT (I, 2E, A) 
is equivalent to the format 

FORMAT (17 , 2E15 . 7 ,A8) 

See also the discussion on short field termination (Section 7-7.15). 

Note that in the case of E, F, G and D formats, the entire w.d portion of the format 
is defaulted; partial defaulting is not allowed. 

Legal defaults: E, F 

Illegal defaults: E10, F.2, F. 
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7.15 SHORT FIELD TERMINATION ON FORMATTED INPUT 

When using formatted input, either a comma in an input field or the end of the 
input record will override the field width specified by the FORMAT statement. This 
provides a simple method of free form input that is particularly useful when 
entering data from a keyboard, short field termination may be used with E, F, G, 
I, L, and formats. For example, 

READ (6,100) I, J, A,B 
100 FORMAT (216, 2F10. 2) 

and input data record 

1,-2,1.0,35 

v 
will result in the following assignments being performed by the READ statement. 

1 = 1 

J - -2 
A = 1.0 
B = 0.35 

Note that the comma must be included within the field width specified by the FORMAT 
statement. The only effect of the comma is to adjust the field width used by the 
conversion routines. A zero length field is allowed. In the case of a null 
(zero length field) the value read will be 0. for D, E, F or G formats, and I or 
formats and .FALSE, for L format. 

A comma may not be used to terminate a field being read with A format. However, if 
the end of record occurs within the A field width, the available characters will be 
transmitted successfully with trailing blanks appended as necessary to fill the 
data element. 

7.16 RECORD LAYOUT SPECIFICATIONS 

Spacing control within a record is performed through the use of the X specification 
to generate a given number of blanks or skip a given number of characters, or the 
T specification which allows the user to specify "tabulation" to a given position 
within a record. 

A specification of the form 
wX 
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provides for the transmission of w spaces on output or the skipping of w characters 
on input. The value of w must be greater than zero, For example: 

WRITE(5,75) I,X 
7 5 FORMAT ( 5H STEPI 5 , 10X 2HY«F7 . 3 ) . 

could be used to output the folloiwng line (A indicates a blank) . 

STEP AAA28AAAAAAAAAAY=b- 3 . 8 7 2 

The first blank would not be printed on the terminal or line printer as it would be 
considered the carriage control character. Notice that the comma following an X 
specification is optional. 

A specification of the form 

Tw 

specifies the character position in a FORTRAN record where the data transfer is to begin. 
The value of w must be an unsigned integer constant and indicates the (w)th character 
of the record on input and on nonprinted output, and the (w-l)th character of a printed 
output record (where the first character is a carriage control character). For example, 

2 FORMAT (T50,'BLACK/T30, 'WHITE 'J 
would cause the following line to be printed. 

Print Position 29 Print Position 49 

WHITE BLACK 

The statements 

1 FORMAT (T35,"MONTH') 
READ (3,1) 

causes the first 34 characters of the input data to be skipped, and the next five 
characters would replace the characters M, 0, N, T, and H in storage. If an input 
record containing 

ABCJrftftfXYZ 
is read with the format specification 

10 FORMAT (T7,A3,T1,A3) 
then the characters XYZ and ABC are read, in that order. 

7.17 VARIABLE FORMAT EXPRESSIONS 

An integer valued expression may be used in a FORMAT statement any place an integer 
may be used (except as a Hollerith count) by enclosing it in angle brackets, 
<expression>. For example, 

FORMAT (I<J+1>) 

will cause an I conversion with a field width one greater than the value of J at the 

time the format is scanned. The expression will be re-evaluated each time it 
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is encountered in the normal format scan. If the expression is not of type integer 
it will be automatically converted to integer prior to use. Any valid FORTRAN 
expression may be used including function calls and references to dummy parameters. 

A complete example is shown in Figure 7-4. 

The value Of a variable format expression must be valid for its use in the format 
or an error will occur at runtime. 



The FORTRAN program 

DIMENSION A (5) 
DATA A/1. ,2. ,3. ,4., 5./ 
C 

DO 10 I = 1,10 

Write (5,100) i 

100 FORMAT (I<MAX0 (1,5) >) 
10 CONTINUE 

C 

DO 20 I =1,5 

WRITE (5,101) (A(J),J=1,I) 

101 FORMAT (<I>F10.<I-1>) 
20 CONTINUE 

END 

produces the following output when executed. 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
1. 

2.0 2.0 
3.00 3.00 3.00 
4.000 4.000 4.000 4.000 
5.0000 5.00000 5.0000 5.0000 5.0000 



Figure 7-4 
Variable Format Expression Example 
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7.18 CARRIAGE CONTROL 

The first character of each ASCII record controls the vertical spacing of the line 
printer and terminal. When an output line (record) is generated, it is constructed 
within the system exactly as described by the combined WRITE and FORMAT statement 
designations. Before output to either of the printing devices, the first character 
of the record is interpreted as the carriage control character and is deleted from 
the output line. 

The carriage control characters are specified in Table 7-6 

Table 7-6 
Carrieige Control Characters 



Character 


Effect 


A 


blank 


carriage return/line feed 

(advance printing position to beginning of next line) 





zero 


carriage return/line feed/line feed 

(advance printing position to beginning of the line 

after the next) 


1 


one 


carriage return/form feed 

(advance printing position to the beginning of the next 

line and perform a form feed where the hardware capability 

exists) 


+ 


plus 


carriage return 

(advance printing position to the beginning of the current 

line, allows current line to be overprinted) 


$ 


dollar 


inhibit carrieige return/line feed 

(do not advance printing position, used when creating output 

lines to interact with the keyboard) 



Any character other than those described in Table 7-6 that appears first in the line 
is deleted from the output line and treated as though it were a blank carriage 
control character. 
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7.19 REPETITION OF FIELDS, GROUPS, AND MULTIPLE RECORDS 

A field specification can be repeated by preceding that specification with an 
unsigned integer constant denoting the number of repetitions. For example, 

2BT FORMAT (2E12. 4, 315) 

is equivalent to 

20 F0RMAT(E12.4,E12.4,I5,I5,I5) " 

Similarly, a group of field specifications can be repeated by enclosing the group 
in parentheses and preceding the whole with an unsigned integer denoting the 
number of repetitions. For example, 

3jZf FORMAT(2I8,2(El5.5,2F8.3)) 

is equivalent to 

3/ F0RMAT(2I8,E15.5,F8.3,F8.3,E15.5,F8.3,F8.3) 

Where several input/output records are specified in a single format specification, 
and those different records have different field specifications, a slash is used 
to indicate a new record. For example, 

WRITE(6,40)K,L,M,N,O,P 

40 FORMAT (308/16, 2F 8. 4) 

is equivalent to the following. * 

WRITE(6,40)K,L,M 
WRITE(6,40)N,O,P 
4j6 , FORMAT (308) 

41 FORMAT (16, 2F8. 4) 

A single record may be up to 133 characters long; this provides for a 132-character 
line printer plus a carriage control character. 

The separating comma can be omitted when a slash is used. The presence of n slashes 
at the beginning, or end of a FORMAT statement causes n blank records to be output 
or to be skipped on input. When n slashes appear between the first and last 
FORMAT specifications, (n-1) blank records are output or skipped on input. 
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7.20 FORMATS STORED AS DATA 

The ASCII character string that makes up a format specification may be stored as 
the values of an array. Input/output statements may refer to the format by 
giving the array name instead of the statement number of a FORMAT statement. The 
stored format has the same form as a FORMAT statement, except for the word 
FORMAT and the statement number. The enclosing parentheses are included, however. 

As an example, consider the following sequence. 

DIMENSION SKEL(2) 
READ(4,1)SKEL 
1 FORMAT (2A4) 

READ(5,SKEL)K,X 

The first READ statement enters the ASCII string into the array SKEL, In the 
second READ statement, SKEL is referenced as the format-governing conversion of 
K and X. 

As another example, the format created by the DATA statement below 

REAL FORM1 (3) 

DATA FORfil/ 1 (16' , V/I6' , ' ,13) '/ 

WRITE(6,F0RM1) I,J,K 

is equivalent to 



WRITE (6, 50) I,J,K 
50 FORMAT (16 //I6,I3) 



One caution should be noted with respect to the storing of formats within arrays. 
Such storage does not result in the correct linkage of the I/O conversion 
routines from the OTS Library. For example: 



REAL FORM1 (2) 

INTEGER II 

LOGICAL LL 

DATA F0RM1 / ' (15, ' , 'L5) '/ 

11=0 

LL=.TRUE. 

WRITE(6,F0RM1) II, LL 

END 



If the preceding program is run, it produces the following run- time error 
message. 

FORT008000 LINKAGE ERROR (MISSING FORMAT CONVERSION ROUTINE) 
When this error occurs, execution of the program is terminated. 
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To correct this error, the programmer should include within his FORTRAN program a 
FORMAT statement that contains any format conversions that will be needed at run- 
time. This FORMAT statement need not be used or referenced by any READ or WRITE 
statement; it need only be present. In the preceding example, the following 
statement could be included immediately before the END statement. 

9999 FORMAT (II, LI), 

Also note that the variable format expression notation cannot occur in formats 
stored as data. 
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PART 7 

CHAPTER 8 

I/O AND DEVICE HANDLING STATEMENTS 

8.1 INTRODUCTION TO FORTRAN I/O 

Input of data to a FORTRAN program is performed by the READ statement. Output from 
a FORTRAN program is performed by either the WRITE or PRINT statement. These 
statements are generally used in conjunction with a FORMAT statement or format 
specification contained in an array. Each READ or WRITE statement contains a 
reference to the device to or from which the data transfer is to occur. The 
PRINT statement assumes output to the line printer, although the assumed device can 
be altered. 

There are three categories of READ and WRITE statements. 

a. Formatted READ, WRITE and PRINT statements are used with or without 
an I/O list. Where a list of variables is provided, those values are 
transferred and automatically converted to ASCII on output or internal 
format on input. Where an I/O list is not provided, an alphanumeric 
data transfer is indicated. 

b. Unformatted READ or WRITE statements perform a transfer of binary 
information only. These statements make no reference to FORMAT 
statements or specifications. 

c. Direct access disk READ and WRITE statements perform I/O of fixed-length 
records to a file previously defined by a DEFINE FILE statement. Direct 
access READ and WRITE statements cannot be used with FORMAT statements 
or specifications. 

A feature available with all forms of READ and WRITE statements allows the program 
to transfer control to a specified line number if an end-of-f ile or error condition 
develops during the I/O operation. 

Table 7-7 shows the skeletal formats for each of the three categories described 
above . 

The remaining statements covered in this chapter perform file-oriented or device- 
oriented operations with respect to bulk storage devices such as magnetic tape, 
disk and DEC tape. 
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Table 7-7 
READ and WRITE Statement Summary 



Category 


Input Statements 


Output Statements 


formatted 


READ (u,f) list 


WRITE (u,f) list- 




READ (u,f) 


WRITE (u,f) 




READ f, list 


PRINT f, list 


unformatted 


READ (u) list 


WRITE (u) list 




READ (u) 


WRITE (u) 


direct disk access 


READ (u'r) list 


WRITE (u'r) list 


transfer on END= and 


READ (u,f,END=n) list 


WRITE (u,f,END=n) list 


ERR= can also be 


READ (u,f,ERR=n) list 


WRITE (u,f ,ERR=n) list 


used with unformatted 


READ (u f f ,END=n,ERR=n) 


WRITE (u,f,END=n,ERR=n) list 


and direct access I/O 


list 




statements 






Where: 


u is a logical unit number, either an integer constant or variable 


f is a format reference (line number of a FORMAT statement or an array 


name) 


list is a list of I/O variables 


r is an integer variable or constant designating the record number 


n is a statement number 


* is the character used to separate the symbolic file number and the 


record number 


, is the character used to separate all other arguments of I/O statements 



8.1.1 I/O Devices 

PDP-11 FORTRAN supports all I/O devices supported by the Monitor. The default 
FORTRAN logical device assignments are described in Table 7-8. Device assignments 
•can be altered by means of the ASSIGN routine. (See Section 7-14.2.3.) 
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Table 7-8 
FORTRAN Logical Device Assignments 



Logical Unit Number 


Device 


1 
2 
3 
4 
"5 
6 
7 
8 


system disk, SY: 

system disk, SY: 

system disk, SY: 

high-speed paper tape reader, PR: 

line printer, LP: 

terminal, KB: 

system disk, SY: 

BATCH system input device, BI: 



8.1.2 I/O Records 

All data is transmitted by input/output statements in terms of records. The maximum 
amount of information in one record and the manner of separation between records 
depends upon the medium. For punched cards, each card constitutes one record; on a 
teletypewriter a record is one line; for ASCII records, the amount of information is 
specified by the FORMAT reference and the I/O list; for unformatted binary records, 
the amount of information is specified by the I/O list; DEFINE FILE is used to speci- 
fy record length for direct-access files. 

Each execution of an input or output statement initiates the transmission of a 
new data record. Thus, the statement 

READ 2,FIRST,,SEC0ND„THIRD 

is not necessarily equivalent to the following statements. 



READ 2, FIRST 
READ 2, SECOND 
READ 2, THIRD 



In this case at least three separate records are required, whereas the single 
statement 

READ 2, FIRST, SECOND, THIRD 
may require one, two, three, or more records depending upon FORMAT statement 2. 



7-88 



If an input/output statement requests less than a full record of information, the 
unrequested part of the record is lost and cannot be recovered by another input/ 
output statement without repositioning the record. Repositioning is not possible 
on all devices (see Section 7-8.6). If an input list requires more than one ASCII 
record of information, successive records are read. 

8.1.3 I/O Lists 

An input/output list contains the names of variables and array elements whose 
values are assigned on input or written on output. Input and output lists have 
the same form, consisting of individual elements separated by commas. 

The simplest form of I/O list consists of a list of variable names or array element 
names. For example, 

A,B,C 

indicates three simple variable names. 

An array name used alone in an I/O list specifies all of the array elements defined 
when the array is declared. Elements are read into core or output in the order 
specified for array storage (see Section 7-5.3.1). Thus, where the maximum element 
of the array X is X (400) , an I/O list to include the entire array X could be 
written as 

A, B ,C ,X 

The elements of a list are specified in the order of their occurrence from left to 
right and may be enclosed in parentheses. The following list is identical in 
meaning to the previous list. 

(A,B,C),X 

An implied DO loop within the I/O list can be used to indicate transmission of only 
part of an array. This is performed by including a parenthesized quantity as a 
list element as follows: 

One-Dimensional Array Two-Dimensional Array 

(A(I) ,I-n,m,k) ((A (I, J) ,I»n,m,k) ,J«n' ,m' ,k' ) 



or 



(A(l),I«n,m) ((ft(I,J),I«n',i) ,J*n l ,iti') 
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where A is an array name, I and J are array subscripts, and n, m, and k represent 
the initial, terminal and step values of the implied loop. When the step parameter 
is omitted it is presumed to be +1. For example, 

(X(K),K=1,4) 
is equivalent to 

X(1),X(2),X(3),X(4) 
The implied DO loop 

((Y(I,J) f I-l f 2),J-l,3) 
is equivalent to 

Y(1,1),Y(2,1) ,Y(1,2),Y(2,2),Y(1,3) ,Y(2,3) 

As shown above, nesting is used to provide for various levels of indexing, up to three 
dimensions . 

8.2 FORMATTED I/O 

8.2.1 Formatted BEAD Statement 

A formatted READ statement is of the form 

READ (u,f) list 
READ (u,f) 



READ f, list [equivalent to READ (4,f)] 

where u indicates a logical unit number and f indicates a format reference (either 
a FORMAT statement or an array name whose elements contain a format specification) . 

A formatted READ statement causes information to be read from the specified device 
and placed in memory. The data is converted from external to internal format as 
Indicated by the format specification. If an I/O list is provided, data items are 
stored as the values of the listed variables. If a list is not provided the data is 
stored as an alphanumeric string in the format itself. For example, 
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READ (1,100) A,B,C 

reads three values from device 1 according to the format in the FORMAT statement 
at line 100 and assigns the values to the three variables named. 

READ (1,100) 
100 FORMAT (5HHELL0) 

When the above READ statement is executed, and the characters ABODE (for example) 
are given as input, the FORMAT statement at line 100 becomes 

100 FORMAT (5HABCDE) 

A READ statement of the form ' 

READ 110,X,Y,Z 

is similar to the first form discussed, except that input is assumed to come from 
logical unit 4 (default to PR:). The FORMAT statement associated with the operation 
is at line 110. 

Example: 

DIMENSION B(32) 
50 READ (6,200) B 



100 READ (6,B) I,J,K,L 



200 FORMAT (32A4) 

When statement 50 is executed in this example, the value of the array B is read 
from device 6. If the first two elements of B contain the following 

B(l) - '(4I6 1 
B(2) = •)• 

and the remaining elements of B, B(3) through B(32) , are blank, the array B contains 
the format specification (416) . (The format statement in line 200 provides for- 
matting information for the user to input enough data to completely fill the array B.) 
When statement 100 is executed, the contents of the array B become the format 
reference for that READ statement. The variables I, J, K and L are read from 
device 6 according to the format 416 (that is, four 6-digit integers). This type of 
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READ statement (statement 100) is said to contain variable read-in formats and 
allows the user to specify an I/O format at run time (see Section 7-7.16). 

Any format conversions and I/O control routines not needed in the resident section 
but required by overlay sections must be forcibly loaded into the resident section. 
This can be done by declaring the appropriate globals in an assembly language 
routine or inserting dummy FORMAT and input/output statements in the resident main 
program for all those routines needed in the overlays and not required in the 
resident- section. See Section 7-9.10 for further details. 

8.2.2 Formatted WRITE Statement 

A formatted WRITE statement is of the form 

WRITE (u,f) list 



WRITE (u,f) 

where u indicates a logical unit number and f indicates a format reference. 

A formatted WRITE statement causes one or more logical records to be constructed 
and written on the designated device in ASCII format. A logical record contains 
up to 133 characters, including the carriage control character. The data is 
converted to external form as specified in the designated format reference. If 
an I/O list is provided, the values of the variables indicated are output. If a 
list is not provided, the information is read directly from the format reference 
and written on the device designated in ASCII form (this is generally done with 
Hollerith output) . For example: 

WRITE (5,115) IERR,INET,ILOW 
115 FORMAT (314) 

The above combination outputs three integer values to device 5. 

WRITE (5,120) 
120 FORMAT (25H THIS IS HOLLERITH OUTPUT) 

The above combination outputs 25 characters to form, the message 
THIS IS HOLLERITH OUTPUT 



7-92 



on device 5 (the first character, blank, is not printed, but is used as a carriage 
control character) . 

8.2.3 PRINT Statement 

The PRINT statement is of the form 

PRINT f, list 

where f indicates a format reference. The PRINT statement is a formatted WRITE 
statement whose output is always sent to the same device, the line printer. PRINT 
performs a WRITE (5,f). The default device can be changed by using the ASSIGN 
subroutine (see Section 7-14.2.4) or the ASSIGN Monitor command with respect to 
logical unit 5. 

As an example, 

PRINT 115, NA,NB,NC 

outputs three values according to the format described in line 115. This statement 
is analogous to the third type of READ statement described in Section 7-8.2.1. 

8.3 UNFORMATTED I/O 

8.3.1 Unformatted READ Statement 

An unformatted READ statement is of the form 

READ (u) list 
or 

READ (u) 

where u indicates a logical unit number. 

An unformatted READ statement causes binary information to be read from the unit 
designated and stored in memory as the values of variables in the I/O list, if any. 
For example: 

READ (7) A,B,C 
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This statement reads enough binary values from device 7 to fill three real 
variables (2 words each). If A, B, and C have been previously defined as double- 
precision variables, then three real values (4 words each) are read, and so on. 

READ (7) 

allows the user to skip over one record on device 7. One such statement must 
appear for each record to be skipped. 

Each READ statement reads one unformatted binary record. If the I/O list does not 
use all values in the record, the remaining values are discarded. If more values are 
required than are contained in the record, a run-time error message results. 

Unformatted binary I/O generally is used to enable the storage of intermediate 
results from one program to another without the loss of precision involved in for- 
matting data output. Unformatted I/O is performed between memory and such peri- 
pherals as paper tape, DECtape, magnetic tape or disk, but hot to an ASCII device 
where such I/O is generally meaningless. 

8.3.2 Unformatted WRITE Statement 

An unformatted WRITE statement is of the form 

WRITE (u) list 



WRITE (u) 

where u indicates a logical unit number. 

An unformatted WRITE statement causes binary information to be read from memory 
and written on the device indicated in binary form. For example: 

WRITE (7) A,B,C 

This statement writes the binary values of the variables A, B and C onto device 7. 

WRITE (7) 

allows the user to write one null record on device 7. Execution of such a state- 
ment must appear for each null record to be output. 
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8.4 DIRECT ACCESS I/O 

Direct access I/O can be performed on any directory-structured device. 

8.4.1 DEFINE FILE Statement 

The DEFINE FILE statement defines the record structure of a disk or DEC tape file. 
However, DEFINE FILE automatically allocates a file if none exists. The DEFINE FILE 
statement can be used as either a declarative (nonexecutable) or executable state- 
ment in the program. Where DEFINE FILE is executable, arguments can be computed 
within the program. 

The DEFINE FILE statement is formatted as follows : 

DEFINE FILE u (m,n,U,ivar) 

where the arguments are as follows: 

u is an integer constant or variable used as the logical unit number 
for this file specification. 

m is an integer constant or variable defining the number of records 
in the disk file . The records in the file are numbered from 
1 to m. 

n is an integer constant or variable defining the length (in words) 
of each file record. 

U is a fixed argument designating that the disk file is unformatted 
(i.e., binary). No other characters are legal. 

ivar is an integer variable, called the associated variable, which is set 

to point to the next record at the conclusion of an I/O operation on the 
file. If the associated variable is to be used by more than one 
subprogram or in an overlay environment, it should be placed in a resi- 
dent common block. 

The statement 

DEFINE FILE 1 (1000, 100, U, IVAR) 

specifies a 1000-record file on unit number 1, each record of which is 100 words 
long. 

The variable IVAR maintains an index of records processed, providing a pointer to the 
next record. 
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With respect to overlays, the DEFINE FILE routine uses storage space that must be 
allotted in such a way that the overlay system is not corrupted. In the statement 

DEFINE FILE 1(10,10,U, INDX) 

the variable INDX (the associated variable) is used to point to the record to be 
processed for statements such as the following ones. 

READ (1* INDX) 
WRITE (l'INDX) 
FIND (l'INDX) 

Every READ, WRITE or FIND done on unit 1 modifies INDX. INDX must be placed in the 
local storage of the resident section or in a resident common block. If it is not, 
the FORTRAN I/O system may destroy an unpredictable location in core when attempting 
to store into INDX. 

8.4.2 FIND Statement 

The format of the FIND statement is 
FIND (u'r) 

where u is the file specification as assigned by the DEFINE FILE statement, and r is 
an integer variable or constant specifying the record number. The FIND statement 
is included for compatibility with other FORTRAN systems. The Monitor does not 
support disk head positioning; therefore, the only effect of executing this statement 
is to update the associated variable for the specified unit. 

8.4.3 Direct Access READ Statement 

A direct access READ statement is of the form 

READ (u'r) list 

where u is a logical unit number as established in a prior DEFINE FILE statement 
and r is an integer variable or constant specifying the record number. An 
apostrophe is used to separate the two elements to differentiate this form of 
READ statement from a formatted READ statement. The list specifies the number 
of the items to be read. 
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A direct access READ statement provides random access to fixed-length records in a 
disk or DECtape file. For example, 

READ il'50) ADC, ADX, AOD 

reads record number 50 of the file on unit number 1 to obtain values for the 
variables ADC, ADX and AOD. Three 2-word values are transmitted since the variables 
indicated are real variables. If the variables were byte, integer, complex, or 
double precision, a corresponding number of bytes or words of data would be read for 
each variable. 

As in other I/O operations, only the number of items in the list are transmitted 
regardless of the number of items in the record. If a list requires more values 
than are provided in the record, an error message is printed. 

Implied DO loops can be included in direct access I/O statements. 

8.4.4 Direct Access WRITE Statement 

A direct access WRITE statement is of the form 

WRITE (u'r) list 

where u is a logical unit number and r is an integer variable or constant specifying 
the record number. An apostrophe is used to separate the two elements to differentiate 
this form of WRITE statement from a formatted WRITE statement. The list specifies 
the number of items to be output. 

The direct access WRITE statement outputs a fixed-length record directed into a 
file and is analogous to the direct access READ statement described in the 
previous section. 

8.5 TRANSFER OF CONTROL ON ERROR CONDITION 

Both READ and WRITE statements can be written so as to transfer control to a given 
statement number if an error condition is detected. These variations can be 
specified in one of the following formats: 

READ(u,f,END=n) list 
READ(u,f ,ERR=n) list 
READ(u,f,END=,ERR=n) list 
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WRITE (u,f,END=n) list 
WRITE (u,f,ERR=n) list 
WRITE (u,f,END=n,ERR=n) list 

where u is a logical unit designation, f is a format reference, and n is a line number 
within the program. The list is optional. END= and ERR= can also be used with 
unformatted and direct-access I/O. 

The arguments END=n and ERR=n can appear separately, or together in the order 
illustrated. If an end-of-file condition is encountered during a READ, control 
transfers to the statement having the line number associated with the END parameter. 
If an END parameter, is not specified, I/O on that device terminates and the program 
terminates with an error indication. 

If an error is encountered on input or output, control transfers to the statement 
having the line number associated with the ERR parameter. If an ERR parameter. is 
not specified, the program terminates with an error message. 

If an error variable has been associated with a logical unit as a result of a 
CALL ASSIGN or CALL SETFIL statement, then when the ERR transfer is taken, the 
error variable will contain a code designating the type of the error. This code 
is the number of the error message within the Class 1 errors as described in 
Appendix K. 

8.6 • DEVICE CONTROL STATEMENTS 

Refer to Table 7-13 which defines the type of I/O and the specific operations that 
can be performed on each of the peripherals. 

8.6.1 BACKSPACE Statement 

The BACKSPACE statement is of the form 

BACKSPACE u 

where u is an unsigned integer constant or integer variable specifying the logical 
unit number of a device. The unit is repositioned to the beginning of the file 
opened on that unit, and the read/write mechanism spaces forward n-1 records, where 
n is the number of the record processed before the BACKSPACE. The effect is to 
backspace one ASCII record or one logical binary record, depending upon the device. 
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8.6.2 REWIND Statement 

The REWIND statement is of the form 

REWIND u 

where U is an unsigned integer constant or integer variable specifying the logical 
unit number of a device. This statement repositions the designated unit to the 
beginning of the currently opened file. 

8.6.3 END FILE Statement 

The END FILE statement is of the form 

END FILE u 

where u is an unsigned integer constant or integer variable specifying the logical 
unit number of a device on which input or output operations are being performed. 
The END FILE statement causes the currently open file to be closed. END FILE also 
closes any open random-access files to allow another DEFINE FILE to be performed. 

8.7 ENCODE AND DECODE 

ENCODE and DECODE statements transfer data, according to format specifications, 
from one section of memory to another. DECODE changes data in ASCII format to the 
specified format. ENCODE changes data of the specified format into ASCII format. 
The two statements are of the form 

ENCODE (c,f,v)list 
DECODE (c,f,v)list 

where the arguments are as follows: 

c the number of characters (bytes) to be converted. 

f the FORMAT statement number or array specifying the format. 

v the name of the array that contains, or is to contain, the 
ASCII character string? it must be an array name. 

list the list of variables to be changed from or into ASCII format; it is 
the same as any other I/O list (see Section 7-8.1.3). 
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Example of the use of the DECODE statement: 



DIMENSION A (3) ,1(3) 
DATA A/«1234 , , , 5678 , , , 9012'/ 
DECODE (12, 100, A) I 
100 FORMAT (314) 



This code causes the integer array I to take on the numeric values of the three 
strings stored in the ASCII array, as follows. 

1(1) = 1234 
1(2) = 5678 
1(3) = 9012 

DECODE is analogous to a READ statement, since it causes conversion from ASCII to 
internal format. ENCODE is analogous to a WRITE statement, causing conversion 
from internal ..data formats to ASCII. 

The DECODE statement is used when it is necessary to process records having different 
formats. Figure 7-5 contains a segment of FORTRAN code that reads 80-character 
records from logical unit 1. The first two characters of each record contain an 
integer between 1 and 3 indicating the format of the remainder of the record. 
DECODE is used twice in the example: first to scan the record type field and then 
to read and translate the remainder of the record. 

Note that the character data processed by ENCODE/DECODE must be in adjacent bytes 
in core. This means that if character data is to be stored in an integer array, 
the /ON (1-word integers) compile time switch must be used. If 1-word integers 
are undesirable then the character data should be stored in a real array. 
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BYTE BUFF(80),TITLE<12> 
C 

C GET NExT RECORD 
5 READ(1,10,END*900) BUFF 

10 FORMAT(80A1> 

C 

C DECODE T Y PE FIELD AND BRANCH ON IT 

C 

DECODE<80,20»BUFF> IT Y P 
20 F0RMATU2) 

IF <ITYP,GE,1,AND,ITYP,LE„3) GOTO (100,200,300) ITYp 
C 
C 3AQ RECORD 

WRiTE(5,33)BUFF 

30 FORMAT (' #*BAD RECORD**' /1X, 80A1) 

GO TO 5 
C " 

C lTyP«l 

100 QECODE<80,U0,BUFFn»J»K 
110 F0RMAT(T4 l 2I4,T24i!6> 



C 

C ITYP*2 

200 DECODE < 80, 210, BUFF) A, B,C,D, I, J 

210 FORMAT«T20|4F10 l 5i2i5> 



C 

C I T Y P ■ » 3 

300 DECODE (80, 310, BUFF) TITLE, I, J 

310 FORMAT(T3,12A1,T60,2I5) 



C ENb"OF«FILE TRANSFE RS T 900 

900 CONTINyE 



Figure 7-5 
Decode Example 
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PART 7 

CHAPTER 9 

OPERATING PROCEDURES 

9.1 USING THE FORTRAN SYSTEM 

Figure 7-6 outlines the steps required to prepare a FORTRAN source program for 
execution under the DOS /BATCH Monitor: (1) compilation, (2) linking, and 
(3) execution. 



(OTS *\ 

Library J 




Compiler 



/^Object \ 
\^ Module/ 



Linker 



(Load A f Running ^ 
Module _J \^ Program J 















Figure 7-6 
Steps in Compiling and Executing a FORTRAN Program 



Step 1 in Figure 7-6 is initiated by a call to the FORTRAN Compiler, accompanied by 
a command string that describes input and output files, and switch options, if desired, 
to be used by the compiler. Step 2 is initiated by a call to the linker, accompanied 
by a similar command string. Steps 2 and 3 can be made to follow step 1 automatically 
by including the /GO switch option in the command string presented to the compiler. 
The /GO switch causes the linker to be called, the program linked, and execution 
begun in automatic sequence following compilation. 
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9.1.1 Filename Specifications 

The DOS/BATCH FORTRAN compiler accepts a standard DOS command string of the follow- 
ing form. 

output files < input files 

The DOS/BATCH command string syntax is described more completely in Part 3 the 
DOS/BATCH Monitor. 

The FORTRAN Compiler can produce two output files: an object module file and a 
listing file. The single input file permitted is a FORTRAN source file, containing 
one or more FORTRAN program units (a FORTRAN subprogram or main program) . Each 
file specification is of the form 

dev : filename . ext [uic] 

where dev: is any legal device specification code. The usual DOS device specifica- 
tions are shown in Table 7-9. 

Table 7-9 
Device Specification Codes 



dev: 


Device 


SY: 


system disk, assumed default device 


DKn: 


RK11 disk, cartridge unit n; n=0 is the default condition 


DC: 


RC11 disk 


DF: 


RFll disk 


DPn: 


RP03 disk, cartridge unit n; 
n=0 is the default condition 


KB: 


user terminal 


CR: 


card reader 


LP: 


line printer 


PR: 


high-speed paper tape reader 


PP: 


high-speed paper tape punch 


DTn: 


DECtape unit n; n=0 is the default condition 


MTn: 


Magtape unit n; n-fi is the default condition 


BI: 


BATCH input data set 
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The filename may be any 1- to 6~character alphanumeric name. The filename extension 
may be any 1- to 3-character alphanumeric sequence. The filename extensions assumed 
or supplied on default for the FORTRAN Compiler (and for LINK) are shown in 
Table 7-10. 

Table 7-10 
Filename Extensions 



File 


Assumed Extension 
on Input File 


Default Extension 
on Output File 


FORTRAN Compiler 

object file 
listing file 
source file 

LINK 

load module 
object module 


.FTN 
• OBJ 


.OBJ 
.LST 

• LDA 



When a source file of the given filename (having no extension) cannot be found on 
the device specified using the defeiult extension, the system attempts to find 
the filename with a null extension. 

9.1.2 Compilation and Linking Procedures 

To activate the FORTRAN Compiler under the DOS Monitor, the command 

$RUN FORTRN 

is given. The FORTRAN Compiler then prints FORTRAN and the compiler version 
number, and a # to indicate that it is ready to accept a command string. 

A sample FORTRAN command sequence is shown below. 



$RUN FORTRN 
FORTRAN Vxxx 



#OBJECT , LIST<FILE1 
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This command string directs the compiler to take the source file FILEl.FTN from the 
system device and output the files LIST.LST and OBJECT OBJ to the system device. The 
user then calls the DOS/BATCH Linker to link the object module (s) with the FORTRAN 
Library as follows. 

( 1) Link (2) Link and Execute 

$RUN LINK : $RU LINK 

LINK Vxxx LINK Vxxx 

#LOAD<OBJECT,FTNLIB/L/E # LOAD<OB JECT , FTNLIB/L/GO 

Input to LINK is one or more compiled object modules and any required library files. 
Output from LINK is a single load module and, optionally, a load map and symbol table 
file. 

In the linkage shown in example (1) above, the linker output is a load module, which 
can then be loaded and run by means of the RUN command. The /L switch indicates that 
FTNLIB is a library file and the /E switch indicates the end of input to the linker. 
In example (2) , the /GO switch causes the linker to output the load module to the 
system device, to load a copy of that load module, and to begin execution of the 
program, as though a 

$RUN LOAD 

command had been given. When the /GO switch is used, the /E switch is omitted from 
the linker command string. 

If any of the FORTRAN source files contains more than a single subroutine or a single 
main program, the resulting object module produced by the FORTRAN Compiler is a con- 
catenated object module file. When such a file is linked by LINK, the linker /CC 
switch must be used following the object module file specification. Failure to use 
the /CC switch with a concatenated object module results in an error message from 
LINK and a defective load module. 

The following sequence might be used to compile, link, and execute a FORTRAN system 
consisting of 

a. the FORTRAN main program MAIN.FTN, 

b. the FORTRAN subroutine SUBR1.FTN, and 

c. several FORTRAN subroutines in the file UTILTY.FTN. 

SRUN FORTRN 

FORTRAN V06,09 

#MA1N,UP|<MA.|N 

#SUBRt,lPKSUBRl 

WlLTYiLPKUTtlTY 

*tC 

' KI 

5RUN LINK . 

LINK Vxxx 

m.LP'KMAWfSUBRliUTlLTY/CCiFTNUB/LVE 

#t0 

'.KI 

SRUN Pi 
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Either of the compiler output files can be eliminated by omitting its file specifica- 
tion from the command string. For example, 

$RUN FORTRN 
FORTRAN V06,09 
#FlLEi<FILE:l 

produces FILE1.0BJ on the system device but no listing file, while 
#_,LP:<FILE1 

produces a listing on the line printer, but no object module output. 

9.1.3 Compile-Load-and-Go Operation 

Rather than give separate command strings to both the FORTRAN Compiler and LINK, a 
single command string to the FORTRAN Compiler can be used to cause compilation, 
linking and execution. This one-step compile- load-and-go sequence is performed as 
follows s 

$RUN FORTRN 
FORTRAN Vxxx 
£FILE1 . LST<FILE1/G0 

in which case the source program FILE1.FTN on the system device is compiled and 
linked to the FORTRAN Library (FTNLIB) . A listing file (FILEl.LST) , an object 
file (FILE1.0BJ) and a load module file (FILEl.LDA) are created on the system 
device. The load module file is then loaded and execution begun. 

9.1.4 FORTRAN Library Usage 

By means of the DOS Librarian, the \iser can construct his own libraries of machine 
language and FORTRAN routines, which he then searches at link time. The user is 
constrained to search all of his own libraries before searching the. FORTRAN System 
Library . 

Users should not add subroutines or functions to the DEC- supplied FORTRAN Library. 
Many routines in the OTS Library are order-dependent and the insertion of new 
routines could result in undefined global references when linking FORTRAN programs. 
Similarly, the deletion or rearrangement of routines in FTNLIB is likely to cause 
similar problems unless the user is familiar with the ordering dependencies of the 
FORTRAN Library. Therefore, users should create separate library files, using the 
Librarian. Thus, if MATLIB is a user library containing matrix manipulation 
routines, and the user writes a program (PROG) that uses routines from MATLIB, a 
command string to the Linker might be: 

#PROG , LP : <PROG ,MATLIB/L ,FTNLIB/L/E 
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9.2 FORTRAN OUTPUT LISTING FORMAT 

When a listing file is requested, each page contains a header with the following 
information : 



Compile version identification 
current time 
current date 
listing page number 

Using the /LI:n switch option, the user can vary the type of listing file obtained; 
from a mere listing of error diagnostics to a complete listing of source, assembly 
language and symbol table. All listing files include block summaries describing 
the current program unit. 

The block summary or descriptor block is a collection of information about a 
FORTRAN main program or subroutine, printed following the END statement on the source 
listing (before the assembly listing, if one is requested) . Figures 7-7 and 7-8 
show sample FORTRAN programs and their descriptor blocks. A descriptor block 
provides the following data: 

a. A listing of the routines called from the program. 

b. A listing of the switches specified in the command string. 

c. A listing of the size of the program in both decimal words and 
octal bytes. Sizes are also provided for each common block declared 
within the program. Since common blocks may be shared among several 
programs, their sizes are of individual interest. 

The information is formatted into three columns under two headings: 
BLOCK and LENGTH. The name of the program or common block is 
listed under BLOCK, followed by the size in decimal words, followed 
by the size in octal bytes (in parentheses ). These last two entries 
are under the LENGTH heading. 

The line containing the program name is followed by an asterisk. 
Where no program name is supplied, the name MAIN, is shown. The 
name .$$$$. denotes blank common; the name DATA, indicates a block 
data subprogram. 

d. Following the Assembly listing, where one is requested, a final data 
item is supplied, entitled: 

♦♦COMPILER CORE** 

It has three subheadings: PHASE, USED, and FREE. 

These terms head columns that contain, respectively, the phase of the 
compiler (three phases) . The number of decimal words of core storage 
used in that phase, and the number of decimal words of core storage 
free .in that phase. 
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FORTRAN 


V06.12 


17:16:59 23-MAY-73 


PAGE 


1 


0001 




SUBROUTINE SUB2X(I) 






0002 




DOUBLE PRECISION DOS 






0003 




COMMON Tl, Z, DOS /L2/ J2 (10) , V2 , W2 






0004 




Tl-1 






0005 




Z— 1.0 






0006 




DOS=-l .0D0 






0007 




W2 = SQRT (EXP (J2 (I) -1.0)) 






0008 




V2= -1.0 






0009 




RETURN 






0010 




END 








ROUTINES CALLED: 








SQRT , 

OPTIONS 

BLOCK 
SUB2X 
.$$$$. 
L2 


EXP 

== /ON,/CK,/OP:l 

LENGTH 
94 ( (000274)* 
7 (000016) 
24 (000050) 








**COMPILER CORE** 








PHASE USED FREE 








DECLARATIVES 00366 18276 








EXECUTABLES 00543 18099 








ASSEMBLY 00965 20594 







Figure 7-7 
Block Summary Example #1 



FORTRAN V06.12 

0001 
0002 
0003 
0004 
0005 

0006 

BLOCK 
DATA. 
LI 

L2 



17:17:18 



-23 -MAY -7 3 



PAGE 



BLOCK DATA 

COMMON /LI/ XBD, YBD, CBD, DBD /L2/ UBD, VBD, WBD 
COMPLEX CBD 
DOUBLE PRECISION DBD 

DATA XBD, YBD, CBD, DBD, UBD, VBD, WBD / 2*1.0, 
1 (1.0, 1.0), 1.0D+0,3*1.0 / 
END 





12 

6 



LENGTH 

(000000) 
(000030) 
(000014) 



**COMPILER CORE** 

PHASE USED FREE 
DECLARATIVES 00366 19276 
EXECUTABLES 00550 18092 
ASSEMBLY 00845 20714 



Figure 7-8 
Block Summary Example #2 



7-108 



These sizes describe the symbol table space used in each phase and 
provide the user with a ratio of used to unused space in the 
compiler. 

If the FORTRAN program was compiled without subroutine calls or switches, the 
information described in (a) or (b) , respectively, is not provided. The information 
in (c) and (d) is provided as part of the listing for every FORTRAN program or 
subroutine. 

9.3 COMPILE-TIME MEMORY REQUIREMENTS 

The DOS /BATCH FORTRAN Compiler uses all available memory space for the program 
symbol table. A simple variable entry in the symbol table is eight words long. 
A constant entry is eight words plus the size of the constant. An array entry 
is ten words plus one word for each dimension. For example, a complex constant 
entry is twelve words long and a 3-dimensional array entry is thirteen words long. 
The user can reduce the amount of memory used by the compiler by his choice of 
input/output files in command strings and by reducing the number of continuation 
lines permitted. Memory space saved in this way is available as additional symbol 
table space. « ■ . 

The /CO switch, which controls the number of allowable continuation lines on a 

given FORTRAN statement, takes space from the symbol table area. Thus, /CO: 99 

takes the maximum amount of space from the symbol table while /CO:0 makes 

additional space available to the symbol table. (The default value is /CO: 5). 

# 

The devices used for I/O affect the compiler core requirements. The minimum core 
required for I/O operations would be used with the command sequence 

$RUN FORTRN 
FORTRAN Vxxx 
#FILE<FILE 

where source and object files are on the system device. Memory requirements are 
increased, and symbol table space correspondingly reduced, when the user requests 
listings or uses devices other than the system device. For example, 

#FILE,KB:<DT:FILE 

requires considerably more memory space than the previous command string. 
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PART 7 

CHAPTER 10 

FORTRAN OPERATING ENVIRONMENT 

10.1 FORTRAN OBJECT TIME SYSTEM 

The FORTRAN Object Time System (OTS) is composed of the following: 

a. Mathematics routines, including all standard FORTRAN functions plus other 
arithmetic routines needed to do arithmetic operations (e.g. , floating 
point) . 

b. Miscellaneous utility routines (PDUMP, SETERR and SETFIL) . 

c. I/O Routines, which handle the various types of FORTRAN I/O. 

d. Error handling routines, which process arithmetic errors, I/O errors and 
system errors. 

e. Miscellaneous (Polish) routines required by the compiled code ($SBS, 
$DOEND, $POP, etc.). 

The library is designed as a large number of small pieces so that unnecessary rou- 
tines can be omitted at link time. Thus, if the user performs only sequential 
formatted I/O, none of the random access routines are linked to his program. 

10.2 THREADED CODE 

Most FORTRAN compiled statements generate calls to library subprograms. These calls 
are based on the technique of evaluating expressions in Polish notation, which 
breaks down expressions into a large number of simple operations performed in a 
linear sequence. These operations use the PDP stack for evaluating all expressions 
and subexpressions. 

The implementation of Polish notation makes several assumptions. 

a. The first operation done in a Polish sequence is always a "push". 

b. It is not necessary to place breakpoints (as in ODT) in the middle of an 
arithmetic statement. 

c. Speed does not suffer by assignment of a register (R4) for special purposes. 

As an example, the FORTRAN statements , 

A=l.p 
B=l.j0 
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generate the following code. 



;A=1.0 

$Ppppi 

.GLOBL $P0P3 

$P0P3,A 

;B=l.p 

$PJ00J01 

$P0P3,B 



Routines such as those above, referenced directly or indirectly, are inserted into 
the user program and are found at the end of the assembly listing; other routines 
are linked to the user program from the FORTRAN Library. .The routines inserted 
into the user program as a result of the above code follow. 



$Pjdj0 ) 01 : MOV 
BR 

$R0€00 : 04020J& 
00000 

$Fj6f6fil: MOV 
MOV 
JMP 



#$R0j0je)j0+4,R0 
$Fl6/6j6l 



-(Rtf),-(SP) 
-(Rtf),-(SP) 
@(R4) + 



;GET ADDRESS OF VALUE 

? FLOATING-POINT CONSTANT 1.0 

;PUSH 2-WORD VALUE ONTO 

; STACK 

;GO TO NEXT ROUTINE 



The routine $P0P3 is in the FORTRAN Library. $P0P3 pops a real value off the stack 
into the memory location whose address follows the call to $P0P3 in the threaded 
code. 

$P0P3,A 

pops the value on top of the stack into the two memory words reserved for A. 
Similarly, 



$P0P3,B 



places the 2-word value found on top of the stack into B. As another example, the 
statement 



A =B+C 



generates code similar to the following 
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$P0)601 

$ptf02 

.GLOBL $ADR 
$ADR 
$POP3,A 
.+2 



; PUSHES THE VARIABLE B ONTO THE STACK. 
;EACH OPERATION CONSISTS OF THE ADDRESS OF A 
; ROUTINE! TO BE EXECUTED. A PUSH PLACES A 
; VALUE ON THE STACK; A POP REMOVES A VALUE. 

;THIS LINE CAUSES AN EXIT FROM POLISH MODE; 
; NORMAL EXECUTION IS RESUMED. 



$Pjeftfjefl and $J?tf02 push the values of B and C onto the stack; $ADR is a FORTRAN 
Library routine to add two floating-joint numbers; $POP3 saves the result in the 
variable A. The subroutine calls for the above sequence follow. 



$pje^i: 


MOV 


#B+4,R0 




BR 


$Fje&l?fl 


$VjtM02 : 


MOV 


#C+4 f R0f 


MOV 


-(R0),-(SP) 




MOV 


-(RjZf),-(SP) 




JMP 


@(R4) + 


$POP3 : 


MOV 


(R4) + ,R3 




MOV 


<SP)+, (R3) + 




MOV 


(SP)+ f (R3)+ 




JMP 


@(R4) + 



GET THE ADDRESS OF B. 
JUMP TO PUSH. 
GET THE ADDRESS OF C. 
PUSH TWO 
WORDS ON STACK. 
JUMP TO NEXT. ROUTINE. 
GET ADDRESS OF 
VARIABLE DESTINATION. 
POP A VALUE 
FOR THE VARIABLE. 
;GO TO NEXT ROUTINE. 



Note that the instruction 



JMP 



@(R4) + 



jumps to the next routine in the list, as well as incrementing R4 over that item in 

the threaded code. All internal functions are called in this manner and exit using 

such a jump instruction and must clear any stack space used (except for the return 
value left on top of the stack) . 

All routines explicitly called by the user (i.e., subroutines and external func- 
tions) are called using the PDP-11 subroutine calling convention (see Chapter 16) . 



1#.2.1 Entry to Polish Mode 



Entry to Polish mode is made via a call to the $POLSH routine, as follows. 



JSR R4 , $POLSH 



This invokes the following routine. 
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$POLSH: TST (SP)+ ; DELETE OLD VALUE OP 

;R4 PUSHED ON ENTRY 
JMP . @(R4) + ; EXECUTE NEXT THREADED CODE CALL 

The word following the JSR to $POLSH is the first word of Polish code to be executed. 

ltf.2.2 Exit from Polish Mode 

An explicit exit from Polish mode is made by means of a word containing the address 
of the following word. For example: 



;IN POLISH MODE 
.WORD .+2 ; LEAVE POLISH MODE 

/•CONTROL PASSES TO THIS LOCATION 



ljtf.2.3 Polish Mode Subroutine Calls 

The PDP-11 FORTRAN calling sequence convention is described in Chapter 7-15. The 
PDP-11 FORTRAN Compiler (V0fc) implements this convention by means of Polish Mode 
operators described in this section. 

A subroutine or function call is performed using the $CALL service routine followed 
by one argument (the routine to be called) . The argument list follows immediately 
in memory. The routine called is entered in direct execution mode. (This is true 
of FORTRAN compiled subprograms as well as assembly language programs.) 

Thus, it is necessary to re-establish Polish mode after each subroutine or function 
call using one of the library routines $RPOL0 or $RPOLN. These routines are similar 
to $POLSH but $RPOLN additionally adjusts the processor stack register to delete 
any temporary arguments created for the call. For example the statement 

X - FNC(A+B,C) 

will result in the following code. 
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$Fjef0#l : 



$P0001 

$Pj0J0#2 

$ADR 

$SVSP,$Fj0^0'l 

$CALL,FNC 

BR .+6 

i 

+C 



JSR 
+4 



%4 , $RPOLN 



;PUSH VALUE OF A ON STACK 

;PUSH VALUE OF B ON STACK 

;ADD REAL - RESULT LEFT ON STACK 

;COPY ADDRESS OF A+B INTO ARG LIST 

;TWO ARGUMENTS 

; FILLED IN AT EXECUTION TIME 

ASSEMBLED IN ADDRESS OF C 

/ARRIVE HERE IN NON-POLISH MODE 

; RETURN TO POLISH MODE AND 

• DELETE TEMPORARY VALUE A+B 



l£f.3 FORTRAN RUN-TIME MEMORY ORGANIZATION 



The memory map during execution of a FORTRAN program is shown in Figure 7-9, 



xx 



75JZ#: 



R6 



tf: L 



OTS Routines 

and 

FORTRAN Program 



Stack 



I/O Buffers and 
Device Drivers 



Resident 
Monitor 



Interrupt 
Vectors 



) Monitor Free Core Area 



Figure 7-9 
Memory Map Organization 
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As files are opened, space is taken from the Monitor free-core area to accommodate 
additional I/O buffers and device drivers. The FORTRAN program and OTS routines 
use the processor R6 stack for storage of temporary results. Overflow of the 
stack into the areas dynamically allocated by the Monitor may cause DOS fatal 
errors or system crashes. 

I/O buffer space and device drivers are released and the memory reclaimed when the 
logical unit is closed by an END FILE statement. END FILE may be used on both in- 
put and output files. 

DOS I/O uses space from the free core area as follows: 

1. One device driver routine for each device that has been initialized 
(.INIT) but not released (.RLSE). (The END FILE statement performs 
a .CLOSE and .RLSE) . The approximate DOS driver sizes are given in 
Table 7-11. 

2. One I/O buffer for each open dataset (FORTRAN logical unit in use). 
The buffer size depends upon the device. 

3. One Device Data Block (DDB) for each open dataset (16 words). 

4. One File Information Block (FIB) for each open file (16 words). 

5. One bit map for each unit on which a file is open for output (64 
words) . 

Table 7-11 
Approximate DOS Device Driver Sizes 





Device Driver 


I/O Buffer 


Device Specifier 


Size (Words) 


Size (Words) 


DK 


130 


256 


DF 


70 


64 


DC 


70 


64 


DP 


260 


512 


DT 


19jrf 


256 


MT 


46J0 


256 


LP 


160 


48 


CR (ASCII only) 


210 


48 


KB 


600 


32 


PP 


8jtf 


32 


PR 


70 


32 
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UZf.4 DEBUGGING FORTRAN PROGRAMS 

The use of traditional FORTRAN debugging techniques such as PDUMP is recommended for 
development of FORTRAN applications (see Chapter 7-14) . 

System error reporting and traceback information is of significant value in debug- 
ging. Figure 7-10 is an example of the traceback feature in the FORTRAN system. 

At run time, diagnostics are printed by OTS with a trace of the flow of control 
within the user-written code. Following each error code are printed the error 
headings NAME and SEQ, below which are the names of the routines through which the 
call is being traced and the sequence number of the specific line in which the 
error occurred (or that from which that subprogram was called) . The first name and 
sequence number at the top of the list is the error location; subsequent names and 
numbers refer to the path through which the program reached that point. 

The example in Figure 7-20 shows the same error message being printed twice. Though 
the two errors were generated by the same line in the subprogram, they are traced to 
different lines in the main program., 

It requires considerable experience to successfully use ODT with a FORTRAN program. 
ODT was not intended to handle the problems of FORTRAN debugging; it is especially 
hard to use when trying to debug threaded, code. 

The threaded code generated by the compiler is a sequence of addresses, rather than 

machine instructions. ODT breakpoints can be placed only on machine instructions. 

The user is therefore constrained to use breakpoints only in places where the code 

leaves Polish Mode; e.g., a subroutine or function call. Alternatively, breakpoints 
can be placed in the Polish routines themselves. 

There is one significant inconvenience associated with putting breakpoints in Polish 
routines. Polish routines are usually called from several places in a program and 
when the breakpoint in a Polish routine is encountered, the user must look at R4 to 
find the address from which the routine was called. 

When executing a FORTRAN program, that stores data outside the memory area allocated 
for an array, the Monitor messages F342 and F344 may occur and the system may halt. 
The compiler /CK switch (Appendix K) causes all subscript references to be checked 
for upper and lower bounds. Use of /CK is recommended while debugging FORTRAN 
programs . 
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0j6tii 

tf#02 
J6003 
00JZ4 
0JS0S 



BLOCK 
MAIN. 



J=2 

CALL MUL (J) 

J=4 

CALL MUL(J) 

END 

LENGTH 
51 (000146)* 



Main Program 



J 



0001 
0002 

0004 
i 0005 



BLOCK 
MUL 



SUBROUTINE MUL (J) 

DO 10 1=1,5 

J=J*J 

RETURN 

END 

LENGTH 
46 (j%0i34)* 



"\ 



Subprogram 



$RU EXAMPLE 

FORH003014 PRODUCT OUTSIDE OF RANGE ON INTEGER MULT. 
NAME SEQ 
MUL 003 
MAIN. 0002 , 



1 



FORT003j814 PRODUCT OUTSIDE OF RANGE ON INTEGER MULT. 
NAME SEQ ^ 

MUL #5903 
MAIN. 0004 



Run Time 
Diagnostics 



Figure 7-10 
Example of Run Time Diagnostics 
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10. 5 FORTRAN OTS ERROR PROCESSING 

The FORTRAN OTS error diagnostics are divided into nine classes (Classes &-Q) . The 
messages associated with each diagnostic are given in Appendix K. 

The FORTRAN OTS maintains an error class table, which governs the actions taken for 
each error class. The user can modify the values in the error class table (and 
hence the error processing) by means of the SETERR subroutine call (Chapter 7-14) . 

Each error class entry contains values for the following variables: 

COUNT a count of the number of errors that have occurred for this 
class. 

FLAG set to 1 on any occurrence of an error in the class. May be 

tested and reset, by means of CALL TSTERR (Chapter 7-14) . 

MAX an integer specifying the action to be taken on an error in 

the class. 

The actions taken by the OTS when an error condition occurs are determined by the 
value of MAX as follows: 

MAX>jef Set FLAG, log error message on output device, increment COUNT, 
and call EXIT if COUNT>MAX. 

MAX=£f Set FLAG, log error message on output device, and continue 
program execution. 

MAX=-1 Set FLAG, do not log message, ignore count, and continue program 
execution. 

MAX=-2 Do not log message, EXIT to Monitor, and print 1351 message. 

MAX=-3 Immediate exit to Monitor with a message of the form: 

F,03# J0CQ0NN 

where CC and NN are the error class and number (in octal) of 
the condition causing the exit (see Appendix K) . 

The error classes and default MAX values are detailed in Table 7-12. 
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Table 7-12 
FORTRAN OTS Error Classes 





Default 




Class 


MAX Value 


Meaning 





-3 


Fatal errors; it is recommended that the 
occurrence count (-3) not be changed. 


1 


ff 


Physical I/O errors. 


2 


1 


Errors in FORMAT statements. 


3 


3 


Arithmetic overflow or division by 0. 


4 


4 


Incorrect arguments to library functions or 
subroutines . 


5 


-1 


Arithmetic underflow errors. 


6 





Conversion errors. 


7 


1 


Subscripting errors. 


8 


1 


Errors involving incorrect linkages to OTS 
routines; the maximum occurrence count (1) 
cannot be changed , although no diagnostic 
is given. 



For example, arithmetic underflows (class 5) are ordinarily ignored by the FORTRAN 
OTS; arithmetic overflows (class 3) are logged and the default error limit is three 
occurrences. 
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?ART 7 

CHAPTER 11 

FORTRAN MONITOR I/O CONSIDERATIONS 

Input/output operations for a FORTRAN-compiled program are performed by 1:he FORTRAN 
Object Time System. The FORTRAN language statements related to I/O are described 
in Chapters 7 and 8. It is the purpose of this Chapter to describe the .processing 
of I/O requests in greater detail, with emphasis on the FORTRAN OTS and DOS /BATCH 
Monitor mechanisms and representations. File formats and interactions between DOS 
and the FORTRAN OTS are also described. 

The FORTRAN OTS accomplishes all physical I/O through the DOS Monitor. FORTRAN I/O 
is therefore device-independent to the extent that such independence is supported 
by the DOS /BATCH Monitor. For example FORTRAN formatted WRITE statements can be 
directed to any ASCII output device, but direct-access I/O can only be done on 
devices having fixed-length addressable blocks (disks and DECtape, but not industry- 
standard magnetic tape) . 

11.1 INPUT/OUTPUT OVERVIEW 

As described in Chapter 8, all FORTRAN I/O operations are done in terms of logical 
records. Each READ or WRITE steitement transmits one or more logical records. The 
overall flow of control in processing a single FORTRAN READ or WRITE statement 
appears in Figure 7-11. 



I/O call from FORTRAN program 



I 



Initialize I/O Operation 



I 



Transmit Items Specified in 
FORTRAN I/O List 



Finish I/O Operation 



Return Control to 
FORTRAN Program 



Figure 7-11 
FORTRAN I/O Flow 
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The operations accomplished follow. 

1. Initialize I/O Operation: 

a. Verify legal unit number, obtain $DEVTB entry. 

b. Test to see if the file is open; if it is not, open the file. 

c. If the file is open, check for compatibility of device and file for 
the type of I/O requested. 

If a contiguous file is to be allocated, this is done as a part of the 
file-opening process in step (b) above. 

2. Item Transmission: 

a. If the operation is a READ, obtain a record if needed; check for I/O 
errors. 

b. Transmit items according to the I/O list. If the record length is 
exceeded, issue a diagnostic message. 

3. Finish I/O Operation: 

a. Complete the current record with padding and record separators as 
required . 

b. If the operation is a WRITE, transmit the record and wait for 
completion. 

c. Evaluate physical I/O error returns. 

In the above descriptions, a file is not opened until a READ or WRITE statement for 
that unit is executed. Thus, errors such as nonexistent files or inability to 
locate a contiguous file are not detected until an I/O operation is attempted. 

All FORTRAN I/O is single-buffered. A DOS .WAIT macro is issued after each Monitor 
I/O call. All transfers of control on END= and ERR= precisely identify the FORTRAN 
I/O statement causing the END= or ERR= transfer. 

The principal data structure used by the FORTRAN OTS for controlling I/O is the 
device table, $DEVTB, described in Sections 7-11.4 and 7-11.5. The unit number 
supplied in a READ or WRITE statement is used to select a device table entry. The 
device* table entry contains information for each logical unit, including: 

Filename. 

Physical device specification and unit number. 

User identification code (UIC) . 

Protection code. 

Open/close status. 
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FORTRAN logical unit defaults are specified by the values assembled into ?DEVTB in 
the version supplied with the FORTRAM OTS. CALL SETFIL and CALL ASSIGN (Chapter 
7-14) can be used to modify $DEVTB entries under program control. The DOS ASSIGN 
command can be used to override some of the $DEVTB attributes. 

The physical devices supported by FORTRAN and DOS, and the types of FORTRAN opera- 
tions permitted on each device, are specified in Table 7-13. 



Table 7-13 
DOS/BATCH FORTRAN Standard Peripheral Devices 





Device 


FORTRAN 


FORTRAN 




Name 


Specification 


I/O Type 


I/O Operations 








Direct- 


READ WRITE ENDFILE 


BACKSPACE 






Formatted Unformatted Access 




REWIND 


DC: ^ 












DF: 










Disk* 


\ or Sy: 


X X- X 


X X X 


X 




DK: 










DECtape 


DP: 
DTn: 




X XX 


XXX 


X 


Magtape 


MTn: 


X X 


X X X 


X 


Line Printer 


LP: 


X 


X X 




Card Reader 


CR: 


X X 


X X 




Terminal 


KB: 


X 


XX X 




High-speed 
paper tape 
reader 


PR: 


X X 


X X 




High-speed 
paper tape 
punch 


PP: 


X X 


X X 




Low-speed 
paper tape 
reader/punch 


PT: 


X 


XXX 




* The universal mnemonic SY can be used to specify the 




system disk regardless of whether that unit is an 
RC11, RF11, RK11, RP#2, or RP,#3 disk. 




X Indicates that a specific I/O mode or operation is 
supported on the indicated device. 
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11.2 FILE STRUCTURES AND FORMATS 

FORTRAN input/output facilities are provided by three packages of OTS routines: 
formatted I/O, unformatted I/O, and direct-access I/O. General characteristics of 
DOS file structures and the three I/O packages are discussed in following sections. 

The DOS Monitor provides a variety of I/O modes, and two principal file structure 
organizations : 

a. Linked files for sequential access. 

b. Contiguous files for direct or sequential access. 

Linked files consist of a series of blocks, which need not be physically contiguous 
on the device. Each block contains a pointer to the next block of the file. When 
a linked file is initially opened for output, the first block is allocated. As out- 
put requests occur, additional blocks are allocated and linked to the file. The 
length of the file is therefore unknown at the time the file is opened. 

Formatted and unformatted output ordinarily uses linked files. A contiguous file 
may be used for output following a file allocation operation through SETFIL. 

Contiguous files are specifically intended for direct-access I/O. A fixed-size 
contiguous file must be initially allocated on disk or DECtape by means of DEFINE 
FILE or CALL SETFIL. Direct-access I/O statements may then read or write individual 
records without having to read all preceding records. 

Table 7-14 summarizes th*e file structures and DOS I/O modes used by the FORTRAN I/O 
packages. 

Table 7-14 
File Structures and I/O Modes 



FORTRAN 
I/O Type 


DOS 
File Structure 


DOS 
I/O Mode 


Formatted 

Unformatted 

Direct-access 


Linked or Contiguous 
Linked or Contiguous 
Contiguous 


Formatted ASCII Normal 
Formatted Binary Normal 
Unformatted Binary 
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11.2.1 FORTRAN Formatted I/O 

The formatted input/output routines read or write variable-length formatted ASCII 
records. A record consists of a maximum of 133 ASCII characters transmitted under 
control of a format specification, followed by record separator characters (typi- 
cally carriage return/line feed). For example, the format specification 133A1 
transmits a maximum length record. 

On output to a nonprinting device (e.g., disk file), the record separator consists 
of the carriage return/line feed sequence. The maximum length output line thus 
requires 135 bytes on a disk file. The file consists simply of a sequence of out- 
put lines; there are no additional headers, checksums, or byte counts. 

On output to a printing device (KB: or LP:) the record separator appended to each 
record consists of a carriage return and a vertical tab. The first character of 
each record is deleted from the record and is interpreted as a carriage control 
character; it may cause output of 0,, 1 or 2 line feed characters or a, form feed 
character. 

On input, records longer than 135 bytes (including separator) are truncated. For 
shorter records, the last character (line feed, form feed or vertical tab) is 
deleted. The next-to-last character is deleted if it is a carriage return. No 
other special interpretations are provided; that is, form feeds, tabs and other 
special ASCII characters are passed to the program as single characters. In par- 
ticular, tab characters are not converted to blanks., 

Note that interpretation of the carriage control character is provided by the 
FORTRAN OTS, not by the Monitor or utility programs. If output containing FORTRAN 
carriage control characters is directed to a disk or tape file for later printing , 
that printing must be done under control of a FORTRAN program in order to obtain 
the desired carriage control. 

11.2.2 FORTRAN Unformatted I/O 

The unformatted input/output routines read or write variable length binary records 
using the DOS/BATCH formatted binary mode. Each record contains checksums, for 
parity checking, and a byte count. 

A FORTRAN unformatted WRITE statement outputs a single logical record. A, logical 
record consists of one or more segments, each segment containing a maximum of 61 
words of data from the user program and a segment control word supplied by the 
FORTRAN OTS. The segment control word is inserted as the first data word of the 
segment and may take one of the values shown in Table 7-15. 
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Value 




1 
2 
3 



Table 7-15 
Segment Control Word 



Meaning 



Neither first nor last segment 
First segment of logical record 
Last segment of logical record 
First and last segment of logical record 
(Single- segment record) 



Each segment of the logical record is output as a DOS formatted binary line. The 
format of each segment is shown in Figure 7-12. 



word 0: 
word Is 
word 2: 



^ 

y 



byte count 



segment control 
. word 



user data 
maximum of 
122 bytes 



Checksum 
byte 



supplied by DOS 



supplied by FORTRAN 
OTS 



S 



supplied by FORTRAN 
program I/O list 



supplied by DOS 



DOS-11 (Y08-J3f2 and later) supplies a single null charac- 
ter, if necessary, to end a logical record on a word 
boundary; otherwise, no further characters are appended 
to the file by the system. 



Figure 7-12 
Logical Record Segment Format 

Each segment has four words of overhead; the user must allow for these overhead 
words when allocating a contiguous file for unformatted I/O through CALL SETFIL. 
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11.2.3 FORTRAN Direct-Access I/O 

The direct-access I/O routines read or write fixed-length binary records in a con- 
tiguous file. The DOS unformatted binary I/O mode is used. The maximum record 
length is 32,767 bytes. 

The logical record structure for a direct-access file is determined entirely by the 
DEFINE FILE statement. The records themselves contain no checksums, byte counts, 
or record separators. Different programs can process the same file using different 
logical record structures. In Figure 7-13, Program 1 operates on records consisting 
of 1$ x 10 integer matrices, while Program 2 operates on records consisting of a 
single row of a matrix. 



Program 1 

INTEGER MATRIX (10, 10) 



DEFINE FILE 1 {20, 100, U, INDX) 



WRITE OUT A 10 x ISf MATRIX 

WRITE (l'NREC) ( (MATRIX (I, J) , J=l,l#) ,1=1,10) 



END 

Program 2 

INTEGER ROW (1120 



DEFINE FILE l(20Qf, 10, U, INDEX) 



READ A 10 ELEMENT ROW OF MATRIX 
READ (l'NREC) ROW 



END 



Figure 7-13 
Program Example Using Logical Records 
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If a WRITE statement transmits fewer words than will fit in the record, the remain- 
der of the record is filled with zeros. 

The direct-access record structure is independent of the physical block size of the 
I/O device, but more efficient operation results if the record size is an exact 
divisor or multiple of the physical block size (RF and RC disks = 64 words, RPjZfe or 
RK disk and DECtape = *256 words, and RP03 disk - 512 words). 

11.3 FORTRAN I/O ERROR HANDLING 

Three classes of run-time error conditions are specifically associated with FORTRAN 
I/O operations . 

Class 1 Nonexistent file, physical I/O error, etc. 

Class 2 FORMAT syntax errors 

Class 6 Conversion errors during format conversions. 

A complete list of FORTRAN run-time errors is contained in Appendix K. 

As described in Chapter 8, READ and WRITE statements can specify transfer of control 
on an end-of-file condition (END=) and on certain physical I/O conditions (ERR=) . 
The SETFIL and ASSIGN subroutines (Appendix F) can specify an error variable, IERR, 
which is set to indicate the specific type of I/O error. 

Errors in classes 2 and 6 are processed as described in Section ljrf.5 and Appendix 
K. Class 1 I/O error conditions are processed as follows: 

Assume the error condition is identified by the diagnostic FORTjeflBlptlN. 

1. Set the IERR variable, if specified in the previous SETFIL to NN. 

2. If END=n is specified, and the error condition is an end-of-file 
(FORT#0Oj9J0r4) , transfer control to statement n. 

3. If ERR=n is specified, transfer control to statement n. 

4. Otherwise, process the error according to Section 1-10.5. 

11.4 FORTRAN DEVICE- ASSIGNMENTS 

The default FORTRAN device assignments are shown in Table 7-16. These assignments 
can be temporarily altered through use of a call within the user program to the 
SETFIL or ASSIGN routines or by use of the ASSIGN Monitor command. Section 7-11.5 
describes how to permanently alter the FORTRAN device assignments. 
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Table 7-16 
FORTRAN Logical Device Assignments 



Logical Unit 




Default 


Number 


Device 


Pile Name 


1 


system disk, SY: 


FOR033OL.DAT 


2 


system disk, SY: 


FOR0p2.DAT 


3 


system disk, SY: 


FOR0J8T3 . DAT 


4 


high-speed paper 
tape reader, PR: 


FORB^.DAT 


5 


line printer, LP: 


FORGJ0S.DAT 


6 


terminal keyboard, KB: 


F0RPJy6.DAT 


7 


system disk, SY: 


FOR007.DAT 


8 


*EiATCH system input 
device , BI : 


FOR008.DAT 



Device assignments are determined in one of three ways. These are described below 
in order of increasing priority. 

a. The device table default values described in Appendix D. 

b. A call within a FORTRAN program to the SETFIL (or ASSIGN) subroutine 

(see Chapter 7-14) overrides the default device assignments. 

c. The Monitor ASSIGN command overrides both default assignments and any 
assignments made within the program by a call to the SETFIL routine. 

The SETFIL routine causes entries to be made in the FORTRAN device table, $DEVTB, 
which is described in greater detail in Section 7-11.5. 

An ASSIGN Monitor command for use with FORTRAN is specified in the form: 

$ASSIGN dev:file.ext[uic] , n 

where the $ character is printed by the system; dev:file.ext[uic] represents a file 
specification, and n represents a FORTRAN unit number. For example, 

$ASSIGN DT:,2 

forces a search on DECtape unit for the filename (and extension) associated with 
logical unit 2 in the device table (or in a SETFIL call) . 
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$ASSIGN DT3:FILE.TMP,3 

causes the file FILE.TMP on DECtape unit 3 to be associated with logical unit 3. 

Alterations in the device table made by a call to the SETFIL subroutine are only 
valid during the program in which they appear. Device assignments made with DOS 
ASSIGN command vary in scope depending upon when the command is given. 

a. When no program is in core, assignments made through an ASSIGN command 
are valid through several runs of the same or different programs; they 
can be removed by an ASSIGN command with no arguments. For example: 

$KI 

$AS DT:,1 
$RU X 
$RU Y 

During programs X and Y, DECtape unit ft is treated as logical unit 1. 

b. When a program is in core, assignments made through an ASSIGN command 
are valid only through the running of that program. For example: 

$GE X 
$AS DT:,1 



$KI 

During the run of program X, DECtape unit is treated as logical unit 1. 
The KILL command returns the device assignments to their original default 
conditions. 

c. When a program is in core and executing, an assignment made by an ASSIGN 
command is valid only through the running of that program. For example: 

$RU X 

PLEASE ASSIGN INPUT ON UNIT 4 

h00S 100 ^| (PAUSE message) 

$AS SY:FILE1.XX,4 

$C0 

During the run of program X, the FORTRAN program wrote a message requesting 
a device assignment and then executed a PAUSE statement. The user then 
indicated that the desired file was to be read from the system disk. 
Following the execution of program X, the device* table reverts to its de- 
fault assignments. . . . 

11.5 THE FORTRAN DEVICE TABLE, $DEVTB 

Logical device assignment is governed by the device table, which contains entries 
for eight devices. Additional entries can be inserted, default entries changed, or 
entries deleted by reassembling the device table file. To allow for such changes, 
the source file for the device table (DVB. MAC) is supplied with each FORTRAN system. 
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When the source has been altered and assembled/ the new device table object module 
can replace the supplied module in the FORTRAN Library or be linked as an object 
module (to individual programs) preceding the library file. The device table con- 
sists of a 12-word header followed by a l"6-word entry for each device. 



Table 7-17 
Device Table Entry 



$DEVTB : 



Header 
Word 

Word 1 
Word 2 
Word 3 
Words 4-11 



Word 


1 


Word 


2 


Word 


3 


Words 4-5 


Word 


6 


Word 


7 


Word 


8 


Word 


9 


Words 1^-14 


Word 


15 


Word 


16 



Address of device table entry for logical unit -3 
(the message logging unit) 

Address of entry for error routine message file 

Number of entries in device vector table 

Device number of message logging file 

Addresses of device table entries for each of the 
devices 1 through 8 



Link pointer (from Link Block, after .INIT) 

Physical device name stored in Radix-Sj^ format 

Bits 15-8 = unit number (Default #) ; Bits 1-JX 
- DOS Open Code, taken from File Block-2 

Filename (stored as .RAD507FOR/,/NNN/; NNN = 
Entry number) 

Filename extension (stored as .RAD50/DAT/) 

Switches* and protection code (Default = 233) 

Status/Mode (from Line Buffer Header) 

Count of I/O operations for this device 

Unused for formatted and unformatted I/O 

User Identification Code (UIC) (default=#, 
indicating the current UIC) 

Address of error variable (IERR from CALL 
SETFIL) 
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Table 7-17 (Cont.) 



For Random I/O, Words 


8- 


14 


are as follows: 


Word 8 






Function word 


Word 9 






Block number 


Word 1$ 






Buffer address 


Word 11 






Buffer length * 


Word 12 






Associated variable address (from DEFINE FILE) 


Word 13 






Maximum number of records (from DEFINE FILE) 


Word 14 






Record length (from DEFINE FILE) 


*Switches are as follows 


: 




Bit 


Setting Meaning 


0-1 
2 




1 
2 
3 

1 


Closed 

Open formatted 

Open unformatted 

Open random 

DEFINE FILE done on this device 



Table 7-17 describes the device table header and individual entry formats. The 
supplied device table is listed in Chapter 7-17. 

In order to add a device to the table , 

1. Alter the word at $DEVTB to reflect the number of entries in the device 
table (the number of devices available to FORTRAN) . The supplied value 
at $DEVTB is 8. 

2. Append the address of the device entry to the 8-word header sequence of 
device table entry addresses (known as the device table entry vector) . 

3. Append the new 16-word entry describing the additional device. 
To delete a device from the table, 

1. Alter the word at $DEVTB to reflect the entries in the device table. 

2. Set the address of the device entry in the device table entry vector to 
2ero. 

3 Delete the entry for that device from the device table. 
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The default physical device for a table entry can be changed by modifying the second 
word of the entry. This allows the user to alter the associations between devices 
and logical unit numbers shown in Table 7-17. 

The default filename and extension for a table entry can be changed by modifying 
words 4 through 6 of the table entry. 

11.6 I/O EXAMPLE 

This section presents an example using the FORTRAN I/O system. The program FORDGN.- 
FTN on the following pages is used in the FORTRAN system to read an ASCII source ' 
file and produce a contiguous file of specially formatted records containing the 
English language diagnostic messages for the FORTRAN Compiler and OTS. A listing 
of the message file is also printed on the line printer in the example following. 

The program reads an input file from logical unit 4. The first record of the input 
file specifies the number of error messages in the output file and the name of the 
output file. CALL SETFIL sets the name of the output file and the DEFINE FILE 
statement specifies the record structure. 

The source listing of FORDGN.FTN follows. 

FORTRAN V06. 13 14:52:12- 01-APR-74 PAGE 1 

C FORTRAN SYSTEM DIAGNOSTIC MESSAGE FILE BUILDER 

C THIS PROGRAM CAN BE USED TO BUILD FILES FOR 

C EITHER DOS-il OR RSX-11D AND IT CAN BE 

C RUN UNDER EITHER SYSTEM BY MAKING THE APPROPRIATE 

C CHOICE OF "ASSIGN" OR "SETFIL" BELOW 
C 

C CREATES AND THEN PRINTS A FILE OF MESSAGES 

C FOR ACCESS BY THE FORTRAN COMPILER OR OTS 
C INPUT: 

C FILE - AS SPECIFIED BY KEYBOARD TYPE- IN 

C LUN - 4 

C 1ST RECORD— 13, 40A1 I3=#0F 64 CHARACTER MESSAGES 
C TO BE ALLOCATED. 40A1=FILE SPECIFICATION 

C OTHER RECORDS- 12, 64A1 I2=P0SITI0N OF CURRENT 
C MESSAGE IN THIS FILE. 64A1=CURRENT MESSAGE 

• C LAST RECORD— 13, 64A1 I3=NEGATIVE INTEGER. 64A1=IGN0RED 
C OUTPUT : 

C FILE- FILE SPECIFICATION READ INTO FILSPC 

C LUN - 1 

C 64 CHARACTER FIXED LENGTH RECORDS 
C 

C LUN - 5 <: NORMALLY LP:> 

C PRINTED OUTPUT 
C 
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0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 



0009 
0010 



0011 
0012 



0013 
0014 



0015 
0016 



1000 

10001 

C 

C 

c 
c 



C 
C 

1001 

C 
C 
C 
C 

c 

c 
c 

10 

c 



INTEGER COUNT, INDEX 

BYTE T0DAV<9), FILSPCC40), A<64), NULL < 64) 

DATA NULL / 64 * / 

COUNT=0 

WRITE <6, 1000) 

FORMATS *SPEC I FV INPUT FILE>' > 

READ<6, 10001) 18, FILSPC 

FORMAT < G, 40AD 

SET UP THE INPUT FILE NAME ASSIGNMENT 

CALL ASSIGN<4, FILSPC, IO) 

CALL SETFIL<4, FILSPC, IERR, 'SV',0) 

ENDFILE 6 

READ FIRST RECORD AND ALLOCATE CONTIGUOUS FILE 
READ < 4, 1001) IBLOK, FILSPC 
FORMAT <I 3, 40A1) 

SET UP THE OUTPUT FILE NAME ASSIGNMENT 

CALL ASSIGN<i, FILSPC, 40) 

CALL SETFIL<1, FILSPC, I ERR, ' SV , 0) 

DEFINE FILE KIBLOK, 32, U, INDEX) 

INITIALIZE ALL RECORDS OF OUTPUT FILE TO ZERO*: 
DO 10 1=1, IBLOK 
WRITE < I'D NULL 



FORTRAN V06. i: 



14 ; 53 : 12 



01-APR-74 



PAGE 



0017 
0018 
0019 
0020 



C READ INPUT AND WRITE EACH RECORD OF OUTPUT FILE 
40 READ <4, 1002, END=50)IO, II, A 
1002 FORMAT<G, 13, 64AD 

IF <I1 . LT. 0) GO TO 50 

COUNT=CGUNT+l 



0021 

0022 
0023 



10=10-3 

IF <I0 . LE. O) GO TO 40 
IF <I0 . LT. 64)A<I0+1)=© 



0024 
0025 
0026 



0027 
0028 
0029 

0030 
0031 

0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 

0040 



C 
C 
C 
50 



1003 

C 



51 

1004 
1005 



11=11+1 
WRITE Cl'IDA 
GOTO 40 

READ THE CONTIGUOUS FILE AND 

PRINT A LISTING OF THE COMPLETED FILE 

ENDFILE 4 

CALL DATE<TODAV) 

WRITE <5, 1003) FILSPC, TODAV 

FORMAT <1H1, 28X, 40A1/1HO, ' MSG MESSAGE' , 17X, 9A1/' NUM' /> 



DO 51 I 
READ<1' 
12=11-1 
WRITE < 
FORMAT < 
WRITE < 
FORMAT < 
ENDFILE 
ENDFILE 
END 



1=1, IBLOK 
IDA 

5.1004) 12, A 
X, 13, X, 64 AD 

5.1005) IBLOK, COUNT 

/X, 13, ' MESSAGES ALLOCATED. ',16, 
1 
5 



MESSAGES INPUT. ' ) 
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ROUTINES CfiLLED: 
SETFIL, DATE 

OPTIONS =/OP:l 

BLOCK LENGTH 

MR IN. 584 < OO2220 >* 

**CGMPILER CORE** 

PHASE USED FREE 
DECLARATIVES ©0622 06070 
EXECUTABLES 00:363 05829 
ASSEMBLV 01355 09977 



The input file, ABC. SRC: 

018ABC. TST 

002 FORT0O0002 SUBROUTINE DIRECTLV < I ND I RECTLV) REFERENCES ITSELF 

003 FORT000003 ILLEGAL FLOATING POINT INSTRUCION 

004 FORT000084 SVSTEM ERROR NO DIAGNOSTIC MESSAGE ASSIGNED 

000 FORT00O000 INVALID CALL TO ERROR 

001 FORT000001 NO SPACE TO DO I/O 

005 FORT000005 SVSTEM ERROR NO DIAGNOSTIC ASSIGNED 

006 FORT001000 SVSTEM ERROR NO DIAGNOSTIC ASSIGNED 

007 FORTO01001 DEVICE PAR I TV ERROR 

008 FORT801002 CHECKSUM/PAR ![ TV ERROR-END OF DATA ERROR (RANDOM) 
010 FORT001004 END OF FILE OR END OF MEDIUM 

012 FORT0010O6 DEFINE FILE NOT DONE (RANDOM ACCESS) 

013 FORT001087 DEFINE FILE DONE <NOT RANDOM ACCESS) 

015 FORT001009 FILE DOES NOT EXIST / OR IS ALREADV OPEN 

016 FORT001010 UNABLE TO OPEN FILE 



Before execution of the program, the disk contains the following. 

#/DI 

DIRECTORY DK0 : C 50,51 3 

01-APR-74 

FORDGN. FTN 6 01-APR-74 <233> 

ABC . BAK 2 01-APR-74 <233> 

FORDGN. OBJ 8 01-APR-74 <233> 

FORDGN. LDA 20 01-APR-74 <233> 

ABC . SRC 2 ©l-APR-74 <233> 

TOTL BLKS: 38 
TOTL FILES: 5 
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The program is run and the output file ABC.TST is created. 

* ASSIGN SV:RBC. SRC, 4 
*RUN RZiRDGN. LDfl 

SPECIF V INPUT FILE>flBC. SRC 

*RUM PIP . 

PIP V10-92 
. #LP:</DI 

# 

DIRECTORY DK0| C 50,5* ] 
03»APR*74 



FORDGN.FTN 
ABC ,BAK 
PORDGN.QBJ 
FORDGN.tDA 
ABC ,SRC 
ABC ,TST 


6 
2 

8 
20 
2 
3C 


0l-»APR-»74 
01*APR*74 
0t»APR-74 
01*APR»74 
01*APR*74 
01*APR«74 


<233> 
<233> 
<233> 
<233> 
<233> 
<233> 


TOTL BLKSl . 
TOIL FILES* 


41 

6 







The listing of the file produced on the line printer: 



ABC. TST 



MSG MESSAGE 
NUM 



©l-APR-74 



1 
2 
3: 
4 
5 
6 

8 
3 
18 
11 
12 
12 
14 
15 
16 
17 



FORT000000 

FORT000001 
FORT000002 
FORT00000S' 
FORT000004 
FORT000005 
FORT001000 
FORT001001 
FORT001002 



INVALID CALL TO ERROR 

NO SPACE TO DO I/O 

SUBROUTINE DIRECTLV < INDIRECTLV)REFERENCES ITSELF 

ILLEGAL FLOATING POINT INSTRUCTION 

SVSTEM ERROR NO DIAGNOSTIC MESSAGE ASSIGNED 

SYSTEM ERROR NO DIAGNOSTIC ASSIGNED 

SVSTEM ERROR NO DIAGNOSTIC ASSIGNED 

DEVICE PARITY ERROR 

CHECKSUM/PARITY ERROR-END OF DATA ERROR < RANDOM) 



FORT0O1004 END OF FILE OR END OF MEDIUM 

FORTO01006 DEFINE FILE NOT DONE (RANDOM ACCESS) 
FORTO0100? DEFINE FILE DONE <NOT RANDOM ACCESS) 

FORT0O1009 FILE DOES NOT EXIST / OR IS ALREADY OPEN 
FORT001O10 UNABLE TO OPEN FILE 



18 MESSAGES ALLOCATED. 



% 



14 MESSAGES INPUT. 
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PART 7 

CHAPTER 12 

FORTRAN WORD FORMATS 



12.1 INTEGER FORMAT 



n+2 



Sign 



0=+ 

1— 


binary number 




15 14 










15 



In a 2-word format, an integer is assigned two storage words, although only the 
high-order word (i.e., the word having the lower address) is significant. Use of 
the /ON switch (see Appendix J) causes integers to be assigned a single storage 
word. Negative integers are stored in a two's complement representation. Integer 
data must lie in the range -32768 to +32767, For example: 



+22 = 000026 
-7 = 177771. 



8 



12.2 FLOATING-POINT FORMATS 

All floating-point data, in both 2-word and 4-word formats, is stored as a 1-bit 
sign, an 8-bit exponent (characteristic) , and a fractional mantissa. The exponent 
is stored in excess 128 notation. That is, a binary zero represents an exponent 
of -128; the binary equivalent of 255 (377.) represents an exponent of +127. The 
mantissa is stored as a binary fraction; that is, the most significant bit repre- 
sents 0.5, the next bit represents .0.25, the third, '0.125, and so on. Because all 
numbers are presumed to be normalized, the most significant bit is redundant and 

is therefore not stored (this is called hidden bit normalization) . If the exponent 

—128 v 

is zero (corresponding to 2 ) , the most significant bit is assumed to be 0; for 

all other exponent values, it is assumed to be 1. A value of zero is represented 

by two or four words completely filled with zeros. As an example of floating-point 

storage, the value +1.0 is represented in the 2-word format by 
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or in the 4-word format by the following. 



04^2 #0 



This is equivalent to an exponent of 129 (corresponding to 2 ) multiplied by a 
mantissa of zero (assumed to have a value of #. 5). The value -5.0 in the 2-word 
format is 



140640 



or in the 4-word format, the following. 



140640 
000000 
000000 
000000 



12.2.1 Real Format (2-Word Floating-Point) 



word n: 



word n+2: 



Sign 




1=- 


Binary excess High-order 
128 exponent mantissa 


15 14 


7 6 Qr ■ 


Low-order mantissa 



15 







Real floating-point format gives an effective precision of 24 bits, or approximately 
7 digits of accuracy. The magnitude range lies between approximately #.14 x llf 3Q 
and 1.7 x 10? 8 . 

12.2.2 Double Precision Format (4-Word Floating-Point) 



word n: 



Sign 



$=+ Binary excess High-order 
1=- 128 exponent mantissa 



15 14 



7 6 
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word n+2: 



word n+4; 



word n+6: 



3'jow-order mantissa 


15 




? 




Lower-Border mantissa 




15 




a 


Lowest-order mantissa 



15 



y 



The effective precision is 56 bits, or approximately 16 decimal digits, or accuracy. 

—38 -*38 

The magnitude range lies between 0. 14 x 10 and 1.7 x 10 . 



12.3 COMPLEX FORMAT 



word ns 



word n+2: 



Sign 



J2f~+ Binary excess High-order 
1~- 128 exponent mantissa 



15 14 



7 6 



T 



Low- order mantissa 



15 



T 



Real 
Part 



word n+4i 



word n+6: 



Sign 



0~+ Binary excess High-order 
1«- 128 exponent mantissa 



15 14 



7 6 



J* 



Low- order mantissa 



15 



Imaginary 
Part 



T J 
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12.4 BYTE FORMAT 



Unspecified 



Data Item 



15 



8 7 



.? 



The range of numbers from +127 to -128 can be represented in BYTE (LOGICAL *1) 
format. BYTE format array elements are stored in adjacent bytes. 

12.5 HOLLERITH FORMAT 



word 0\ 



word 2: 





char 2 




char 1 




15 




8 


7 







char 4 




char 3 




15 




8 


7 


? . 




blank=4j3f 
8 




char n (n<255) 





15 



8 7 



jy 



Hollerith constants are stored internally as one character per byte, filling up to 
word boundaries. An odd number of characters causes a blank to be appended to the 
constant. A Hollerith constant included in an arithmetic expression is treated as 
a 1-word integer quantity where only the first two characters of the constant are 
significant. 



12.6 LOGICAL FORMAT 



true: 



false: 



.15 



1 


7 


7 


7 


7 


7 


15 
















V 





t 


f 






t 



Logical format data items are treated as 1-word integer values for use with arith- 
metic and logical operators. Any nonzero value is considered to have a logical 
value of .TRUE, when tested by a logical IF statement. 

12.7 RADIX-SJZI' FORMAT ' 

See Appendix B. 
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PART 7 
CHAPTER 13 
FORTRAN LIBRARY FUNCTIONS 

This Chapter contains a brief outline of the OTS library of FORTRAN functions that 
involve approximations. Floating-point means single-precision, 2-word, floating- 
point format with a 24-bit fraction and an 8-bit binary exponent. Double-precision 
means 4-word, floating-point format with a 56-bit fraction and an 8-bit binary ex- 
ponent. The values of the coefficients used in the various approximations may be 
found at the cited parts of the following references. 

a. Computer Approximations , by J.F. Hart et al, John Wiley & Sons, 1968. 

*>• Approximation for Digital Computers , by C. Hastings et al, Princeton 
University Press, 1955. 

c. PDP-11 Paper Tape Software Programming Handbook , DEC-11-GGPC-D, 
Digital Equipment Corporation. 

In the descriptions of the various functions, the relative error values given are 

for the approximating polynomials in the cited intervals. These error bounds 

assume exact arithmetic. There are two additional sources of errors in the function 
calculations that are not considered in the error bounds. 

a. Rounding and truncation errors can occur in reducing the given . argument 
to the range in which the polynomial or rational fraction approximations 
are valid. 

b. Rounding errors can occur as a result of using finite precision floating- 
point arithmetic in the polynomial or rational fraction computations. 

All OTS FORTRAN functions are called using one of the following standard sequences 
(see Chapter 7-15) . 

MOV #LIST,RE 

JSR R5,NAME or JSR PC, NAME 

LIST: BR RTN LIST: BR RTN 

.WORD ARGl,...,ARGn .WORD ARG1, . . . ,ARGn 

RTN: RTN: 

The result is returned in R0-R1 for floating-point and RJ2C-R3 for a double-precision 
function. 
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Some FORTRAN functions call other single argument functions in the course of their 
computation. In order for them to be re-entrant, these calls are made via the 
routine $FCALL. 

MOV ARGUMENT ADDRESS, R5 
MOV #FUNCTION NAME, R4 
JSR PC,$FCALL 

$FCALL calls the FORTRAN function whose address is in R4 with the argument whose 
address is in R5. Control is returned to the instruction following the JSR with 
the function result in R0-R1 for floating-point or R0-R3 for double-precision. 

Here are the FORTRAN Library approximation functions. 

1. ALOG(X), Floating-point Natural Logarithm 

If X<0 call error 

Let X=y*2 a where l/2<y<l 

Let Q=(y*/2- l)/(y*/2~+ 1) 

Then ln(X) = a * In (2) + ln{y) 

ALOG(X) = a* In (2) = In (/~2) + Qlj? c.Q 21 

—9 9 

where the c. are drawn from Hart #2662. The relative error is <10 ' . 

2. ALOG10'(X) , Floating point Common Logarithm 
Computed as log .(e)*ALOG(X) . 

3. ATAN(X), Floating-point Arctangent 

If X<0, ATAN(X) = -ATAN (-X) 
If |x|>l, ATAN(|x|) = TT/2 - ATAN (l/|x|) 
• If |x|>tan TT /12, ATAN(X) = TT /6 + ATAN ( (X /T-1)/(X + /Tj") 
For |x|<tan tt/12, ATAN(X) = XE c.*X 21 

—9 5 

where the c. are drawn from Hart #4941. The relative error is <10 . 

4. ATAN2 (X,Y) ,Two Argument Floating-point Arctangent 

If Y = tf, or X/Y >2 , ATAN2(X,Y) = TT/2 (sign X). 

If Y > ,0V and X/Y<2 2 , ATAN2(X,Y) = ATAN (X/Y). 

25 
If Y < 0, and X/Y<2 , ATAN2(X,Y) = TT*signX + ATAN (X/Y). 

5. DATAN(X), Double-precision Arctangent 

The analysis is the same as in that for ATAN(X) except that the polynomial 

approximent is of degree 8. The coefficients are drawn from Hart #4945. 

-16 8 
The -relative error is <\ft ' . 
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6. DATAN2(X,Y), Two Argument Double-precision Arctangent 

The rules for DATAN2 are the same as those for ATAN2 except that the 
DATAN is used in all computations. 

7 . DLOG (X) , Double-precision Natural Logarithm 

The analysis for DLOG is the same as that for ALOG except that the poly- 

2 
nomial in Q is of degree 6. The c. are drawn from Hart #2665. The rela- 

. „w-16.5 1 

tive error xs <lp 

« 

8 * DLOGlj^(X) , Double-precision Common Logarithm 

Computed as log ^ (e)*DLOG(X). 

9 * DSQRT(X), Double-precision Square Root 

If X<0 call, error 
Let X = A*2 B where 1/2 < A <1 
Let Y , = 2 B/2 * (1/2 + A/2) if B is even 
or 
Y^= 2 (B+1)/2 * (1/4 + A/2) if B is odd 

a transformation requiring only two instructions. Starting with Y_,, four 
Newton-Raphson iterations are performed. y 

Vl =:1/2(y n + x/ V- 

-17 
The relative error is <10 

10. DSIN(X), Double-precision Sine 

Let y - Integer (4 * fraction (X/2ir) ) 
Let V => Fraction (4 *fraction (X/27T) ) 
Then DSIN(X) ■ P(Vir/2) if ytf 

- P((l-V)iT/2) if y=l 

- P(-W/2) if y*=2 

- P((V-l)7T/2) if y=3 

8 2i 
where sin (W/2) *?m/2) - vEc.V for -1 < V < 1 

* " ' -18 6 

The c. are drawn from Hart #3345. The relative error is <10 * • 

11. DCOS (X) , Double-precision Cosine 
Computed as DSXN(X + TT/2). 

12. DEXP (X) , Double-precision Exponential 

If X > 87 call overflow 
If |x| < 2~ 6 ^ j DEXP(X) - 1 
If X < -88 „ 7 j DEXP (X) ■ $ 
Let y - Integer (X*log_ (e) ) 



7-142 



Let V = 16 * Fraction (X • log (e) ) 

Let w = 1/16 * Fraction (V) 

DEXP = 2 y * 2 W * 2 Inte ^ er (V) / 16 where ff < w < 1/16. 

1/16 
Powers of 2 are obtained from a table. 

2 2 
o w z P(v )+wQ(w ) „ ■ J « *• ^ 
2 ~ ^ * — r— where P and Q are first 

P(w )-wQ(w ) 

2 
degree polynomials in w . 

The coefficients of P and Q are drawn from Hart #1121. 

-16 4 
The relative error is <lfi " . 

13. EXP (X) , Floating-point Exponential 

If X > 87; call overflow 
If |x| < 2 ~ 28 ; EXP(X) =1. 

If X < -88.7, (EXP(X) = fl. 

Let y = Integer (X * log (e)) 

Let V = Fraction (X *log_(e)) 

Let w = 1/2 ln(2) * V where |x| < ln(2)/2 

Then EXP(X) = 2 Y * (e W ) 2 

where e ~ 1 + 2 * w 



c - w - c 

^ -2 



c 3 + w 
The c are drawn from DEC-11-GGPC-D , page 7-23. The relative error is 

14. RAN (II ,12) and CALL RANDU (11,12 ,F) , Random Number Generator 

RAN and RANDU use the same multiplicative congruential algorithm for 
generating uniformly distributed pseudo-random numbers in the interval 

If 11=0, l2=jZC set generator base 



otherwise 



X .. - (2 16 +3) X mod 2 32 
n+1 n 



Store generator base X , in 11,12. 

n+l 

Result is X .. scaled to a real value Y ,,, &<Y <1. 
n+l n+l ^ n+l 
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15 . SIN(X) , Floating-point Sine 

The analysis is the same as that for DSIN(X). The polynomial approximant 

used is of degree 4 and the coefficients are drawn from Hastings, sheet 

-8 
16. The relative error is <_ % *\tf . 

16 • COS (X) , Floating-point Cosin e 

Computed as SIN(X + TT/2) . 

17* SQRT (X) , Floating-point Square Root 

The analysis is the same as that for DSQRT(X) except that only three 
iterations are performed. The relative error is <^ 10 . 

18. TANH (X) , Floating point Hyperbolic Tangent 

If |x| > 16, TANH(X) * 1 * sign(X) 

otherwise 

let y - EXP (2 * X) 

TANH(X) = (y-l)/(y+l). 
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PART 7 

CHAPTER 14 

SYSTEM SUBPROGRAMS 



14.1 LIBRARY ARITHMETIC FUNCTIONS 

Refer to Table 7-3 for the FORTRAN Library arithmetic functions. 

14.2 SYSTEM SUBROUTINES AND FUNCTIONS 

The FORTRAN library contains subroutines that may be called in the same manner as 
user-written subroutines, similar to the functions intrinsic to the FORTRAN system. 
These subroutines appear in Table 7-18. 

Table 7-18 
FORTRAN Library Subroutines 



System 
Subroutine 


Function 


PDUMP 


performs core dumps of specified sections of 
core. 


SETPDU 


changes logical unit to which PDUMP output is 
written. 


SETFIL 


overrides default value for a logical unit at 
run time. 


ASSIGN 


allows specification at run-time of filename 
or device and filename to be associated with 
a FORTRAN logical unit number. 


SETERR 


allows the user to reset the maximum occurrence 
count for any class of error except £f. 


TSTERR 


returns an indication of whether an error of 
the specified class has occurred. 


RANDU,RAN 


returns a single random number with uniform 
distribution between & and 1. 


EXIT 


terminates the execution of a program and re- 
turns control to the Monitor (also invoked 
automatically by the error processing routines 
and STOP statement) . 


DATE 


returns a 9-byte string containing the ASCII 
representation of the current date. 


TIME 

- - 


returns an 8-byte string containing the ASCII 
representation of the current time in hours, 
minutes and seconds. 
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Table 7-18 (Cont) 
FORTRAN Library Subroutines 



System 
Subroutine 


Function 


SECNDS 

IRAD50 
RAD50 

R50ASC 

SSWTCH 


provides system time of day or elapsed time as a 
floating-point value in seconds. 

performs conversion of Hollerith strings to 
Radix- 50 representation. 

converts Radix-50 strings to. Hollerith strings. 

tests specified bits of the console switch 
register. 



14.2.1 PDUMP 

The PDUMP subroutine causes specified portions of core to be dumped. Control re- 
turns to the calling program following execution of the dump. The call to PDUMP 
is performed as follows: 

CALL PDUMP (L, ,U, ,F„ , .. . . ,L ,U ,F ) 
111 n n n 

where L is the lower limit (a variable name) and U is the upper limit of the area 
to be dumped. The order in which the two variables appear may be reversed; regard- 
less of the order of specification, memory is always dumped from the lower limit 
to the upper limit, inclusive. For example, 

DIMENSION N(10) 



CALL PDUMP (N ( 1 ) , N ( 1/tf ) , 3 ) 

causes the array N to be dumped on the logical unit 5 (generally LP:) as a string 
of 2-word integers from N(l) to H(ltf). The following line performs exactly the 
same dump. 

CALL PDUMP (N(1^),N(1),3) 
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The value of F is an integer constant or variable indicating the dump format, as 
follows. 



code 


Description 


Format 


i> 


1-word octal 


808 


1 


1-word decimal 


818 


2 


2 -word octal 


808 


3 


2-word decimal 


818 


4 


real 


4G15.7 


5 


double precision 


3D22.14 



Codes jl and 1 should be used if the /ON (1-word integer) switch is set; codes 2 
and 3 should be used if the /ON switch is not set. Codes ]/f and 1 print every word 
while codes 2 and 3 print every other word, since the alternate words are unused. 
No default dump mode is assumed; a code must be specified or a diagnostic will be 
generated. Where several sections of core are to be dumped, each set of three 
values (L. ,U,,F.) specifies a separate core area and dump format. 

No spacing is performed between dumps of individual sections, nor are addresses 
printed on the PDUMP output. 

As an example, the output shown below is produced using the call: 

CALL PDUMP (B(3) ,B(3) ,jeJ,R(l) ,R(7) /y 0,F (1) ,F(4) ,5,S (1) ,S (24) ,f6) 

where the program contains: 



DOUBLE PRECISION F(4) 
INTEGER B(3),R(7),S(24) 



and the /ON (pne-word integer) switch is used at compilation time. 



1404 



14 



#. 49303812265972D-31 
22060 34770 17756 
4467 12672 
327jSfe 1 20024 



12746 

# j6 20124 20130 17776 

0. 0.197215245087#LD-30 0.98607624310639D-31 



50002 


30 


21172 


406 


20130 


32240 








50561 


55740 


4567 


1412 


400 


4467 


12640 
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14.2.2 SETPDU 

The SETPDU subroutine allows PDUMP output to be directed to any legal logical unit 
number. The call to SETPDU is performed as follows: 

CALL SETPDU (IUNIT) 

where IUNIT is a constant, variable or expression of type INTEGER designating the 
logical unit desired. 

When used in conjunction with the ASSIGN subroutine (see Section 7-15.2.4), PDUMP 
output may be directed to any file or device on the system. 

14 . 2 . 3 SETFIL 

The SETFIL subroutine overrides the default value for a logical device assignment 
at run time. The SETFIL call must occur before the file in question has been 
opened for I/O operations (by a READ or WRITE statement) . The new logical device 
assignment remains in effect until the end of the current program run or until the 
file is closed by END FILE and a new device assignment made. The call to SETFIL 
appears as shown, below: 

CALL SETFIL (N,FILE, IERR,DEV,NU,UIC,PC,CS,LREC,NREC) 

SETFIL accepts a. variable number of arguments; nonessential arguments may be 
omitted. The complete list of arguments to this call appears in Table 7-19. 



Table 7-19 
FORTRAN SETFIL Arguments 



Argument 



Significance 



N 
FILE 



Logical unit number, expressed as an integer value. 

Character string containing the file name and extension 
specification, of the form: 

file.ext 

where "file" may be up to six alphanumeric* characters, of which 
the first is alphabetic. Optionally, ".ext" may be specified 
to indicate an extension. The entire file specification is an 
ASCII string that can be expressed as a Hollerith string or 
enclosed in single quotes. For example: 
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Table 7-19 
FORTRAN SETFIL Arguments (Cont.) 



Argument 


Significance 






•MAT.TMP' 
8HOBJ2.DAT 




See also the discussion of string argument storage in Section 
F.2.13. 


IERR 


Integer variable into which the error returns from this routine 
or from I/O operations on logical unit "N" are placed. IERR is 
set to -1 if CS is neither 1 nor 2; it is set to 2 if CS=1 and 
neither LREC nor NREC is specified. 


DEV 


A 2- or 3-character device specification, expressed as a 
Hollerith string or enclosed in single quotes. For example: 




•DT' 
2HLP 


NU 


Unit number of the device specified, if appropriate. This is 
an integer constant or variable (e.g., 1 if the device is DEC- 
tape unit 1, or if a line printer or other single-unit device). 


UIC 


User identification code (UIC) , expressed as an integer value. 
A UIC specification of $ indicates the current UIC (the usual 
specification) . A simple way of preparing the UIC entry is 
shown below (assume the desired UIC is 120,15j2f) 




INTEGER UIC 
LOGICAL*! UICB (2) 
EQUIVALENCE (UICB (1) , UIC) 
UICB(l) = "150 
UICB (2) - "12$ 




CALL SETFIL (2, 'DATA.TMP* ,IERR, 'DT' ,2, UIC) 




which specifies that the file DT2:DATA.TMP [120,15$ is to be 
accessed on logical device 2. 


PC 


Protection code , specif ied as an integer value, which is the 
decimal equivalent of the desired octal protection code (pro- 
tection codes are defined with octal values in Part 3 of this 
handbook) . The default value for a FORTRAN protection code is 
233(8). 


CS 


Integer value used to indicate one of the following: 




CS = 2, allocate a contiguous file for random I/O. The 
DEFINE FILE statement sets the record length and 
number of records in the file. 




CS = 1, allocate a contiguous file for formatted or un- 
formatted I/O. The arguments LREC and NREC are 
used to set the record length and number of re- 
cords in the file. 


LREC 


Logical record length in words, expressed as an integer value; 
required if CS=1, ignored otherwise. 


NREC 


Number of records to allocate, expressed as an integer value;, 
required if CS=1, otherwise ignored. 
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The discussion of file formats in Section 7-11.2 should be consulted as a guide in 
choosing values for LREC and NREC. 

When SETFIL is used in an overlay system, a special restriction exists for the 
variable IERR. In the statement 

CALL SETFIL (7 , 'FILE. DAT' ,IERR, 'DK' ,0) 

the address of the variable IERR is saved and is used by the I/O system for report- 
ing certain kinds of I/O error conditions. This variable name is saved by SETFIL 
in the device table, $DEVTB. Any subsequent I/O errors (during READ or WRITE op- 
erations) attempt to set IERR. Therefore, IERR must be in the local storage of the 
resident section or in a resident common block. If it is not, the FORTRAN I/O 
system may destroy an unpredictable location in memory when attempting to store 
error data in IERR. Further discussion of the use of SETFIL and of the FORTRAN I/O 
system is presented in Chapter 7-ll Two examples of the use of SETFIL are given 
below. These examples assume the default logical unit assignments. 

Example 1: Make logical unit 5 the card reader and logical unit 6 the 
line printer. 



CALL SETFIL (5, ' A' , IERR, 'CR' ) 
CALL SETFIL (6, ' A' , IERR, ' LP 1 ) 



Note that the parameters to SETFIL are positional. The CR: and LP: are not file 
structured, so the file name argument is superfluous; a file name argument must be 
included, however. 

Example 2: Allocate a file for direct-access I/O on disk unit 1 under 
the current UIC. 



CALL SETFIL (1, 'FILEl.XYZ' , IERR, 'SY' ,0,0,100,2) 
DEFINE FILE 1 (20, 200, U, INDX1) 



The file FILEl.XYZ is created on the system disk with protection code 100 (read, 
write, execute access; the file cannot be deleted until its protection code is 
changed) . 



7-150 



14.2.4 ASSIGN 



The ASSIGN subroutine assigns, at run time, the device and/or file name to be 
associated with a logical unit number. The ASSIGN call must occur before the 
logical unit is opened for I/O operations (by READ or WRITE) . The assignment 
remains in effect until the end of the program or until the file is closed by 
END FILE and a new CALL ASSIGN is performed. The call to ASSIGN is performed as 
follows. 

CALL ASSIGN (N , NAME , ICNT , IERR) 

The arguments to this routine follow. 

Argument Significance 

N Logical unit number, expresses as an integer value 

NAME Character string containing any device and file name acceptable 

to the Monitor system. 

The maximum acceptable length for a name string is 3$ characters. 

If the device is not specified, the device remains unchanged 
from that specified by default or by a previous ASSIGN or SETFIL 
call. If the device is specified, the filename must also be 
specified in the case of file-structured devices. (If the de- 
vice is not file-structured, a file name, if present, is 
ignored . ) 

ICNT Actual length of the name string in characters (bytes) ; if the 

value is 0, then the name string itself will be assumed to be 
terminated by a zero byte. 

IERR Integer variable into which error code values from this routine, 

or I/O operations on logical unit "N" are to be placed. 

The error variable name is an optional argument and may be omitted. If specified, 
it should not be allocated in an overlay segment of the program. 

14.2.5 SETERR 

The system maintains an error count for each error class. The SETERR subroutine 
allows the user to reset the maximum occurrence count for any class of run.-time 
error (except class 0) to 8. (Following the maximum number of errors in a given 
class, control exits to the Monitor.) The call to SETERR is performed as follows. 



This subroutine is the preferred alternative to the SETFIL subroutine when random 
access file allocation is not required. 
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CALL SETERR (CLASS, MAX) 



The argument CLASS is an integer indicating the error class affected. If CLASS is 
not a valid number, no action is performed. The CLASS argument is shown in rela- 
tion to its associated error messages in Appendix K. The error classes appear here 
in brief. 



Class 

i 

1 
2 
3 
4 

5 
6 
7 



Default 
MAX Value 

-3 



-1 

i 

7 



Meaning 

Fatal errors; it is recommended that the maximum 
occurrence count (-3) not be changed 

Physical I/O errors 

Errors in FORMAT statements 

Arithmetic overflow or division by Jo 

Incorrect arguments to library functions or subrou- 
tines 

Arithmetic underflow "errors 

Conversion errors 

Subscripting errors 



The argument MAX is an integer with the following significance: 



MAX Value 

-1 
-2 
-3 



Meaning 

log message, increment error count, call EXIT if 
error count is > MAX., 

Log messages and ignore error count for the specified 
error class. 

Do not log messages, ignore error count for the 
specified error class. 

Do not log messages, exit to Monitor after printing 
an error message and closing files. 

Immediate exit to Monitor with a message. 



14.2.6 TSTERR 

The TSTERR subroutine allows the user program to monitor the types of errors detec- 
ted during program execution. The call is of the form 

CALL TSTERR (I, J) 
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where I is the error class number (between 4 and 8) ; a value is returned to the 
variable J as follows: 



J»l if an error of class I has occurred. 
J«2 if an error of class I has not occurred, 



The sequence 



CALL TSTERR(3,J) 
GO TO (10,20), J 
2$ CONTINUE 



transfers control to statement 10 if a class 3 error (arithmetic overflow) has 
occurred. See Section 7-14.2.4 for a discussion of error classes and messages. 

The TSTERR routine also resets to $ the error flag for that error class (but not 
the error count used by SETERR) . For example : 



CALL TSTERR (I, J) 
CALL TSTERR (I, J) 



The second call is guaranteed to return J=2. The TSTERR subroutine is independent 
of the SETERR subroutine; neither directly influences the other except that SETERR 
can cause execution to terminate. 

14.2.7 RANDU, RAN 

The random number generator may be called as a subroutine, RANDU, or as an intrin- 
sic function, RAN. The subroutine call appears as 

CALL RANDU (II, 12, X) 

where II and 12 are previously-defined integer variables and X is the real variable 
name in which is returned a random number between <j> and 1. II and 12 should be 
initially set to <jt. They are updated to a new generator base during each call. 

7-153 



Resetting II and 12 to repeats the random number sequence. The values of II and 
12 have a special form; only $ or values supplied by the random number generator 
should be stored in these variables. 

Use of the RAN function is similar to the use of the random number subroutine: 

RAN(I1,I2) 

is the function reference to the random number generator. 

14.2.8 EXIT 

A call to the EXIT subroutine, in the form 

CALL EXIT 
is equivalent to the END statement and causes program termination. 

14.2.9 DATE 

The DATE subroutine can be used in a FORTRAN program to obtain the current date as 
set within the system. The DATE subroutine is called as follows: 

CALL DATE (array) ■ . 

where array is a predefined array able to contain a 9-byte string. The array 
specification in the call may be expressed as the array name alone, 

CALL DATE (A) 

in which the first three elements of the real array A are used to hold the date 
string ; or as 

CALL DATE (A(i)) 
which causes the 9-byte string to begin at the i element of the array A. 
The date is returned as a 9-byte (9-character) string in the form 

dd-mmm-yy 
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where : 



dd is the 2-digit date. 

ramm is the 3-letter month specification. 

yy is the last two digits of the year. 

For example: 

27-JAN-74 

In the case where the array is a real array, 4-1/2 words are used to contain the 
date string; the remaining array storage is untouched. Therefore, the date string 
is stored in the first nine bytes in the elements A(i), A(i+1) , and A(i+2). The 
last three bytes of A(i+2) are untouched and should be made blank by the user if 
he intends to print the date with a 3A4 format. 

14.2.10 TIME 

The TIME subroutine allows the user to obtain the current system time or to per- 
form time conversions within a FORTRAN program (see also SECNDS subprogram, sec- 
tion 7-14.2.11). The TIME subroutine may be called by any of the following three 
statements. 

CALL TIME (A) 
CALL TIME (11,12) 
CALL TIME (A, II, 12) 

The first two forms of the TIME call return the current system time; the third is 
used to perform time conversions. 

The first type of TIME call returns the time as an 8-character ASCII string of 
the form: 

hh:mm:ss 



where : 



hh is the 2-digit hour indication, 
mm is the 2-digit minute indication, 
ss is the 2-digit second indication. 



For example: 

15:45:23 
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A 24-hour clock is used. A 2A4 format specification is generally used to output 
the time value. 

The second type of TIME call returns the time in "clock ticks" (1/60 of a second 
for 60Hz systems, 1/5J9 of a second for 5)ZfHz systems) elapsed since midnight, in 
the integer variables 1(1) and 1(2). 1(1) contains the high-order 15 bits of the 
number and 1(2) the .low-order 15 bits. Elapsed-time calculations can be performed 
as follows (the use of SECNDS is recommended, however; see Section 7-14.2.11). 



CALL TIME (11.12) 



CALL TIME (13,14) 
C COMPUTE ELAPSED TIME IN TICKS: FIRST CONVERT TICKS STORED 
C IN II AND 12 TO A SINGLE FLOATING-POINT VALUE. 

START == 11*32768.0+12 

FINISH =12*32768.0+14 
C NOTE THAT MULTIPLICATION BY REAL CONSTANT CAUSES ENTIRE 
C .CALCULATION TO BE PERFORMED IN REAL MODE. ELAPSED TIME 
C IS NOT CALCULATED IN INTEGER MODE DUE TO THE POSSIBILITY 
C OF INTEGER OVERFLOW. 

11 = ELAPSE/32768.0 

12 = AMOD (ELAPSE, 32768.0) 
WRITE(6,21j2f) 11,12 

210 FORMAT (206) 



The above example prints the elapsed time indication as follows : 

pf 263J2C 

showing the octal representation of the number of clock ticks between the two TIME 
calls. 

The third type of TIME call accepts 1(1) and 1(2) as integer values and returns 
the clock time as represented by 1(1) and 1(2) in an 8-character ASCII string in 
A. For example (the following is a continuation of the preceding example) : 



CALL TIME (A(1),I1,I2,) 
C SCALE ELAPSED TIME TO SECONDS 

ELAPSE = ELAPSE/60.0 

WRITE(6,215) ELAPSE, A(2) 
215 FORMAT (' OELAPSED TIME = I ,F5. 1, 'SECONDS '/ 
1' AS CONVERTED BY "TIME" 1 :',2A4/) 
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outputs the following: 



ELAPSED TIME =23.9 SECONDS 

AS CONVERTED BY "TIME" : 00:00:23 



14.2.11 SECNDS 1 

The SECNDS function returns the system time, minus the value of the argument, as 
a single precision floating-point value. For example, 

TIM = SECNDS (0.) 

will return the number of seconds since midnight; that is, the current time of day. 
It may be called with a non-zero argument for performing elapsed-time computations 
as in 

C START OF TIMED SEQUENCE 

Tl = SECNDS (0. ) 
C ' 

C CODE TO BE TIMED 
C 

DELTA « SECNDS (Tl) 

where DELTA will give the elapsed time. 

The value of SECNDS is accurate to the resolution of the system clock: 0.0166... 
seconds for a 6#Hz clock and 0.^2 seconds for a 50Hz clock. 

With 24 bits of precision for real values, this representation is accurate to the 
clock tick for values up to about two days in duration. 

14.2.12 Radix-50 Conversions 

14.2.12.1 IRAD5JZ( 

The IRAD5# subprogram performs conversions between Hollerith (text) strings and 
Radix-5J2f representation. (See Section 7-2.2.9 for details for the Radix-50 repre- 
sentation.) 

IRAD5^ may be called as a FUNCTION subprogram if the return value is desired, or 
as a SUBROUTINE subprogram if no return value is desired. The form of the call is 



r This function is the preferred alternative to the 2- and 3-argument forms of the 
TIME subroutine call. 
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N = IRAD5j/(ICNT, INPUT , OUTPUT) 



CALL IRAD5^(ICNT, INPUT , OUTPUT) 

where 

ICNT is the maximum number of characters to be converted (integer) . 

INPUT is an ASCII (Hollerith) text string to be converted to Radix-5^. 

OUTPUT is the location for storing the results of the conversion. 

N is the number of characters actually converted. 

Three characters of text are packed into each word of output. The number of output 
words modified is computed by the following expression (in integer mode) . 

(ICNT+2)/3 

Thus, in a count if four is specified, two words of output will be written even 
if only a 1-character input string is given as an argument. 

Scanning of input characters terminates on the*first non-Radix-50' character en- 
countered in the input string. 

14.2.12.2 RAD5J2f 

The RAD5^ function subprogram provides a simplified way of converting ASCII data 
to Radix-50 representation. The form of the call is 

X = RAD5^ (string) 

where 

string is the input ASCII string. Up to six characters are scanned and 
converted . 

X is a real variable to which the converted value is to be assigned. 
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The RAD50 function is equivalent to the following FORTRAN function. 

FUNCTION RAD5^(A) 

CALL IRAD5jZf (6,A,RAD50f) 

RETURN 

END 

14.2.12.3 R50ASC 

The R5^ASC subprogram provides decoding of Radix-5j2f encoded values into ASCII strings. 
The form of the call is 

CALL R5jZiASC(ICNT,IN,0UT) 

where 

ICNT is the number of output characters to be produced. 

IN is the variable or array containing the encoded input. Note that 
(ICNT+2)/3 words will be read for conversion. 

OUT is the variable or array into which ICNT characters (bytes) will be 
placed. 

If undefined Radix-5/J code is detected, or the Radix-Sp' word exceeds the maximum 

value 174777 , question marks will be placed in the output field, 
o 

14.2.13 SSWTCH 

The sense switch (SSWTCH) subroutine can be used in FORTRAN programs to test the 
current status of specified bits of the console switch register. The subroutine 
is called by the statement 

CALL SSWTCH (I, J) 

where: 

I is an integer value designating the console switch register bit to 
be tested; it must be within the range of & through 15. If it is 
outside of this range, no value is returned; no diagnostic is print- 
ed in this case. 

J is an integer variable that is to contain the value returned by the 
subroutine. It is set to 1 if bit I is 1 (the switch is up) or to 
2 if bit I is $ (the switch is down) . 
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14.2.14 Character String Arguments to System Subroutines 

String arguments to system subroutines must be stored in adjacent bytes. String 
values stored in integer arrays are not stored in adjacent bytes unless the /ON 
switch is used in compiling the program (see Section 7-9.2 and Chapter 7-12); such 
character strings are unacceptable as arguments to system subroutines. 

The string values returned by the DATE and TIME subroutines are stored in adjacent 
bytes independent of the use of the /ON switch. Hence, the /ON switch is needed 
to produce correct results. 

Object time format specifications (see Section 7-7. 2p) must also be stored in 
adjacent bytes. 

String arguments to the ASSIGN subroutine must be stored in adjacent bytes and have 
a terminating zero byte* The FORTRAN system automatically appends a zero byte to 
character string arguments specified as Hollerith constants. 
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PART 7 

CHAPTER 15 

FORTRAN CALLING SEQUENCE CONVENTIONS 

The calling sequence convention currently in use with FORTRAN on the PDP-11 family 
of computers is called "R5", and is used with all versions of the DOS/BATCH opera- 
ting system. 

The following points should be noted: 

1. Version V0fe (DOS/BATCH) of the FORTRAN Compiler implement subprogram 
calls $CALL, by means of an OTS threaded code routine, as described in 
Section 7-10.2. For this reason, there is no compilation-time switch 
for specifying the convention to be used. FORTRAN routines compiled by 
the Vfife Compiler may be linked with the R5 forms of the library. 

2. FORTRAN routines compiled by earlier Compilers (V4A and before) may be 
linked with the R5 form of the library. It is recommended, however, that 
these routines be recompiled using the Vj0"6 Compiler. 

15.1 PDP-11 FORTRAN "R5" CALLING SEQUENCE CONVENTION 

The form of the subprogram call under the R5 convention is the same as that used 
with earlier versions of FORTRAN under the QOS/BATCH operating system. 

15.1.1 The Call Site 

The basic form of the call follows. 



JSR 


%5,SUB 


BR 


NEXT 


.WORD 


ADR1 



;CALL SUBPROGRAM 

; BRANCH OVER ARGUMENT LIST 

; FIRST ARGUMENT ADDRESS 



.WORD 



NEXT: 



ADRN 



;NTH ARGUMENT ADDRESS 

;NEXT EXECUTABLE INSTRUCTION 



Note that the low-order byte of the branch instruction contains the count of the 
'number of arguments that follow. A maximum of 127 arguments are permitted. 

15.1.2 Return 

Control is returned from the called program unit to the calling program unit by 
restoring the stack pointer register (SP, register 6) to its value at the time of 
entry, if necessary, and executing the following. 
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RTS 5% 



15.1.3 Return Value Transmission 



FUNCTION subprograms return values in general registers R$ through R3. The number 
of registers used is determined by the data type of the function's returned value, 
as follows. 

Data Type Value Returned in 



BYTE ( LOGICAL* 1) 

LOGICAL 

INTEGER 



} ¥ 



REAL R0, Rl 

DOUBLE PRECISION 

COMPLEX Rgf, Rl , R2 , R3 



Example : 



The FORTRAN statement 



CALL SUB (A1,A2,A3) 
results in a subprogram call equivalent to the following. 



.GLOBAL SUB 

JSR %5 , SUB 

BR $Fnnnn 

.WORD Al 

.WORD A2 

.WORD A3 
$Fnnnn: 



Within the subprogram, arguments can be accessed by means of indexed references, 
using register 5 as the index register. In particular, the instruction 

MOVB §%5, ... 
obtains the number of arguments in the list, 

MOV 2%5, ... 
obtains the address of the first argument, and, in general, the instruction 



7-162 



MOV 2*n(%5), ... 
obtains the address of the n argument. 
15.1.4 Null Arguments 

Null arguments are represented in an argument list by using an address of -1 

(177777 ) . This address can be easily detected and generally assures that an error 
o 

will result at execution time if a null argument is passed to a subprogram that is 
not prepared to accept it. Null arguments are included in the argument count as 
shown in the following examples: 

FORTRAN Statement Resulting Argument List 

CALL SUB 
CALL SUB ( ) 

CALL SUB (A,) 



CALL SUB (,B) 



15.2 MACROS FOR PDP-11 FORTRAN CALLING SEQUENCE 

15.2.1 Introduction 

Three macros have been defined to facilitate the writing of PDP-11 subprograms 
that interface to PDP-11 FORTRAN compiled programs. These provide three functions. 

1) Calling a FORTRAN subprogram, including construction of the appropriate 
argument- list. 

2) Returning to a FORTRAN subprogram, including the transmission of the 
return value (if any) . 

3) Obtaining the value returned by a FORTRAN subprogram and moving it to the 
desired destination. 

The goal of these macros is to. make assembly language routines independent of the 
details of the FORTRAN calling sequence conventions. 

The general form of each macro call is 



BR 


.+2 


BR 


.+4 


.WORD 


-1 


BR 


.+6 


.WORD 


A 


.WORD 


-1 


BR 


.+6 


.WORD 


-1 


.WORD 


B 
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F4CALL Subroutine-name, argument-list, optional-label 
F4RTN Typecode , location 
F4VAL2 Typecode , location . 

These macros are discussed in detail in the following sections. 

In the following discussions reference will be made to the addressing mode of a 
macro argument. If unfamilar with the details of the PDP-11 addressing modes 
please refer to the appropriate processor handbook. Briefly, the following should 
be noted. 



1) In some cases the code generated by these macros is dependent on the 
addressing mode of one or more arguments. The addressing mode is obtained 
by using the .NTYPE directive described in Part 6, MACRO. 

2) Intuitively, the addressing mode is the six bit octal value that would be 
placed in either the source or destination field of a PDP-11 instruction. 



3) The following cases are particularly important. 



Argument For m 

Simple identifier, 
e.g., X 

Push to stack, 
e.g., -(%6) 



Addressing Mode 
67 

46 



15.2.2 Calling a FORTRAN subprogram 



The form of the macro call is 



F4CALL 



NAME , <ARGLIST> , LABEL 



where 



NAME is the name of the subprogram to be called. 

ARGLIST is a list of addresses of the arguments. 

LABEL is an optional identifier to be placed oh the argument list. 
In most cases this label is not needed. 



Examples : 



F4CALL SUB,<A,B,C> 

F4CALL X 

F4CALL X,<A,>,L1 



;CALL SUBROUTINE SUB 
;WITH ARGUMENTS A, B, AND C 

;CALL SUBROUTINE X WITH NO ARGUMENTS 
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There are two parts to calling a subroutine: filling in the argument list (either 
at assembly time or at execution time) and transferring control to the subroutine. 
The latter is performed by 

JSR REG, NAME 

where the REGister used depends on the convention. 

Filling in the argument list is slightly more complex. The argument list itself 
consists, of a block of N+l words where the low byte of the first word contains the 
number of arguments, N, and the next N words are the addresses of the arguments. 
There are four cases to be considered. 

Case 1 - assembly time constants 

In many cases the argument is at a fixed location and the address can be simply 
assembled into the argument list. This is done for all arguments with an addressing 
mode of 67. 

Case 2 - argument value on stack 

Often a result has been computed and left on a stack. In this case the address is 
the sum of the contents of the stack register (which need not be %6) and a constant 
offset. This is represented by an argument of the following form. 

Register , of f set 

For example: 

<%6,0> 
<%1,6> 

For this situation the address of the data is computed by the sequence 

MOV Register, Argument List Entry 

ADD Offset, Argument List Entry 

Note that the address is stored in the appropriate position in the argument list; 
if the offset is zero, no ADD instruction is needed or produced. 
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Case 3 - empty argument 

In rare cases the user may want to fill in an argument address by code unrelated 
to the code that the P4CALL macro produces. Here it is necessary to designate 
the argument list label as part of the F4CALL. (Note that the first argument 
address goes at LABEL+2, the nth at LABEL+2*n / and so on.) Space for such an 
argument is reserved in the argument list and an aiddress of -1 is assembled in. 

Case 4 - general argument expression 

If none of the above cases holds, then the argument address is obtained by 
generating the following instruction. 

MOV Argument , Argument List Entry 

Note that any auto-incrementing or auto -decrementing indicated in a general argu- 
ment expression will be performed as a result of the MOV instruction. 

Putting these all together we illustrate with the example shown in Figure 7-14. 
Figure 7-15 illustrates how the F4CALL macro relates to the FORTRAN source code. 
Assume the following: 



ARG1 is the address of argument 1,, 

R0 contains the address of argument 2, 

Rl contains a pointer to a two word table which contains the 

addresses of arguments 3 and 6, 
R5 points to the argument list of the call that entered this 

subroutine and the first argument of that call is to be the 

4th argument of this call, and 

the address of the fifth argument is the third word on the 

processor stack. 

the seventh argument is filled in elsewhere. 

Then write 

F4CALL NAME, <ARG1„R0, (Rl)+,2 (R5) , <SP,6>, (R1)+,>,L 

which would generate the following code for the PC calling sequence con- 
vention : 

MOV R^,L+4 

MOV (Rl)+,L+6 

MOV 2(R5),I/flrf 

MOV SP,L+12 

ADD #6,L+12 

MOV (R1)+,L+14 

MOV R5,-(SP) 

MOV #L,R5 

JSR PC, NAME 

BR $KPLST 



Figure 7-14 
Argument List Construction 
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L: 


.BYTE 


id 




.WORD 


ARG1 




.WORD 


i 




.WORD 







.WORD 







.WORD 







.WORD 







.WORD 


-1 


$KPLST: 


MOV 


(SP)+,R5 


(Where $KPLST is a 


macro generated 


label.) 



Figure 7*-14 (Cont.) 
Argument List Construction 





SUBROUTINE A(X, 
DIMENSION Z(10) 


Y,Z) 




CALL 


B(I+J,X+C, 


D,Z,(I) ,X) 




END 






might generate the 


following 


code using these macros: 


Compute 
Compute 
Compute 


I + J and place 
X + C and place 
address of Z(I) 


on stack (1 word) 
on stack (2 words) 
and leave in RO 


F4CALL B,«SP 


,4>,<SP,£k>, 


D,R#,2(R5)» 



Figure 7-15 
Relationship Between F4CALL and FORTRAN Source Code 



15.2.3 Returning to a FORTRAN Program 



The form of the macro call is 



F4RTN 



Typecode , Location 



where 



Typecode is a single letter representing the data type as follows. 



B - Byte 

L - Logical 

I - Integer 

J - Double Integer 

R - Real 

D - Double Precision Real 

C - Complex 
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Location is either a simple variable name (addressing mode 67) or a stack pop 
addressing mode (2X) . 

i 

The macro generates the move instruction needed to place the returned value in the 
correct position for recovery by the calling routine. (This is presently in the 
general registers %#-%3 but we wish to retain the option of using the floating 
registers on 11/45 FPP systems at a future date.) Finally the appropriate RTS 
instruction is generated. 

If the macro is written with no arguments then only the RTS instruction is produced. 

Note that the user code is responsible for assuring that the stack pointer (%6) is 
correctly positioned at runtime. 

Several examples are shown in Figure 7-16. 







• 




Example 1 




Macro call: 


F4RTN 




Expanded code- 


RTS 
Example 2 


R5 


Macro call: 


F4RTN 


B, (%1) + 


Expanded code- 


MOVB 


(%1)+,R0 




RTS 


R5 




Example 3 




Macro call: 


F4RTN 


R,L0C+4 


Expanded Code- 


MOV 


LOC+4,R,0' 




MOV 


LOC+4+2,Rl 




RTS 


R5 



Figure 7-16 
Returning to the FORTRAN Program 



15.2.4 Obtaining the Returned Value 



The form of the macro call is 



F4VAL2 Typecode , Location 
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where 



Typecode is a single letter representing the data type as in the F4RTN macro, and 
location is either a simple variable name (mode 67) or a stack push addressing mode 
<4X). 

The expansion produced moves the returned value into the designated location or 
onto the designated stack. 



Examples are shown in Figure 7-17. 



Note that code generated for this macro is not currently call sequence dependent. 
This, macro is included to isolate the entire call sequence mechanism so that, for 
example, floating-point values may some day be returned in FPP registers. 





Example 1 




Macro calls 


F4VAL2 


I,-(SP) 


Generated code- 


MOV 


R0,-(S) 


» 


Example 2 




Macro call: 


F4VAL2 


C,VAR 


Generated code- 


MOV 


R^,VAR 




MOV 


Rl,VAR+2 




MOV 


R2,VAR+4 




MOV 


R3,VAR+6 



Figure 7-17 
Return Value Transmission 

15.2.5 Obtaining and Using the Macros 

The macros are distributed as part of the system macro file SYSMAC.SML. They are 
defined in the user module in the following way. 

The macro .F4DEF must be obtained from the system macro file via a 

.MCALL .F4DEF 

directive. Then the macro .F4DEF must be invoked with an argument specifying the 
calling form to be defined: tf f or the R5 form. Its action is simply to define 
the macros as appropriate for the value of its argument. 
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Thus to obtain the R5 form call the user includes the following. 

.MCALL .F4DEF 
.F4DEF 

As a side effect .F4DEF also defines the variable .F4SEQ to have the value of its 
argument. This symbol may then be used to conditionally assemble sections of code 
which, for reasons not anticipated here, may be call sequence dependent. 

Finally the .F4DEF macro may be called without an argument in which case the value 
of »F4SEQ is assumed to be already defined as 0, and its value is used to define 
the remaining macros. This is useful, for example, where the value of .F4SEQ is 
defined in a separate parameterization file used with a collection of routines. 

15.2.6 Programming Cautions and Notes 

In addition to the macros F4CALL, F4RTN and F4VAL2, several additional macros are 
defined for use by the user level macros and several variables are used for com- 
munication between macros at assembly time. All of these internal macro and 
variable names begin with the period character. Users are urged to not use iden- 
tifiers beginning with the period character. 

It is a general policy that use of identifiers containing the period character is 
to be specified exclusively by DEC. Users are urged to avoid all such names in 
order to assure that name conflicts will not be introduced in future releases of 
existing or new software. 

In the R5 form of calling convention FORTRAN register 5 always points to the argu- 
ment list of the routine calling the current one and this pointer is conveniently 
saved and restored by the JSR R5/RTS R5 pair. Fortunately the new FORTRAN genera- 
ted code will not require that the argument pointer be also available in register 
5 and so does not need to explicitly save/restore R5. 

For those concerned about minimizing core and willing to expend more care in check- 
ing the compatibility, it may be worthwhile to try the following call. 

.F4DEF 2 

This is identical to .F4DEF 1 except the explicit save/restore of R5 is not 
generated . 
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PART 7 

CHAPTER 16 

FORTRAN TRACE PACKAGE 

16.1 TRACE PACKAGE 1 

A FORTRAN trace package is provided as a library-called TRCLIB.OBJ that provides a 
FORTRAN-level trace of program execution. Actions that can be traced include 
assignment to variables, call and return from subprograms, and transfers of control 
within a program unit. 

The library consists of two modules: TRACEX, which is a special version of selected 
OTS threaded code service routines, and TRACEF, which is a DEC-supplied FORTRAN sub- 
routine for performing formatted tracing of program action. If desired, a user 
may supply his own version of the SUBROUTINE TRACEF for altering the trace dump to 
his own needs. 

16.2 TRACE OUTPUT DESCRIPTION 

A summary of the types of trace provided and the trace output generated is shown 
in Table 7-20. 

16.3 SELECTIVE CONTROL OF PROGRAM TRACING 

Trace output for each type of action may be independently 'enabled or disabled 
either under control of the PDP-11 console switches or under program control. 

v 
The normal or default mode is to control tracing using the console switches. In 
this mode each time the trace package is entered the console switches are inter- 
rogated to determine whether a trace output is to be produced. Each type of 
trace is output only if the corresponding console switch is in the up position. 
For example, subprogram entry information is printed only if console switch 1 is 
up, subprogram return only if switch 2 is up, and so on. 

Tracing may be placed under program control by means of the FORTRAN subroutine 
TRCTRL. This subprogram is called as follows. 

CALL TRCTRL (ITYPE,ICODE) 



J This package depends on the implementation of FORTRAN compiled output using 
threaded code techniques such as are found in V#6. It may not be compatible 
with later versions. 
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Table 7-20 
Trace Output Description 



Trace Type 


Description 


Trace Output 


1 


Entry to FORTRAN 


ENTER name 




subprogram 




2 


Return to calling 
program 


RETURN TO name 


3 


Sequence number of 
statement, about to 
be executed 


SEQ . nnnn 


4 


Transfer of control 
resulting from any 
type of GOTO statement 


GOTO nnnn 


5 


Arithmetic assignment 


type = value 


6 


Arithmetic assignment 
to variables specified 
in CALL TRCLST state- 
ment 


type = value 


7 


ASSIGN statement 


ASSIGN = nnnn 


where 






name 


is a subroutine or function 


name 


nnnn 


is a statement sequence number 


type 


is the type of value being 
for arithmetic assignment 


assigned: INTEGER, REAL, etc. 


value 


is the value assigned, printed in the appropriate format 



where 



ITYPE 



ICODE 



is an integer value designating the type of trace output 
to be effected (see Table 7-20) . 

is an integer value designating whether tracing is to be 
enabled or disabled (l=enable, 0=disable) . 



Special cases of the above: 

CALL TRCTRL (6,0) places tracing under program control with all tracing 

disabled. 

CALL TRCTRL (p,l) places tracing under program control with all tracing, 

enabled . 
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CALL TRCTRL (-1,^) places tracing under control of the cfonsole 

switches. 

Note that more than one type of trace may be enabled by multiple calls to TRCTRL. 
For example, to trace only subprogram calls and returns, the two statements 

CALL TRCTRL (1,1) 

CALL TRCTRL (2,1) . 4 

may be executed. An initialization call of CALL TRCTRL (^,/) is unnecessary be- 
cause any call to TRCTRL that changes the tracing mode from switch control to 
program control automatically disables all tracing except that specified by the 
call. 

In many cases it is desirable to trace assignment statements for a designated set 
of variables only. This can be accomplished by the TRCLST entry. The form of 
this call is 

CALL TRCLST (V. , V_ , . . . , V ) 
12 n 

where each V is the name of a variable to be traced. (Note that only single 
elements can be specified. An unsubscripted array name causes only the first ele- 
ment of the array to be traced, not the entire array.) Up to 2$ variable names 
may be specified for tracing in the program (regardless of the number of TRCLST 
calls) ; variables entered in excess of 2jef are simply ignored. Variables may be 
deleted from the trace list by executing the following statement. 

i 

CALL TRCDEL (V , ,V„,...,V ) 
12 n 

An attempt to delete from the trace list a variable that does not exist in the 
list causes no problems. 

16.4 TRACE OUTPUT DEPENDENCE ON COMPILATION OPTIONS . 

If a FORTRAN program or subprogram is compiled with the /SU option in effect, 
sequence number information is not available in the compiled output.- In this 
case all trace output requiring a sequence number prints a / for that number. 
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If a FORTRAN program or subprogram is compiled with the /OP option set to a value 
other than $, many assignment operations will be undetectable by the trace package. 
To ensure that all assignment actions are traceable, use /GPij/in the compiler 
control string. 

16.5 USAGE 

No modifications to FORTRAN source programs are required to obtain tracing under 
switch control. If control of tracing through TRCTRL, TRCLST, or TRCDEL is de- 
sired, however, appropriate calls to those entry points must be inserted in the 
source program where desired. 

To provide the tracing features described above, the library TRCLIB is included 
in the input file portion of the LINK command string following all other user.'s 
files, but before the normal library, FTNLIB. 

Example : 

A complete example of the commands needed to trace a demonstration program, 
DEMO.FTN, and the output provided by the TRACEF subroutine, is presented 
in Figures 7-18, 7-19 and 7-20. 



*,RUN FORTRAN 

FQr?T&ANi V06.12 

#OD1O,L,PKDEMO/QPJ0 

HfC 

'.KILL 

'.RUN LINK 

LI^K Vfeil 

«PFM0<Q£MQ/CC,TRCL1B/L,FTNL1B/L/E 

tf tC 

".KILL 

', ASSIGN LP I i 6 

'.run demo 



Figure 7-18 
Console Commands to trace DEMO.FTN on the Line Printer 
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0001 ' 




0002 




0003 




0004 




0005 




0006 






C 




C 




c 


0007 






c 




c 




c 


0008 






c 




C 




c 


0009 




0012! 






C 




r 




r, 


0011 




0012 


10 




c 




c 




b 


0013 




0014 




0015 




0016 




0017 




001.8 






c 




c 







0019 




0020 




0021 




0022 




0023 




0024 






c 




e* 




c 


0025 






c 




c 




c 


0026 




0027 


40 


0028 


50 


0029 


60 



byte: 8 

LOGICAL L 

INTEGER I 

REAL R 

pOUBtE PRECISION D 

complex c 
trace: eveRvTming 

CALL TRCTRL(0il) 

TyPES 1, 2i AND 3} CALL, RETURN AND SEQUENCE 

CALL OUMMy 

TURN OFF SEQUENCE TRACE 

CALL TRCTRL(3,0) 
CALL DUMMY 

T Y PE 4i GOTO TRACE 

GOTO 10 
CONTINUE 

TyPE 5.1 ARITHMETIC ASSIGNMENT TRACE 

R a 1 

L « ,TRUE, 

I ■ 2 

R s 1,5 

9 2,3D0 

C a (1,1^2,3) 

Ty'PE 6r SELECTED VARIABLE TRACE 

CALL TRCTRL(S,0) 

CALL TRCLST<Y,8) 

Y c 2 

CALL T'RCDEU.(Y,I) 

Y b 3, 

T Y PE 7i ASSIGN STATEMENT TRACE 

ASSIGN 40 TO K 
MORE TyPE 41 GOTO TRACE 

GOTO K 

M ■ 2 

GOfO (50,60) M 

CONTINUE 



0030 CALL.E X IT 

0031 EnO 



Figure 7-19 

Compilation Listing of DEMO. FIN 
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FORTRAN V06.12 1 & * 1 7 ; 3 



PQUTlNfS CALLED! 

?RCT*L, DUMMY , TRCLSTi TRCOEL, EXIT 

OPTIONS a /ONi/QPI0 

block length 

MAIN, 233 (000676)* 
»*rO M P!LrR ----- qORe** 

PHASE USED FREE 

HECLARA.TIVES 03216 01576 

FIXECUTABLES 03438 01334 

ASSEMBLY 03293 04074 



FORTRAN V06.12 15>»17{51 

0001 SUBROUTINE DUMMY 

0002 RETURN 

0003 END 

OPTIONS = /ONi/OP|0 

PLOCK LENGTH 
HUMMY 10 (020024)* 

**r; MpiLCR ----«? CqrE*« 

PHAS£ USED FREE 

DECLARATIVES 00297 01495 

FXFCUTA3LES 00297 01499 

ASSEMBLY 071067 04290 



Figure 7-19 (Cont.) 



7-176 



enter main', 

SrQ, 7 
sEq". a 

ENTER DUMMY 
SEQ, ?. 
RETURN TO MAIN'. 

SEQ, 9 

ENTER DUMMY 
RETURN TO MAIN*. 

GOTO SeQ. 12 

BYTE » 1 

LOGICAL « 177777 

INTEGER s 2 

REAL s 1,5000003 

DOUBLE a 0..2300000000D 01 

COMPLEX ■■ ( l',1000 i*2,3000 ) 

REAL «. 2,0000000 

ASSIGN a 27 

GOTO SEQ, 27 

GOTO SEQ', 29 



Figure 7-20 
Trace Output of Program DEMO 
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PART 7 

CHAPTER 17 

FORTRAN DEVICE TABLE LISTING 



.TITLE SDV8?9 

GL08L uurvTB 

'.CSECT 

S0E VT B V 009A 



J c 0P yr1GHt 1971,1972 digital e qu IpHENt c r PO rAtJqN# mayna r d,ma S s 
j 

ItHEqE Apt THE FORTRAN DEylCE TABLE ENtrIEs 

IWITH THE DEVICE TABLE HEADER AND ENTRY VECTOR 



1AD0R OF ENTRY FOR ERR LOG DEVICE 
jERROR LOG DEVICE IS -3<DEVM3> 
IAODR OF ENTRY FOR ERR MSG ( FILE 

JNUMBER OF ENTRIES IN ENTRY VECTOR 

JDEVICE NUM OF ERROR LOGGING DEVICE 

jDEVICE NUM OF ERROR LOGGING DEVICE 



jADDr qF DEylCE 1 ENTRY 

jADDR OF DEVICE '2 ENTRY 

JADDR OF DEVICE 3 ENTRY 

lADDR OF DEVICE 4 ENTRY 

jADDR OF DEVICE 5 ENTRY 

jADDR OF DEVICE h ENTRY 

lADDR OF DEVICE 1 ENTRY 

jADDR OF DEVICE 9 ENTRY 





', IFNDF 


s$X 
DEVM: 




•JORO 




.WORD 


DEVE 1 




,ENDC 




IDEVTBl .WORD 


s, 




.IFNOT 


RSX 




.WORD 


«3 




.ENDC 






". IFDF 


RSX 




', WORn 


5, 


■ 


.ENDC 




J THE 

8 


DEylCE TAP 


LE EN 


t 


.WORD 


DEyl 




'.WORn 


DEV2 




.wnftn 


HEV3 




,WORn 


DEV4 




", wnRn 


r»£V5 




'.WORD 


DEV6 




'.WORD 


DEV7 




.WORD 


0EV8 



J ENTRY 1 q f DEVICE TABLE 



t 

DEVU 



.WORH 

'. 1FNDE 

.RAD50 

'.ENDC 

.IFDF 

".RAOfB 

'.ENDC 

.3YTE 



*SX 
/SY / 

RSX 
/KO/ 



iLINK BLqCK PTR 

{PHYSICAL DEVICE NAME DEFAULT 

I HOW OPEN SWITCH 



.B^TE a 

',rAQ5>pi /F r/ 

'.RA-DSg /mi/ 

'.RADRa /DAT/ 



i UNIT NUM DEFAULT 
jDEFAuLT FILE NAME 

IDEFAULT EXTENSION 
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",3yTE 


2 


',8YTF 





'.BYTE! 


% 


,8YTF 





,WORn 





".WORD 





, WORD 





.WfJRn. 





,woRn 





', worh 





",W0Rn 





, wo'Rn 






NO AUTO DEL, GROUP & OTHER READ/RUN ONLy 

device; status switch 

mode of i/o n funcn word (random) 

status OF I/O 

RECORD CjOUNT * BLOCK HUM (RANDOM) 
BUFF ADDR (RANDOM) 



BUF L'EN (RANDOM) 

ASS06IATED VAR ADDR 

NUM RECORDS IN FILE 

RECORD LENGTH (FROM 

USER ID CODE 

ERROR VAR ADDR (FROM SETFIL) 



(FROM DEFINE FILE) 
(FROM DFFINE FILE) 
DEFINE FILE) 



ENTRY 2 F OEylCE 'tA*bLE 



DEv^t 



WORO 


'■{ 






IFNDF 


RSX 






RAD50 


/SY / 






ENDC 








IFDF 


QSX 






RAD90 


/Tf/ 






E^OC 






■ ■ 


8YTE 


diti 






RAD50 


/FOR/ 






RAQ t 50 


mzt 






RAD«50 


/DAT/ 






BYTE 


333,0, 





,0 


WORD 


0i0#3i 


21 


0,0,8, 



ENTRY 3 o F DEVICE TABLE 



D E y 3 J 



WORD 


■'<* ■ 


IFNHF 


RSX 


RAO90 


/SY / 


ENQC 




IFDF 


8SX 


RADS0 


/pp/ 


ENDC 




8YTF 


010 


RAD50 


/FOR/ 


RAD^S 


Z00.3/ 


RAD3J? 


/DAT/ 


BYTE 


533,04010 


WORn . 


a # , 2> 1 » , » > 



ENTRY 4 F DEVICE T ABLE 
DEv^J 



'. woro 

.RAD90 
• gYTr 
,«Aq5b 
". RA0«5«3 
,RAD*8 




/PP/ 

0,0 

/FOR/ 
/004/ 
/DAT/ 
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,8yT£ 233,3,0,0 

.word a, 0,0,0, 0,0, 0,0 



IEN TR Y 5 r DEVICE: T ABLE 



DE V 5J 



•WqrB 
.RAD50 

'.byte: 

.RAD50 

,RAD'50 
'.RAD^E 

.BYTE 

,WORn 





/LP / 

/FOR/ 
/005/ 
/DAT/ 
233,9, 



0,0 



0,0,0,0,0,0,0,0 



EN TRY 5 r nE V 1CE T ABLE 



DEy5l 



".WOR D 





'.RAD50 


/KB / 


'.byte: 


0,0 


'.RAD50 


/FOR/ 


'.RAD50 


/006/ 


'.RAD50 


/DAT/ 


'.3YTE 


233,0,0,0 


'.WORD 


1,0,0,0,0,0,0,0 



ENTRY 7 F D-EylCE T a BLE 



DEy7! 



WOR D 


$ 


IFNOF 


RSX 


RAD50 


/SY / 


ENDO 




IFDF 


c?SX 


RAQ50 


./UO/ 


EMDC 




BYTE 


0, 


RAD50 


/FOR/ 


RAQ30 


/007/ 


RAQ50 


/DAT/ 


BYTE 


233,0,0,0 


WORD 


0,0,0,0,0,0,0,0 



JENTRY 8 qF DEVICE t a BLX 

J 

I 

DEv3l 



• WoR? 

', IFNDF 

'.RAD52 

»E N DC 

•IFDF 

,RA D 50 

.ENDC 



RSX 
/Bl - 

/A D / 
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'.BYTE 

'.RAD50 

'.RAD58 

',RAD^0 

.BYTE' 

'.WORD 

', IFNOF 



/FOR/ 
/008/ 
/DAT/ 

233,0, 
Oi 0t 3, 
RSX 



0i0 
0,0 



0,0,0 



J 

JENt rY m 
{SPECIAL 
{LOGICAL 
J 

DE V M35 



3 n r OEylCE T AB LE< E RR 
ENTRY USED AS ~" ' ' 
UNIT NUMBER = 



LqGGING DEwlCE) 
DO.Seli "ERROR LOGGING DEVICE 



LOGICAL 



DEVICE NAME * CMO 



,RAD9g 
.BYTE 

,RAD?0 

',RAQ5g 
'.RAQ50 
'.BYTE 
.WORD 



/KB / 

/FOR/ 
/CMO/ 
/DAT/ 

233,0, 



{LINK BL n CK , 
{DEFAULT PHY! 



lAL DEVICE NAME 



{DEFAULT FfLE NAME 
}FQRCMO,DAT 




0,0,0 r0 



SpECIAL ENTRY /oR Ef? ROR pRQ^SSORS M S& PIUE 





• WOR D 


c> 


DEVERRl 


.WORD 







'.RAD50 


/ERR/ 




'.BYTE 


1 




'.BYTE 


tl ' 




'.RAD50 


/SY / 




'.WORD 


v*< 




'.BYTE 


4 




.BYTE 






'.RAD50 


/FOR/ 




'.RADSg 


/RUN/ 




',RAD^0 


/DGN/ 




'.BYTE 


1 




'.BYTE 


1 




.BYTE 


322,0 




.WORD 


4 




.WORD 


£ 




'.WORD 





{ 


'.WORD 


;-) 


{ 


.ENDC 
'.END 





SET NAME 

OS NAME FOLLOWS 



{LINK BLqCK E rr RT N ADDr 
jLINK PTR 
jLOG DATA 

{PHYSICAL 

{UNIT NUM 

{PHYSICAL DEVICE NAME 

{DEFAULT TO SYSTEM DEVICE 

{FILE BLOCK ERROR RETURN ADDR 

{HOW TO OPEN (OPEN!) 

{ERROR RTN CODE 

{FILE NAME 



{USER ID CODE 

{ALLOW ONLY INPUT ACCESS 
{FUNCTION WORD (READ) 
{BLOCK NUM 
}8L0~CK ADDR 
{BLOCK LENGTH 
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